|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 281, Issue 14, 9552-9559, April 7, 2006
Identification of a Novel Region of the Cytoplasmic Dynein Intermediate Chain Important for Dimerization in the Absence of the Light Chains*![]() ![]() 1
From the
Received for publication, October 31, 2005 , and in revised form, January 17, 2006.
Cytoplasmic dynein is the multisubunit protein complex responsible for many microtubule-based intracellular movements. Its cargo binding domain consists of dimers of five subunits: the intermediate chains, the light intermediate chains, and the Tctex1, Roadblock, and LC8 light chains. The intermediate chains have a key role in the dynein complex. They bind the three light chains and the heavy chains, which contain the motor domains, but little is known about how the two intermediate chains interact. There are six intermediate chain isoforms, and it has been hypothesized that different isoforms may regulate specific dynein functions. However, there are little data on the potential combinations of the intermediate chain isoforms in the dynein complexes. We used co-immunoprecipitation analyses to demonstrate that all combinations of homo- and heterodimers of the six intermediate chains are possible. Therefore the formation of dynein complexes with different combinations of isoforms is not limited by interaction between the various intermediate chains. We further sought to identify the domain necessary for the dimerization of the intermediate chains. Analysis of a series of truncation and deletion mutants showed that a 61-amino-acid region is necessary for dimerization of the intermediate chain. This region does not include the N-terminal coiled-coil, the C-terminal WD repeat domain, or the three different binding sites for the Tctex1, LC8, and Roadblock light chains. Analytical gel filtration and covalent cross-linking of purified recombinant polypeptides further demonstrated that the intermediate chains can dimerize in vitro in the absence of the light chains.
Cytoplasmic dynein 12 is a motor complex responsible for the transport of membranous vesicles and different cargo proteins toward the minus ends of microtubules including mitosis, nuclear migration, Golgi and centrosome localization, organelle and viral transport, and axonal transport (1-4). It is a large multisubunit assembly ( 1.5 MDa) containing two copies each of DYNC1H1 (the heavy chain), DYNC1I (the intermediate chains), DYNC1LI (the light intermediate chains), and two copies of each of three distinct light chain families (1, 5, 6). The globular heads of the cytoplasmic dynein heavy chains generate the force for microtubule-based motility (7-9). The heavy chains dimerize via their N-terminal stalks, and the stalks contain both the light intermediate chain and the intermediate chain binding regions (10, 11).
The intermediate chains, light intermediate chains, and three pairs of light chains are at the base of the cytoplasmic dynein 1 complex and form its cargo binding domain (2, 12-15). The two intermediate chains are essential components of the cytoplasmic dynein complex, serving as platforms that interact with the heavy chains, the three light chain dimers, as well as the p150 subunit of dynactin, a cargo adaptor complex, and other protein cargoes (12, 13, 16, 17). The C-terminal half of the protein contains seven WD repeats (18, 19), which are presumed to form a Biochemical analyses showed that cytoplasmic dynein 1 can be fractionated into two complexes, one containing the heavy chains and the light intermediate chains and the second composed of the intermediate chains and three light chain families (25-27), but it is not known how the intermediate chains might interact directly or whether their interaction requires one or more of the light chain families. A role for the LC8 light chain as a molecular glue has been suggested by its presence in complexes unrelated to dynein (3, 28). Supporting this hypothesis are the observations that an N-terminal fragment of the intermediate chain becomes more ordered upon binding of the LC8 or Tctex1 light chains (29, 30). In addition, there are two regions of alternative splicing in the N terminus of the two intermediate chain genes that generate at least six unique isoforms (1, 12, 21, 31, 32). Although the functional significance of these intermediate chain isoforms is unknown, there are data supporting the hypothesis that the isoforms of the light intermediate chains and isoforms of the light chain families contribute to the specificity of dynein cargo binding, reviewed in Refs. 1, 2, and 33. It is known that the two light intermediate chain family members form only homo-oligomers producing two distinct cytoplasmic dynein complexes in cells (11). It has been suggested that the two Tctex1 family members also only form homodimers (34). However, homo- and heterodimerization has been observed between the two members of the Roadblock light chain family (35, 36). Thus, it is of great interest to determine whether different intermediate chain isoforms are able to associate with each other. In this report, we investigated the functional basis for the interaction of the intermediate chains. We found that the six intermediate chain isoforms can interact in all combinations of homo- and heterodimers when expressed in cultured cells. In addition, the region of the intermediate chain necessary for dimerization is shown to be a previously uncharacterized region of 61 amino acids. This region is immediately N-terminal to the Roadblock binding region, and it does not overlap with the previously identified intermediate chain subdomains: the N-terminal coiled-coil domain, the C-terminal WD repeat domain, and the three distinct light chain binding domains. Furthermore, in vitro analysis of purified intermediate chain fragments confirmed that none of the light chains are necessary for intermediate chain dimerization.
PlasmidsMammalian expression constructs of the intermediate chain isoforms (intermediate chain 1A, 1B, 1C, 2A, 2B, and 2C) were constructed using standard PCR and DNA ligation methods. Briefly, the full-length intermediate chain 1 isoforms were inserted into the XhoI/NotI sites of both pCMV-HA3 and pCMV-Myc vectors (Clontech), whereas the full-length intermediate chain 2 isoforms were inserted into the KpnI/NotI sites of the pCMV-HA and pCMV-Myc vectors. Truncation and deletion constructs of intermediate chain 2C were also cloned into both pCMV-HA and pCMV-Myc vectors using standard recombinant DNA methods. Hexahistidine-tagged versions of the intermediate chain 2C truncations 1-150 and 151-250 were generated by cloning the coding region into the pET14b expression vector (Qiagen Inc.). The coding region for mouse roadblock-1 was obtained from expressed sequence tag using the PCR. A GST-tagged version of roadblock was generated by cloning the coding region into the BamHI/XhoI sites of the pGEX4T-1 expression vector (Amersham Biosciences). All constructs were sequenced to confirm gene sequence and correct reading frame.
Expression and Purification of Fusion ProteinsBacterial expression of recombinant proteins was carried out using Escherichia coli BL21 (DE3) as the host cells. To express GST-Roadblock-1, or GST, host cells containing the expression plasmid were grown in LB medium at 37 °C to an A600 of
GST Pull-down ExperimentsGST-Roadblock-1 or GST alone ( Polyacrylamide Gel Electrophoresis and ImmunoblottingProtein samples were resolved by SDS-polyacrylamide gel electrophoresis on 12.5% polyacrylamide-Tris gels and transferred electrophoretically onto polyvinylidene difluoride membranes (Bio-Rad). Membranes were blocked 1 h at room temperature in blocking solution containing 10% (w/v) skim milk powder and 0.1% Tween 20 in Tris-buffered saline. Primary antibodies were diluted in blocking solution containing 5% bovine serum albumin in Tris-buffered saline and incubated with the membrane for 1-2 h at room temperature. Membranes were washed in Tris-buffered saline and then incubated with horseradish peroxidase-conjugated secondary antibodies for a further hour at room temperature. Labeled bands were visualized using enhanced chemiluminescence (ECL, Pierce) according to the manufacturer's instructions.
Cell Culture, Transfection, and Co-immunoprecipitations293T cells were cultured in Dulbecco's modified essential medium containing 10% calf serum (HyClone) at 37 °C in a humidified, 5% CO2 incubator. For transient transfection, cells were plated on a 6-cm Petri dish to Covalent Chemical Cross-linking ExperimentsCross-linking of the bacterially expressed intermediate chain 2C truncations were performed in 20 mM Hepes, pH 8.0, 150 mM NaCl, 5 mM EDTA and incubated with 0.1 mM disuccinimidyl glutarate (DSG, Pierce) for 30 min at room temperature. Reactions were terminated by the addition of gel sample buffer, and the samples were resolved by SDS-PAGE. DSG was dissolved in N,N-dimethylformamide, and the final N,N-dimethylformamide concentration in the cross-linking reaction mixture was 5% (v/v). Cross-linking of the mammalian expressed full-length 2C was performed in a similar manner except that the cross-linking buffer was 20 mM Hepes, pH 8.0, 150 mM NaCl, 1 mM EDTA, 1% Triton X-100, with protease inhibitors.
Analytical Gel Filtration ChromatographyThe homogeneity of the complexes formed by the truncated forms of intermediate chain 2C was determined by analytical gel filtration chromatography using a Superose 6 column (Amersham Biosciences), which was calibrated with the following standard proteins (Sigma) with diffusion coefficients, D20,w, (x107): sweet potato Sequence AnalysisTo compare the homology of the dimerization regions from various species, BLAST searches were used to obtain cytoplasmic dynein intermediate chain sequences. The accession numbers of the species used are: human (Homo sapiens) IC-1 AAH22540 [GenBank] , IC-2 AAC33445 [GenBank] ; mouse (Mus musculus) IC-1 AAC33445 [GenBank] , IC-2 O88487 [GenBank] ; rat (Rattus norvegicus) IC-1A X66845 [GenBank] , IC-2C U39046 [GenBank] ; red jungle fowl, (Gallus gallus) IC-1 XP_418672 [GenBank] ; dog (Canis familiaris), IC-2 XP_860043 [GenBank] ; African clawed frog (Xenopus laevis) IC-2 AAK28507 [GenBank] ; zebrafish (Danio rerio) IC-1 XP_688516 [GenBank] ; fruit fly (Drosophila melanogaster) AAX33607 [GenBank] ; fungi (Cryptococcus neoformans var. neoformans JEC21) AAW45478 [GenBank] ; nematode worm (Caenorhabditis elegans) AAC02580 [GenBank] ; filamentous fungi (Aspergillus nidulans) AAL67574 [GenBank] ; slime mold (Dictyostelium discoideum) XP_640973 [GenBank] .
The Six Intermediate Chain Isoforms Form Homo-oligomers and Hetero-oligomers in VivoInitial cloning of the intermediate chains found in rat indicated that there were two alternative splicing sites in gene 2 and that there was one alternative splicing site in gene 1, which aligned with the second alternative splicing site of gene 2 isoforms (12, 21, 31). Those data thus suggested that there were two intermediate chain 1 isoforms and three intermediate chain 2 isoforms. However, in a more recent PCR analysis of cDNAs obtained from rat brain mRNA, a second alternative splicing site was identified in intermediate chain gene 1, at a position that aligns with the first alternative splicing site of gene 2 (data not shown). Re-examination of the expression patterns of the intermediate chain 1 alternatively spliced mRNAs indicates that the brain intermediate chain isoforms originally identified as 1A and 1B (31, 32) actually correspond to the newly defined 1B and 1C isoforms, respectively (data not shown). The protein products of the two intermediate chain genes are closely related, with 79% similarity and 69% identity. The longest intermediate chains (the A isoforms) have no regions removed by splicing. There are also two genes for the cytoplasmic dynein 1 light intermediate chains in mammals. Co-overexpression and co-immunoprecipitation studies demonstrated that their gene products only form homo-oligomers (33). We therefore used a similar approach to test the ability of the six intermediate chain isoforms to form homo- or hetero-oligomers. Pairs of the intermediate chain isoforms tagged with the HA or Myc epitopes were co-expressed in cultured cells. The HA-tagged intermediate chain was then immunoprecipitated with an anti-HA antibody, and blots of the co-immunoprecipitates were probed with an anti-Myc antibody. The data presented in Fig. 1, A and B, demonstrate that both intermediate chain 1 and intermediate chain 2 isoforms can associate in all possible combinations of homo- and hetero-oligomers. In contrast, there were no detectable immunocomplexes when cells were singly transfected with Myc-tagged intermediate chains (Fig. 1C), demonstrating the specificity of the immunoprecipitating anti-HA antibody. Furthermore, no gross differences in the strength of interaction between homo-oligomers and the hetero-oligomers were observed. Region of the Intermediate Chain Involved in OligomerizationTo investigate the region of the intermediate chain necessary for oligomerization, we made intermediate chain 2C truncations tagged with the HA and Myc epitopes. The intermediate chain 2C isoform was chosen as the study protein because intermediate chain 2C is the ubiquitously expressed isoform (38) and to ensure that the analysis would not be complicated by the presence of regions that can be removed by alternative splicing. The co-immunoprecipitation assay was then used to identify the truncations of the intermediate chain that were able to oligomerize. The mapping data are summarized in Fig. 2A. A truncation containing only the predicted N-terminal coiled-coil domain, 1-65, did not oligomerize, whereas the complementary construct, 66-612, did oligomerize. A larger N-terminal truncation construct, 1-150, which includes the minimal p150 binding domain (12), also did not oligomerize. These data indicate that the N-terminal coiled-coil region is not required for oligomerization. An essential role for the seven WD repeat domains in oligomerization was excluded by the observations that truncation constructs 1-250, and 151-250 oligomerize in the co-immunoprecipitation assay.
We next investigated the requirements for the light chain binding domains in intermediate chain dimerization. Oligomerization of the 151-282 and 151-250 truncations demonstrated that the Tctex1 and LC8 light chain binding regions are not required for oligomerization. Supporting this conclusion is the observation that the 1-150 and 66-150 truncations, both of which contain these two light chain binding domains, fail to oligomerize. Co-immunoprecipitation of the 1-211 truncation suggested that the Roadblock binding region was also not required for oligomerization. In earlier work, we defined a region of the intermediate chain 1A that binds to Roadblock light chains (24). This region partially overlapped with the first WD repeat. However, it was proposed that a 39-residue fragment, corresponding to amino acid residues 211-250 of intermediate chain 2C, that excluded the portion of the first WD repeat was the actual roadblock binding domain (24). We therefore sought to confirm this hypothesis so as to more precisely define the Roadblock light chain binding region. To accomplish this, we synthesized an intermediate chain 2C construct that contains a deletion of the 39-residue fragment (amino acids 211-250) fused to the Myc epitope ( Robl). Fig. 2B demonstrates that this deletion mutant did not bind to the Roadblock light chain. This confirms that the binding site for the Roadblock light chains resides in the stretch of 39 amino acids. As a positive control, it was found that Myc-tagged full-length intermediate chain 2C interacted with GST-Roadblock but not with GST alone. Our data, and recent structural data (24, 39), demonstrate that the Roadblock light chain binds specifically to both intermediate chains and that the 39-residue fragment immediately N-terminal to the first WD-repeat domain is necessary for this binding. The ability of this Robl intermediate chain construct to oligomerize in the co-immunoprecipitation assay confirmed that the Roadblock light chain binding region is not essential for intermediate chain oligomerization (Fig. 2A). Thus, the Roadblock light chain binding domain, like the Tctex1 and LC8 light chain binding domains, is not required for oligomerization. Finally, a construct that lacks all three light chain binding domains, LC, was found competent to oligomerize. These data conclusively demonstrate that the intermediate chains are competent to oligomerize in the absence of the light chains.
Collectively, the data obtained with these intermediate chain truncation and deletion constructs indicated that amino acids 151-211 contained the region of the intermediate chain 2C that was necessary for oligomerization. To test this conclusion, the 151-211 truncation construct was synthesized and expressed in cultured cells, but it did not oligomerize. The smallest continuous intermediate chain region that was expressed in cultured cells, and which was competent to oligomerize, was the 151-250 polypeptide. One possible explanation for the failure of the 151-211 truncation polypeptide to oligomerize was that it was unable to fold into the proper conformation for oligomerization. Therefore to further investigate the role of this region in dimerization, several constructs with a deletion of this region were synthesized. Although we were unable to express a full-length construct with a deletion of that region, an N-terminal construct composed of amino acids 1-282, with a deletion of the 151-211 region ( 151-211), was expressed. The 151-211 construct was unable to oligomerize in the co-immunoprecipitation assay. As a positive control for this result, it was observed that the 1-282 truncation construct was competent to oligomerize. These data confirm the essential role of the domain defined by amino acids 151-211 in intermediate chain oligomerization. We also note that the 151-211 construct, which is unable to oligomerize, contains all three light chain binding domains. In Vitro Analysis of the Intermediate Chain Dimerization DomainTwo in vitro approaches were utilized to further analyze the interactions of the truncation polypeptides. The 151-250 polypeptide, which is the smallest polypeptide competent to oligomerize in vivo, and the 1-150 polypeptide, which does not oligomerize in vivo, were tagged with hexahistidine, expressed in bacteria, and purified. The hydrodynamic behavior of these proteins was then characterized by analytical gel filtration chromatography (Fig. 3). The dimerization-competent 151-250 truncation construct eluted as a single symmetrical peak, and no minor peaks were detected. This indicates that this polypeptide exists in a single molecular state in solution. In contrast, the 1-150 truncation protein eluted as multiple small peaks, suggesting the existence of several molecular states. The formation of these multiple molecular states may be due to the nonspecific association between monomeric units. To further characterize the purified 151-250 protein sample, it was treated with an amine-selective cross-linking reagent, DSG, with a linking spacer length of 7.72 Å (Fig. 4A) to cross-link lysine residues that possessed the appropriate accessibility and spatial orientation in solution. The high pH of the reaction buffer and low protein concentration ensured the amine selectivity of the cross-linking reaction. The native and cross-linked 151-250 polypeptides were then analyzed by SDS-PAGE and Coomassie Blue staining (Fig. 4B). A higher mass reaction product was found in the lane with the cross-linked sample. This result indicates that the purified protein also self-associates, consistent with the results obtained in the co-immunoprecipitation assay. The elution pattern of the cross-linked 151-250 protein sample was then analyzed by gel filtration chromatography (Fig. 3). The elution profile of the cross-linked protein was similar to that of the uncross-linked sample, indicating that the cross-linking reagent stabilized an existing conformation and that it did not generate higher order artifacts. Since the 151-250 protein oligomerizes in vivo, the elution of the uncross-linked protein in a single discrete peak suggests a high affinity between the monomers. The slight shift in elution pattern that was observed in the cross-linked protein sample may be the effect of the cross-linking reagent. These results with purified protein further confirm that the light chains are not required for association of the intermediate chains.
It was observed that the uncross-linked 151-250 truncation polypeptide, which has a calculated mass of
Although the uncross-linked truncation polypeptides migrate slower on SDS-PAGE than expected for their calculated molecular weights, the higher mass product found in the cross-linked 151-250 polypeptide sample (Fig. 4B, arrow) migrates close to the position predicted for the mass of a 151-250 dimer. In agreement with this observation, we found that when lysates of cells expressing Myc-tagged full-length intermediate chain 2C were treated with the cross-linking reagent, a higher mass product approximately twice that of the intermediate chain was obtained (Fig. 4C). Species Comparison of the Dimerization Domain Amino Acid SequencesThe identified dimerization domain is in the middle of the intermediate chain and C-terminal to the alternative splicing sites. Thus, it is present in the three alternative splicing isoforms derived from the intermediate chain gene 2. Since the intermediate chain polypeptides encoded by rat genes Dync1i1 and Dync1i2 are not identical, but do heterodimerize in vivo (Fig. 1), we compared the sequences of the 61 residues from the rat intermediate chains 2C and 1A (Fig. 5A). As would be expected for a region important for dimerization, there is exceptional conservation between the sequences found in the two proteins. Of the C-terminal-most 44 amino acids of the two rat proteins, all but 3 either are identical or contain conserved and semiconserved substitutions. Although the N-terminal portion is not as highly conserved, only 10 of the remaining N-terminal amino acids are non-conserved. We also compared the cytoplasmic dynein intermediate chain dimerization domain sequences of other vertebrate species. There was exceptional conservation in the dimerization region, comparable with that observed between the rat gene 1 and gene 2 sequences. We note, for example, that of the C-terminal 35 amino acids, 19 are identical in the 10 different sequences, and that there would be eight other identical columns of amino acids in the C terminus except for amino acid substitutions found in one of the analyzed sequences. When the invertebrate protein sequences are compared, there is significant conservation in the region despite their evolutionary divergence (Fig. 5B). Although there is less conservation in the N-terminal portion of the domain than was observed in the vertebrates, 20 of the C-terminal-most amino acids are either identical or conserved. We particularly note that there is near identity in the C-terminal 7 amino acids of the 15 vertebrate and invertebrate sequences.
In this study, we demonstrate that all the cytoplasmic dynein intermediate chain isoforms are competent to form hetero- and homodimers in vivo. We further identified a 61-amino-acid region of the intermediate chain upstream of the Roadblock light chain binding region that is necessary for intermediate chain dimerization in vivo. This portion of the intermediate chain does not contain the Tctex1, the LC8, or the Roadblock light chain binding domains, and it does not include the N-terminal coiled-coil domain. In vitro biochemical analyses with purified recombinant intermediate chain fragments further demonstrate that none of the three light chain families are required for the dimerization of the intermediate chains. Different intermediate chain isoforms are expressed in different tissues and cell types; for example, intermediate chain 2C is the only intermediate chain isoform found in many cultured cells, including astrocytic glia, whereas most of the isoforms are found in cultured neurons (31, 40). The expression levels of the intermediate chain isoforms are regulated during brain and neuronal development (31, 32, 41). Most interestingly, the fast and slow components of axonal transport have dynein complexes with different intermediate chain and light intermediate chain isoforms (32, 42, 43). The data reported here indicate that there is no structural impediment to the formation of dynein complexes with different combinations of the intermediate chain isoforms and provide a mechanism for generating diverse cargo binding domains. These diverse domains may create specific cargo attachment sites on the individual dynein motor units. Alternatively, the diverse domains may be permissive for differential regulation of cytoplasmic dynein by kinases and phosphatases (16, 41, 44-47). Interestingly, under experimental conditions similar to ours, the isoforms of the light intermediate chain subunit formed only homo-oligomers (33).
Prior to this work, little was known about the association of the cytoplasmic dynein intermediate chains in the dynein complex. The WD repeat domain of the intermediate chains was shown to be important for their binding to the heavy chains, and it was known that the three light chains bound to separate locations on the intermediate chain (20, 24, 26). Also, the cytoplasmic dynein complex could be dissociated by urea and other chaotropic agents into two subcomplexes. One complex contained the heavy chains and light intermediate chains, and a second complex was composed of the intermediate chains and the light chains (25, 26). The intermediate chain-light chain complex could be further fractionated into intermediate chain monomers, with and without light chains, and intermediate chain dimers, also with and without light chains (27).
The in vivo domain mapping experiments presented here clearly demonstrate that an evolutionarily conserved region, defined by residues 151-211 of intermediate chain 2C, is necessary for intermediate chain dimerization (Fig. 6). An intermediate chain construct lacking this region,
Since the
* This study was supported by grants from the NINDS National Institutes of Health and by the University of Virginia Cancer Center. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. 1 To whom correspondence should be addressed: K. Kevin Pfister, Cell Biology Dept., University of Virginia, School of Medicine, P. O. Box 800732, Charlottesville, VA 22908-0732. Tel.: 434-924-1912; Fax: 434-982-3912; E-mail: kkp9w{at}virginia.edu.
2 Cytoplasmic dynein 1 subunit nomenclature: There are two distinct cytoplasmic dynein complexes. The designations of three of the subunits of cytoplasmic dynein 1, the cytoplasmic dynein analyzed in this study, start with DYNC1. They are then grouped into families of related proteins, and an additional letter(s) was added to their names to describe their relative size as follows: DYNC1H, the heavy chain; DYNC1I, the intermediate chain family; and DYNC1LI, the light intermediate chain family. There are also three functionally distinct light chain families as follows: DYNLT, the Tctex1 family; DYNLRB, the Roadblock family; and DYNLL, the LC8 family. The light chains are not exclusive to cytoplasmic dynein 1, so they do not use the C1 designation. A dynein polypeptide subunit was identified at first mention with its formal name followed by the common name in parentheses, with most subsequent mentions using the common name: for example DYNC1LI, (light intermediate chain). The subunit families are encoded by two or more genes, and their products are distinguished by adding numbers to the name. Alternative splicing isoforms are referred to with letters, and for clarity, hyphens may be used; for example, DYNC1-I2C is the cytoplasmic dynein 1 intermediate chain 2, alternative splice isoform C. This nomenclature has been endorsed by the Human Genome Organization Nomenclature Committee (HGNC) and the International Committee on Standardized Nomenclature for Mice (5)
3 The abbreviations used are: HA, hemagglutinin; GST, glutathione S-transferase; DSG, disuccinimidyl glutarate; LC, light chain.
We thank R. John Lye, Stephen J. Susalka, and Kim Yasutis for identifying and characterizing the expression of the intermediate chain gene 1 alternatively spliced mRNAs.
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||