|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 281, Issue 15, 9901-9908, April 14, 2006
CpG Sites Preferentially Methylated by Dnmt3a in Vivo*![]() ![]() ![]() ![]() ![]() 1
From the
Received for publication, October 12, 2005 , and in revised form, January 25, 2006.
Dnmt3a and Dnmt3b are two major de novo DNA methyltransferases essential for embryonic development in mammals. It has been shown that Dnmt3a and Dnmt3b have distinct substrate preferences for certain genomic loci, including major and minor satellite repeats. However, the exact target CpG sites where Dnmt3a and Dnmt3b catalyze DNA methylation remains largely unknown. To identify a CpG site that is specifically methylated by Dnmt3a or Dnmt3b, we screened methylated genomic loci by methylation sensitive restriction fingerprinting using genomic DNA from wild-type, Dnmt3a null, Dnmt3b null, and Dnmt3a-Dnmt3b double null ES cells. Interestingly, one of the CpG sites was preferentially methylated in wild-type and Dnmt3b null ES cells but not in Dnmt3a null or Dnmt3a-Dnmt3b double null ES cells, suggesting that the site-specific methylation was Dnmt3a-dependent. Sequencing results revealed that the isolated CpG site is located within the 1st exon of the G isoform of fibroblast growth factor (Fgf-1.G) on mouse chromosome 18. Exogenous expression of Dnmt3a but not Dnmt3b in the double null ES cells restored DNA methylation of this CpG site. When we examined alternative transcription initiation sites, we determined that another CpG site in the 5'-flanking region of the Fgf-1.A isoform was also methylated specifically by Dnmt3a. Using chimeric constructs between Dnmt3a and Dnmt3b, we further determined that the NH2-terminal regulatory domain of Dnmt3a was responsible for establishing its substrate specificity. These results indicate that certain CpG sites within the Fgf-1 gene locus are preferentially methylated by Dnmt3a but not by Dnmt3b. Selective methylation by a specific member of Dnmt3 may therefore play a role in the orchestration of gene expression during embryonic development.
In mammals, global DNA methylation is catalyzed mainly by three DNA methyltransferases, Dnmt1, Dnmt3a, and Dnmt3b. Dnmt1 has a high preference for hemimethylated DNA and is essential for maintaining methylation patterns during DNA replication (1, 2). On the other hand, Dnmt3a and Dnmt3b have a preference for both unmethylated and hemimethylated DNA (3, 4) and are responsible for de novo methylation during early development. Although Dnmt3a and Dnmt3b contain highly conserved amino acid sequences in their PWWP domains, cysteine-rich domains, and carboxyl-terminal catalytic domains, the overall homology of these proteins is relatively low (3, 5). Gene disruption studies revealed that Dnmt3a null mice and Dnmt3b null mice showed different developmental defects (6). Conditional knock-outs of Dnmt3a or Dnmt3b in germ cells also exhibited distinct phenotypes (7). Dnmt3a disruption in male germ cells showed impaired spermatogenesis, while offspring from the Dnmt3a conditional mutant females died in utero. In both cases, Dnmt3a null cells displayed a lack of methylation at imprinting loci (7). In contrast, no apparent phenotype was observed in germ cell development when Dnmt3b was conditionally disrupted (7). In in vitro cell culture, a disruption of Dnmt3b, but not Dnmt3a, resulted in a partial loss of DNA methylation on endogenous C-type retroviral DNA in mouse embryonic fibroblasts (8). These results suggest that Dnmt3a and Dnmt3b have different target sites for DNA methylation depending on the cell types and the stage of development. It is known that the minor satellite DNA regions of the genome, which consist of 50,000100,000 of DNA repeats at the centromeric region, are specifically methylated by Dnmt3b but not by Dnmt3a (6, 9). In contrast, major satellite repeats, located in the pericentromeric region, are preferentially methylated by Dnmt3a (9). Furthermore, a CpG site in the Xist promoter region is revealed as a Dnmt3a-specific site of DNA methylation (9). In addition, Hsieh (10) found that the cloned EBNA region of the Epstein-Barr virus genome on the episome was a better substrate for Dnmt3a than the control luciferase coding region, suggesting that there is a substrate preference for the enzyme regardless of its chromosomal location or structure. In fact, it has been demonstrated that purified Dnmt3a protein has a DNA sequence preference when either a plasmid or an oligonucleotide is used as the methylation substrate in vitro (11, 12). In the present study, we attempted to identify the genomic DNA loci that are preferentially methylated by Dnmt3a or Dnmt3b in vivo. This was initially accomplished using methylation-sensitive restriction fingerprinting (MSRF),2 which has been successfully utilized for the identification of differentially methylated genes in cancer cells (13) or during ES cell differentiation (14). Using the MSRF method, we identified a CpG site in the G isoform of the fibroblast growth factor 1 (Fgf-1) gene locus that is preferentially methylated by Dnmt3a. This paper will further discuss the prospect that the flanking DNA sequence surrounding the CpG sites may play a role in this selective methylation by Dnmt3a.
In Vitro ES Cell DifferentiationES cells (J1, Dnmt3a null, Dnmt3b null, Dnmt3a-Dnmt3b null) were maintained on gelatin-coated dishes in ES cell maintenance medium containing 1000 units/ml recombinant mouse leukemia inhibitory factor (LIF) (ESGRO; Chemicon International) as described previously (15). ES cells were differentiated using a hanging drop method as described (15). After 48 h, embryoid bodies (EBs) were collected and attached to tissue culture plates. Medium was then replaced every other day.
MSRFThe MSRF method was performed as described previously (14). Briefly, ES cells and differentiated EBs were harvested, and genomic DNA was extracted using the Wizard genomic DNA purification kit (Promega). Extracted DNA was digested with MseI alone (New England Biolabs), or in combination with BstUI (New England Biolabs), at 10 units per 1 µg of DNA for each enzyme. The PCR reaction was performed in a 20-µl reaction mixture containing 100 ng of digested genomic DNA, 0.4 µM primers, 1.25 units of HotMasterTaq DNA polymerase (Eppendorf), 200 µM dNTP, 1 mCi/µl [ Bisulfite Sequencing and Combined Bisulfite Restriction Analysis (COBRA)DNA was extracted from ES cells, EBs, and adult tissues using the DNA Wizard genomic DNA purification kit (Promega). A bisulfite reaction was performed using the EZ DNA methylation kit (Zymo Research). Two µg of genomic DNA were used for conversion with bisulfite. Approximately 80 ng of bisulfite-converted DNA were used as a template for each PCR analysis. Primers used for COBRA analysis and bisulfite sequencing were as follows: for Fgf-1.G, 5'-TAGAAGGTAGGTTTGTAGAGGATAT-3' and 5'-ACCACACACACCAAATCAAATATTA-3'; for Fgf-1.B, 5'-TATAGTTTTTTGGGATTGGTATGGT-3' and 5'-AAAAACTCTCTTCACTCCACTTTATC-3'; for Fgf-1.A, 5'-GTTGGATAGTTGTTAGGTGAAAGGA-3' and 5'-TATTCAATAAAACCCCACAAAAAAA-3'. For bisulfite sequencing, the PCR fragments were cloned into pCRII-TOPO cloning vector (Invitrogen), and individual clones were sequenced.
Lentiviral Vector Construction and InfectionLentiviral vector, pTYF-EF-FLAG-Dnmt3a2, and FLAG-Dnmt3b1 that contains Dnmt3 under the control of the EF1- Reverse Transcriptase Polymerase Chain Reaction (RT-PCR)Total RNA was isolated from adult mouse tissues using TRIzol reagent (Invitrogen). Two µg of RNA were used as a template for RT reactions using SuperScript synthesis for the first strand synthesis kit (Invitrogen). Sequences of forward and reverse primer pairs were as follows: Fgf-1.G (5'-CTAGGAAGTAGAAGGCAGGTT-3', 5'-ACAGCTCCCGTTCTTCTTGA-3'), Fgf-1.B (5'-GTCTTGTGATAAAGTGGAGTGA-3', 5'-ACAGCTCCCGTTCTTCTTGA-3'), Fgf-1.A (5'-GCCAAGAAGCCACCTGTTAA-3', 5'-ACAGCTCCCGTTCTTCTTGA-3'). Southern BlottingDNA from wild-type ES cells, Dnmt3a-Dnmt3b null ES cells, and lentivirus-infected Dnnmt3a-Dnmt3b null ES cells were isolated using the DNA Wizard genomic DNA purification kit (Promega). Fifteen µg of DNA were digested with HpyCH4IV or HpaII (New England Biolabs), separated on 1% agarose gel, and analyzed by Southern hybridization. Probes were amplified by PCR using the following primers: for major satellite, 5'-CTGTAGGACGTGGAATATGGC-3' and 5'-CCGTGATTTTCAGTTTTCTCGC-3'; for minor satellite, 5'-TTCGAAATCGGCATTGTAGA-3' and 5'-CGGTTTCCAACGAATGTGT-3'. Amplified DNA was cloned into pCRII-TOPO cloning vector (Invitrogen) and verified by sequencing.
Identification of a CpG-rich DNA Sequence That Undergoes de Novo DNA Methylation by Dnmt3a during ES Cell DifferentiationTo isolate a DNA locus that is specifically methylated by Dnmt3a or Dnmt3b, we utilized MSRF using genomic DNA from wild-type, Dnmt3a null, Dnmt3b null, and Dnmt3a-Dnmt3b double null ES cells as the substrate. Isolated DNA was digested by a combination of methylation-sensitive (BstUI) and -insensitive (MseI) restriction enzymes and amplified by PCR. Using various combinations of four different 10-mer primers, we identified a band that is amplified less in Dnmt3a null ES cells and Dnmt3a-Dnmt3b double null ES cells (Fig. 1A). The bands were excised from the gels, and the sequences were determined. Nine out of 14 clones had the same sequence, which mapped to the first exon of the Fgf-1.G isoform on mouse chromosome 18 (Chr18:3909025239090501) and contained the BstUI site (Fig. 1B). We confirmed the results using COBRA (Fig. 2A). Bisulfite treatment converts unmethylated cytosines to thymidine residues. Therefore, both CpG sites within the BstUI recognition site (CGCG) must be methylated to be digested by BstUI in COBRA. As a template for the bisulfite-PCR, we used genomic DNA from undifferentiated ES cells or from differentiated EBs (for 2 or 4 days). Differentiation of ES cells was confirmed by monitoring the expression of marker genes, as described previously (Ref. 18 and data not shown). It has been reported that the global level of DNA methylation increased after ES cell differentiation (19), which supposedly reflects the wave of de novo methylation during early embryogenesis (20). Digestion of the bisulfite-treated PCR fragment with BstUI showed that Dnmt3a null and Dnmt3a-Dnmt3b double null ES cells had only slight DNA methylation in the Fgf-1.G-BstUI locus when they were either undifferentiated or differentiated. In contrast, wild-type and Dnmt3b null ES cells showed substantial DNA methylation when they were undifferentiated, and their methylation level was further increased after cell differentiation. We further examined more detailed DNA methylation patterns by bisulfite sequencing (Fig. 2B). As expected from the COBRA results, two CpG sequences within the BstUI site in the first exon of Fgf-1.G (+276 and +278) were methylated in wild-type cells. Similarly, both CpG sites in the BstUI site were highly methylated in Dnmt3b null ES cells. In contrast, a significant decrease in methylation of CpG within the BstUI site was observed in Dnmt3a null ES cells. Interestingly, the second CpG site within the BstUI site (+278) was more dominantly hypomethylated (87.5% of clones) than the first CpG (+276, 50% of clones) in Dnmt3a null ES cells. In COBRA, BstUI could digest bisulfite-treated PCR products only if both of the CpG sites are methylated. Therefore, only 6% of the fragments for Dnmt3a null ES cells are expected to be digested with BstUI, which appears consistent with the COBRA data above (Fig. 2A). In Dnmt3a-Dnmt3b double null ES cells, CpG sites were totally unmethylated. We then analyzed the CpG sites neighboring the BstUI site. As shown in Fig. 2B, two CpG sites upstream of the BstUI site (+232 and +247) were highly methylated in wild-type cells. These upstream CpG sites were relatively highly methylated in Dnmt3a null ES cells. These results suggest that some CpG sites within the Fgf-1 gene are preferentially methylated by Dnmt3a in vivo.
DNA Methylation in Other Isoforms of Fgf-1 Gene during ES Cell DifferentiationSince the murine Fgf-1 gene has multiple isoforms resulting from different transcription initiation sites (21), we analyzed CpG methylation within other Fgf-1 gene isoforms using COBRA (Fig. 1). The Fgf-1.A isoform is predominantly expressed in the heart, and the Fgf-1.B isoform is mainly expressed in the brain (21). A HpyCH4IV site within the 5'-upstream region of Fgf-1.A isoform (Fgf-1.A-HpyCH4IV: Chr18:3915336739153370) and a HinfI site within the first exon of Fgf-1.B isoform (Fgf-1.B-HinfI: Chr18:3914239139142395) were identified as sites suitable for COBRA. Interestingly, the Fgf-1.A-HpyCH4IV site was also hypomethylated in Dnmt3a null ES cells as well as the Fgf-1.G-BstUI site. In contrast, methylation of the Fgf-1.B-HinfI site was slightly higher in Dnmt3a null ES cells when compared with Dnmt3b null ES cells (Fig. 2A). Fgf-1 DNA Methylation Patterns in Adult TissuesWe also examined the gene expression and DNA methylation status of Fgf-1 isoforms in mouse adult tissues (Fig. 3, A and B), since DNA methylation is generally associated with gene inactivation. The RT-PCR results revealed that the Fgf-1.G gene was highly expressed in liver and moderately expressed in kidney, as reported previously (22). DNA methylation analysis by COBRA showed that the Fgf-1.G-BstUI site was hypomethylated in liver, where Fgf-1.G gene expression was high. In kidney, the Fgf-1.G-BstUI site methylation was lower than in other tissues that do not express Fgf-1.G. Bisulfite sequencing further confirmed the COBRA data (Fig. 3C). Overall, there was a tendency for DNA methylation levels to be higher in tissues where the Fgf-1.G gene was not expressed. Similarly, the Fgf-1.A-HpyCH4IV site was hypomethylated in heart and kidney, where the Fgf-1.A gene was expressed. The correlation was less clear with Fgf-1.B. The Fgf-1.B-HinfI locus was highly methylated in testis and only slightly methylated in other somatic organs.
Reconstitution of DNA Methylation by Ectopic Expression of Dnmt3aTo further demonstrate that Dnmt3a is responsible for site-specific DNA methylation within the Fgf-1 gene, we reconstituted the expression of Dnmt3a in Dnmt3a-Dnmt3b double null ES cells and looked for a recovery in methylation (Fig. 4A). To do this, we created stable clones expressing either Dnmt3a or Dnmt3b in Dnmt3a-Dnmt3b double null ES cells. We confirmed that these clones could be maintained as undifferentiated ES cells and could differentiate as EBs (data not shown). The status of DNA methylation in the ES cell clones was then analyzed by COBRA. As expected, DNA methylation of the Fgf-1.G-BstUI site and the Fgf-1.A-HpyCH4IV site was preferentially restored by ectopic expression of Dnmt3a but not Dnmt3b. The Fgf-1.G-BstUI site DNA methylation was not observed when transduced ES cells were maintained in the undifferentiated condition. In contrast to these two sites, the DNA methylation of the Fgf-1.B-HinfI site was increased by the ectopic expression of either Dnmt3a or Dnmt3b. To rule out the possibility that the expressed Dnmt3b is not or less functional, we examined DNA methylation of minor satellite repeats in those clones (Fig. 4B). Methylation of minor satellite repeats in double null ES cells was more efficiently restored by Dnmt3b than Dnmt3a reconstitution. Therefore, we conclude that both Dnmt3a and Dnmt3b used in this study are functional and exhibit distinct substrate specificity. The Regulatory Domain of Dnmt3a Is Important for the Site-specific DNA MethylationDnmt3a and Dnmt3b each contain a highly conserved catalytic domain at their carboxyl termini. In contrast, the amino-terminal halves of the proteins are more divergent and show low homology. To examine which domain is more important to the substrate specificity, we made a chimeric construct of Dnmt3a and Dnmt3b, namely, Dnmt3a-3b, which contains the Dnmt3a regulatory domain with the Dnmt3b catalytic domain, and Dnmt3b-3a, which contains the Dnmt3b regulatory domain with the Dnmt3a catalytic domain (Fig. 5A). These chimeric genes were ectopically expressed in Dnmt3a-Dnmt3b double null ES cells using lentiviral vectors, and their DNA methylation activity was analyzed. As shown in Fig. 5B, DNA methylation in major satellite DNA repeats was preferentially recovered by ectopic expression of Dnmt3a as reported previously (9). Dnmt3a-3b, but not Dnmt3b-3a, partially recovered DNA methylation in major satellite DNA. Therefore, the regulatory domain of Dnmt3a is more important than the catalytic domain to determine the substrate specificity of Dnmt3a to major satellite repeats. Using these constructs, we examined the DNA methylation of Fgf-1 sites (Fig. 5C). The Fgf-1.A-HpyCH4IV site was methylated by the transduction of Dnmt3a-3b in Dnmt3a-Dnmt3b double null ES cells regardless of differentiation status. In contrast, Dnmt3b-3a expression did not restore DNA methylation. These results suggest that the regulatory domain, rather than the catalytic domain, of Dnmt3a is important for the site-specific DNA methylation at the Fgf-1.A site. It should be noted that methylation of the Fgf-1.G-BstUI site was not observed by the expression of either one of the chimeric constructs. Fgf-1.B-HinfI site was similarly methylated by Dnmt3a, Dnmt3b, and Dnmt3a-3b (data not shown).
Using MSRF, we have identified a BstUI recognition site within the first exon of the Fgf-1.G gene locus as a CpG site that is methylated preferentially by Dnmt3a. Similarly, we found that a HpyCH4IV site in the 5'-flanking region of the Fgf-1.A gene was predominantly methylated by Dnmt3a. Furthermore, reconstitution of Dnmt3a in double null ES cells restored the methylation of these sites, whereas the expression of Dnmt3b did not. These results demonstrated that these CpG sites in the Fgf-1 gene are specific target sites for Dnmt3a.
Genome-wide methylation analysis by restriction landmark genomic scanning, using the NotI restriction enzyme, suggested that methylation patterns at CpG islands are quite similar between wild-type, Dnmt3a null, and Dnmt3b null ES cells (23). However, by MSRF screening, we identified the genomic locus that is specifically methylated by Dnmt3a. Even though genome-wide screening has a clear merit, restriction landmark genomic scanning restricts its analysis for the CpG sites located within the recognition site of the restriction enzyme used. In contrast, MSRF screening could be useful for the analysis of CpG methylation, regardless of flanking sequences. It remains undetermined how Dnmt3a catalyzes DNA methylation of the specific target sites in the genome. One hypothesis is that Dnmt3a binds to the sequence-specific DNA binding molecules, which brings the complex to the specific genomic locus. In fact, Dnmt3a has been shown to bind to a sequence-specific DNA-binding protein (2426) and act as a co-repressor (25, 26), although it has not been demonstrated whether Dnmt3a did actually catalyze DNA methylation at this specific site. The other hypothesis is that Dnmt3a itself possesses a substrate DNA sequence preference. It has been reported that the cloned EBNA region of the Epstein-Barr virus genome on the episome was a better substrate in vivo for Dnmt3a than a control luciferase coding region on the same episome (10). Handa et al. (11) reported that there is a preference for the flanking sequence surrounding methylated CpG sites for Dnmt3a when they analyzed DNA methylation in vitro using a purified Dnmt3a and oligonucleotide DNA substrates. According to the study, the RCGY consensus sequence was a more preferred substrate for methylation by Dnmt3a. In contrast, YCGR was not a preferred substrate for Dnmt3a. It should be noted that, however, the study has not demonstrated a differential sequence preference between Dnmt3a and Dnmt3b. In our study, all the sites preferentially methylated by Dnmt3a, namely, two CpG sites within Fgf-1.G-BstUI site and one at the Fgf-1.A-HpyCH4IV site, contain RCGY sequences. Furthermore, the Fgf-1.B-HinfI site, which was not selectively methylated by Dnmt3a, contains the YCGR sequence. Thus, our results support the previous in vitro studies that Dnmt3a may possess a substrate preference to RCGY even in vivo. Of interest, the Fgf-1.B-HinfI YCGR site was more highly methylated in Dnmt3a null ES cells than in Dnmt3b null ES cells (Fig. 2A). Furthermore, Dnmt3b reconstitution caused a slightly higher degree of DNA methylation at the Fgf-1.B-HinfI YCGR site than Dnmt3a (Fig. 4A). Therefore, the YCGR site might be a preferred substrate for Dnmt3b in vivo. In fact, Dnmt3b is known to be highly expressed in adult testis (3, 5), where the methylation level of Fgf-1.B-HinfI site was very high (Fig. 3B). The major satellite repeats selectively methylated by Dnmt3a (Fig. 5) include ACGT (HpyCH4IV or MaeII recognition site), which consists of a RCGY sequence. The minor satellite repeats, which are known to be selectively methylated by Dnmt3b, on the other hand, contain a YCGR sequence. Taken together, these data suggest that the substrate preference of Dnmt3a or Dnmt3b may play an important role in determining which genomic loci are differentially methylated in vivo. Dnmt3a and Dnmt3b each consist of an NH2-terminal regulatory domain, which contains the PWWP domain, a cysteine-rich domain, and a COOH-terminal catalytic domain (3, 5). The catalytic domain of either Dnmt3a or Dnmt3b by itself has been shown to possess a similar substrate preference as a full-length Dnmt3a in vitro (11). Therefore, the catalytic domain may not be critical to determine the substrate specificity of Dnmt3a and Dnmt3b in vivo. Indeed, our experiments using chimeric Dnmt3 constructs between Dnmt3a and Dnmt3b revealed that the amino-terminal regulatory domain of Dnmt3a was more important for the site-specific methylation of major satellite repeats and the Fgf-1.A-HpyCH4IV locus (Fig. 5, B and C). The PWWP domain has been shown to have a high affinity for DNA (2729) and is required for the methylation of the major satellite repeats (29). However, considering the non-sequence-specific binding property of the Dnmt3 PWWP domain (29), it may be difficult to assume that this domain is directly involved in enzyme targeting to specific DNA loci, such as the Fgf-1 gene. FGF-1 is a member of the fibroblast growth factor family, which possesses broad mitogenic activities, and is involved in a variety of biological processes, including cell proliferation, migration, differentiation, and survival (3035). The Fgf-1 gene contains multiple promoters associated with distinct first exons, which result in multiple transcript isoforms containing common sequences except for the first exons (36). These promoters are regulated in a tissue-dependent manner both in mouse and human (21). Since ectopic expression of Fgf-1 induces intimal hyperplasia in vivo (37), their strict gene regulation appears to be important in normal embryonic development and tissue homeostasis. In this paper, we also showed that the degree of DNA methylation inversely correlates with their expression pattern in adult mouse tissues. DNA methylation may be one of the mechanisms involved in the silencing of gene expression by preventing the access of transcription factors. Our data imply that a particular Dnmt3 member might play a role in the regulation of isoform determination of the Fgf-1 gene during embryonic development. In conclusion, we have identified CpG sites preferentially methylated by Dnmt3a in vivo. Although the exact molecular mechanisms underlying this selective methylation remain undetermined, the regulatory domain of Dnmt3a likely plays a predominant role. Further, common flanking sequences (RCGY) around CpG may be involved in this selective methylation by Dnmt3a in vivo. Site-specific methylation by a specific member of Dnmt3 may play a role in the orchestration of gene expression during embryonic development.
* This work was supported in part by National Institutes of Health Grants DK59699 and RR17001 (to N. T.). 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: Dept. of Pathology, University of Florida College of Medicine, 1600 SW Archer Rd., Gainesville, FL 32610. Tel.: 352-392-2696; Fax: 352-392-6249; E-mail: terada{at}pathology.ufl.edu.
2 The abbreviations used are: MSRF, methylation-sensitive restriction fingerprinting; EB, embryoid body; COBRA, combined bisulfite restriction analysis; RT, reverse transcriptase.
We thank En Li for providing us with Dnmt null ES cells and Keith Robertson, Amy Meacham, Sarah Kehoe, Jeffrey Brower, and Catherine Ketcham for helpful discussions and critical reading of the manuscript.
This article has been cited by other articles:
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||