Cloning and Characterization of a Novel Endothelial Promoter of the Human CYP19 (Aromatase P450) Gene that Is Up-Regulated in Breast Cancer Tissue
Siby Sebastian1,
Kazuto Takayama,
Makio Shozu and
Serdar E. Bulun
Departments of Obstetrics and Gynecology and Molecular Genetics (S.S., S.E.B.), University of Illinois at Chicago, Chicago, Illinois 60612; Department of Obstetrics and Gynecology (K.T.), Tohoku University School of Medicine, Sendai 980-8574, Japan; and Department of Obstetrics and Gynecology (M.S.), School of Medicine, Kanazawa University, 13-1 Takara-machi, Kanazawa, Japan 920-0934
Address all correspondence and requests for reprints to: Serdar E. Bulun, M.D. Departments of Obstetrics and Gynecology and Molecular Genetics, The University of Illinois at Chicago, 820 South Wood Street, M/C 808, Chicago, Illinois 60612. E-mail: sbulun{at}uic.edu.
 |
ABSTRACT
|
---|
Intratumoral expression of aromatase P450 (P450arom) promotes the growth of breast tumors via increased local estrogen concentration. We cloned a novel 101-bp untranslated first exon (I.7) that comprises the 5'-end of 2954% of P450arom transcripts isolated from breast cancer tissues (n = 7). The levels of P450arom transcripts with exon I.7 were significantly increased in breast tumor tissues and adipose tissue adjacent to tumors. We identified a promoter immediately upstream of exon I.7 and mapped this to about 36 kb upstream of ATG translation start site of the CYP19 (aromatase cytochrome P450) gene. Sequence analysis of I.7 revealed a TATA-less promoter containing an initiator, two consensus GATA sites, and cis-regulatory elements found in megakaryocytes and endothelial type promoters. Luciferase activity directed by the promoter I.7 sequence (-299/+81 bp) was 4-fold greater than a minimum length promoter sequence (-35/+81 bp) in human microvascular endothelial cells (HMEC-1), but only 2-fold greater in MCF-7 breast malignant epithelial cells. There was no promoter activity in primary breast adipose fibroblasts. Site-directed mutations demonstrated that maximal basal promoter activity required two GATA motifs at -146/-141 bp and -196/-191 bp. Gel shift and deoxyribonuclease I footprinting assays demonstrated the binding of GATA-2 transcription factor but not GATA-1 to the -196/-191-bp region. Overexpression of GATA-2 in HMEC-1 cells increases promoter I.7 activity by 5-fold. In conclusion, promoter I.7 is a GATA-2-regulated endothelial promoter of the human CYP19 gene and may increase estrogen biosynthesis in vascular endothelial cells of breast cancer. The activity of this promoter may also be important for intracrine and paracrine effects of estrogen on blood vessels.
 |
INTRODUCTION
|
---|
LOCAL BIOSYNTHESIS OF estrogen, mediated by the elevated intratumoral expression of aromatase cytochrome P450 (P450arom), promotes the growth of breast tumors (1, 2, 3, 4). The clinical relevance of these observations was exemplified recently by the successful treatment of breast carcinomas with potent aromatase inhibitors (5, 6). In fact, aromatase inhibitors may replace tamoxifen as the first-line endocrine therapy of breast cancer in the near future (7). P450arom, the key protein for estrogen biosynthesis, is encoded by a single copy of the CYP19 gene localized at human chromosome 15q21.2 (8). Thus far, a total of nine distinct tissue-specific promoters and related untranslated first exons designated as I.1 (placenta-major), I.2 (placenta-minor), I.3 (adipose/breast cancer), I.4 (skin/adipose tissue), I.5 (fetal tissue), I.6 (bone), PII (ovary/breast cancer/endometriosis), 2a (placenta-minor), and I.f (brain) have been identified within the approximately 90-kb region 5'-upstream of the coding region of the CYP19 gene (8, 9). Alternative use of each promoter, which regulates mature P450arom mRNA levels by splicing of each untranslated first exon [5'-untranslatory region (UTR)] onto a common splice junction 38 bp upstream of the translation start site (ATG) of the first coding exon (exon II), is the underlying molecular mechanism conferring tissue-specific regulation of aromatase expression (9, 10).
Several studies have demonstrated concentration gradients of P450arom expression within the breast bearing a tumor (1, 2, 3). In breast cancer patients, total P450arom transcript levels were significantly higher in adipose tissue proximal to a tumor in comparison with adipose tissue distal to a tumor (1, 2, 3, 11, 12). Previous studies from this laboratory and elsewhere have demonstrated that a switch in the use of CYP19 gene promoters in the breast adipose fibroblasts surrounding malignant epithelial cells is, at least in part, responsible for the excessive P450arom expression noted in breast cancer (3, 12, 13, 14). The current report describes identification of a novel 101-bp first exon that comprises the 5'-UTR of a significant portion of P450arom transcripts using 5'-rapid amplification of cDNA ends (RACE) generated libraries of human breast cancer tissues.
2 Subsequent cloning and functional characterization of the genomic region upstream of exon I.7 revealed a novel promoter regulated by transcription factors (e.g. GATA-2) typically found in endothelial cells. The activity of promoter I.7 located about 36 kb upstream of the human CYP19 (P450arom) gene coding region gives rise to splicing of this 101-bp exon onto a common splice site 38 bp upstream of the translation start codon (ATG) of the P450arom mRNA.
 |
RESULTS
|
---|
Identification of a Novel First Exon of the Human CYP19 Gene
With an objective to identify all promoters of the CYP19 gene (including previously unknown ones) contributing to the increased levels of total P450arom mRNA in breast cancer tissue, we have conducted 5'-RACE experiments using total RNA isolated from seven separate human female breast cancer samples. Direct sequence analysis of randomly selected clones from two 5'-RACE libraries generated from two separate human breast cancer samples revealed a novel untranslated first exon present in both libraries. These novel clones contained a hitherto unknown 101-bp sequence comprising the 5'-untranslated ends of P450arom mRNA transcripts (Fig. 1
). A basic local alignment search tool (BLAST) search analysis mapped this 101-bp sequence onto a region 36,343 bp 5'-upstream of the exon II of the human CYP19 gene (Fig. 1
). Following the traditional nomenclature adopted to identify various untranslated first exons of human CYP19 gene, we named this 101-bp novel exon as exon I.7. In these two cancer RNA samples, sequencing analysis also revealed previously identified 5'-UTRs, namely promoter II-, I.3-, and I.4-specific sequences (see below).

View larger version (16K):
[in this window]
[in a new window]
|
Figure 1. Sequence (101 bp) and Genomic Location of Exon I.7 of the Human CYP19 (P450arom) Gene
The numbers over a specific nucleotide indicate 5'-RACE clones whose 5'-end corresponded to that particular base. The common splice acceptor site is located 38 bp upstream of the translation start site (ATG) in exon II of P450arom transcripts. The genomic locations of exon I.4 and I.1, two major CYP19 gene promoters expressed in adipose tissue and placenta, respectively, are also noted.
|
|
Abundance of P450arom Transcripts Carrying Exon I.7 in Human Breast Tissues
To assess the abundance of exon I.7-specific P450arom transcripts, we performed additional 5'-RACE and colony hybridization experiments. First, we determined the distribution of P450arom transcripts with promoter-specific 5'-UTRs by sequencing 14 P450arom clones in a 5'-RACE library (Table 1
). Four of 14 clones (
29%) contained exon I.7. The other P450arom clones contained promoter II, exon I.3, and exon I.4, which was consistent with previously published data (11, 12, 13, 14). We found only exon I.4-containing P450arom transcripts in a cancer-free breast tissue, used as negative control (Table 1
). Next we prepared replica plates of bacterial colonies representing the library of 5'-RACE clones of human breast cancer RNA samples from four additional patients. As summarized in Table 2
, the colony hybridization experiments demonstrated that of a total of 270 breast cancer tissue P450arom transcripts, 114 (or
40%) transcripts carried exon I.7. The ratio of I.7-specific transcripts to total transcripts varied from 29%54% in samples represented in Table 1
.
Exon I.7-Specific P450arom Transcripts in Breast Cancer and Adjacent Adipose Tissues
We tested many more breast cancer samples and adjacent adipose tissue to understand the pathological significance and distribution of exon I.7-specific P450arom transcripts. First, using RT-PCR employing exon III- and I.7-specific oligonucleotides, we demonstrated the presence of I.7-specific P450arom transcripts in cancer samples from six patients (part of the data not shown). We then employed semiquantitative RT-PCR to amplify I.7-specific P450arom transcripts in 1) cancer tissue, 2) adipose tissue samples biopsied at 1 cm, 2 cm, 7 cm, or 10 cm to cancer, and 3) adipose tissue samples from cancer-free individuals. We determined that there is a positive concentration gradient of I.7-specific transcripts in mastectomy specimens bearing tumors. The highest I.7-specific P450arom transcript levels were found in cancer tissues or adipose tissue adjacent to cancer. Breast adipose tissue from cancer-free women (n = 3) did not contain I.7-specific transcripts. Interestingly, abdominal adipose tissue from a pregnant woman and buttock adipose tissue from a normal 33-yr-old man contained high levels of I.7-specific P450arom transcripts (part of the data illustrated in Fig. 2
). We should note that the minor RT-PCR products represented as faint bands larger than the expected size for exon I.7-specific product probably represent minor splice variants. We have not come across any unexpected sequences although we sequenced a total of 22 exon I.7-specific 5'-RACE clones generated from three separate breast cancer tissues and one male buttock adipose tissue (data not shown; Table 1
and Fig. 1
). Because these faint bands could be detected by Southern hybridization but not on agarose gels, we did not have a chance to isolate and sequence them.

View larger version (47K):
[in this window]
[in a new window]
|
Figure 2. RT-PCR-Southern Blot Expression Analysis of P450arom Transcripts Containing Exon I.7
RT-PCR to amplify P450arom transcripts with the exon I.7-specific sequence was conducted using total RNA from breast tumor, breast adipose tissue adjacent to the tumor, breast adipose tissue from a cancer-free woman, abdominal adipose tissue from a pregnant woman, and buttock adipose tissue from a normal man. A 593-bp region of human GAPDH mRNA was amplified as an experimental control. M, Molecular weight marker lane.
|
|
To examine whether the expression of exon I.7 occurs in any other human tissues, we have conducted exon I.7-specific RT-PCR experiments using first-strand cDNA prepared from normal human brain, lung, kidney, heart, pancreas, placenta, and skeletal muscle as templates (Human Multiple Tissue cDNA panel No. K1420-1, CLONTECH Laboratories, Inc., Palo Alto, CA). No exon I.7-specific PCR product was detected in any of these tissues even after 40 cycles of PCR amplification (data not shown).
Identification of the Transcription Initiation Site
Based on the initial sequence information obtained from the longest exon I.7 clone found in the 5'-RACE-generated library, the most distal 5'-nucleotide was an adenine, defining the total length of exon I.7 as 101 bp. To ascertain whether additional transcription start sites were present, we generated an exon I.7-specific 5'-RACE sublibrary of P450arom mRNA transcripts expressed in male buttock adipose tissue. All but one of the eight randomly selected clones of PCR products subjected to direct sequencing had the previously identified 5'-boundary, i.e. an adenine nucleotide positioned 101 bp upstream of the guanine of the splice acceptor site of coding exon II (Fig. 1
). These results are indeed in agreement with a potential transcription initiator (Inr) noted in the promoter I.7 sequence (Fig. 3
and see below).

View larger version (45K):
[in this window]
[in a new window]
|
Figure 3. Sequence and Structural Features of a 800-bp DNA Fragment Comprising the -719/+81-bp 5'-Flanking Region of Exon I.7
Sequences with similarity to consensus binding motif of known transcriptional factors are boxed. GATA elements potentially involved in transcriptional regulation of the promoter I.7 are shown in bold. All numbering is relative to the transcription initiation site Adenine (indicated by +1), as defined by 5'-RACE experiments described in Materials and Methods. Exon I.7 sequence is underlined.
|
|
Cloning and Characterization of a Genomic Fragment Encompassing 5'-Flanking
Region of Exon I.7
Employing exon I.7-specific oligonucleotides as probes, a human genomic library was screened. A fragment of 5'-flanking sequence of exon I.7 (
2.4 kb) was amplified by PCR using a positive phage (EMBL3) clone containing a genomic sequence (
15 kb) as the template. Sequence analysis confirmed the presence of exon I.7 in this clone. Incidentally, a BLAST search analysis aligned this genomic sequence to a region about 36 kb upstream of the first coding exon (exon II) of the human CYP19 (P450arom) gene. A number of smaller fragments of this clone encompassing the 5'-flanking region of exon I.7 were derived by restriction digestion and/or PCR to generate various deletion mutants fused to luciferase reporter constructs described below.
Structural Features of Promoter I.7
Figure 3
depicts the DNA sequence of the 719-bp 5'-flanking region of exon I.7. A computational analysis (using the TFSEARCH database available at http://www.cbrc.jp/research/db/TFSEARCH.html) of the 800-bp fragment spanning the -719/+81-bp region of exon I.7 revealed the presence of several consensus binding sequences for a number of trans-activating regulatory proteins. The proximal 5'-flanking region was devoid of a canonical TATAA or CAAT element. Notably, an initiator (Inr) or cap signal for transcription initiation at -4/+6 bp and two consensus GATA binding motifs at positions -146/-141 bp and -196/-191 bp were predicted (15, 16). Some other cis-acting motifs predicted for promoter I.7 include two consensus c-Ets sites (GGAA) at -166/-163 and -285/-282 bp, two reverse c-Ets sites (TTCC) at -126/-123 bp and -182/-179 bp, an E47 site at -247/-243 bp, a sterol response element binding protein-1 site at -454/-443 bp, a hepatocyte nuclear factor-3ß site at -601/-590 bp, and an Oct-1 site at -637/-624 bp.
Functional Analysis of Promoter I.7 in Various Cell Types
The transcriptional activity of promoter I.7 was assayed by ligating selected segments of the 5'-flanking region of exon I.7 upstream of the luciferase reporter gene in pGL3-Basic vector. These deletion constructs were then transiently transfected into various cell types, as described in Materials and Methods. Because we observed elevated levels of exon I.7-specific P450arom transcripts in breast cancer tissues, we first examined the basal activity of promoter I.7 in the breast cancer cell line MCF-7. Interestingly, transfection of MCF-7 cells with equal molar amounts of various reporter constructs of promoter I.7 resulted only in a marginal increase (
2-fold) of maximal promoter activity (-299/+81 bp) compared with the minimal promoter construct containing a -35/+81-bp sequence (Fig. 4
). On the other hand, transfection of normal female breast adipose fibroblasts did not show any promoter activity at all (data not shown). We should note that we routinely transfect primary breast adipose fibroblasts with P450arom promoter II luciferase constructs to study regulation of this promoter and observe up to about a 6-fold increase in basal activity in serial deletion constructs (17).

View larger version (15K):
[in this window]
[in a new window]
|
Figure 4. Functional Analysis of Promoter I.7 in MCF-7 Breast Malignant Epithelial Cells
Luciferase reporter plasmids (in equimolar amounts) containing serially truncated putative promoter I.7 fragments (shown schematically on the left) or vector alone were transiently transfected into MCF-7 cells. Firefly luciferase activities were normalized by accompanied activities of cotransfected Renilla luciferase reporter gene (see Materials and Methods). These data are representative of two independent experiments, each performed in triplicate.
|
|
Promoter I.7 Activity in Endothelial Cells
Previously, several investigators have described the presence of functionally active multiple GATA binding motifs common in the promoters that were transcriptionally regulated in megakaryocytic/endothelial cells (18, 19, 20, 21). Additionally, immunolocalization of P450arom in endothelial cells (human coronary arteries, vena cava, and mouse vessel wall) have also been demonstrated (22, 23, 24). All these observations prompted us to examine the possibility that promoter I.7 may be the endothelial-type promoter of the human CYP19 gene.
As shown in Fig. 5A
, transient transfection of equal molar amounts of various deletion constructs of the putative promoter I.7 into the human microvascular endothelial cells (HMEC-1) resulted in substantial increases (up to
4-fold) in the luciferase activity compared with a promoterless pGL3 basic vector or the construct containing a minimal fragment (-35/+81 bp) of promoter I.7. The highest promoter activity was observed in the reporter construct containing the -299/+81-bp region. Deletion up to the -35-bp position resulted in a complete loss of promoter activity, indicating that the sequence from -299 to -35 bp (which included two GATA sites) is critical for the basal promoter activity. On the other hand, addition of 5'-upstream sequences up to position -2344 bp gradually reduced the promoter activity (Fig. 5A
).
The GATA Binding Site at -196/-191 bp Is Critical for Promoter I.7 Activity
The results obtained from transfection experiments described above indicated that the -299/-35-bp region is critical for promoter I.7 activity. Sequence analysis had shown the existence of two consensus GATA binding sites at -196/-191 bp and -146/-141 bp. To investigate the roles of these sequences in directing basal promoter I.7 activity, site-directed mutations (separately or in combination) were introduced into these putative sites, and functional consequences of these mutations were then examined by transient transfection of reporter plasmid constructs into HMEC-1 cells. Mutation of the -146/-141-bp GATA site from the consensus sequence GATA
ctcg did not affect the promoter activity significantly (Fig. 5A
). In contrast, mutation of the -196/-191-bp GATA site (GATA
ctcg) resulted in 55% reduction in luciferase activity (Fig. 5A
). Mutations of both of these GATA sites abolished promoter I.7 activity in HMEC-1 cells. Taken together, these results demonstrated that even though both GATA motifs are required for optimal promoter activity, the GATA site at -196/-191 bp is more critical for sustaining the basal activity of promoter I.7 in endothelial cells.
Human GATA-2 Transactivates Promoter I.7 in a Concentration-Dependent Manner in HMEC-1 Cells
To determine whether GATA proteins could modulate promoter I.7 activity, we introduced the -299/+81 promoter I.7 reporter plasmid construct along with the mammalian expression vector pMT2-hGATA-2 (18), which contained a full-length human GATA-2 cDNA into HMEC-1 cells. As shown in Fig. 5B
, cotransfection of pMT2-hGATA-2 (100 ng) resulted in approximately 5-fold stimulation in the luciferase activity of the reporter plasmid -299/+81 in HMEC-1 cells compared with the empty expression vector used as a control. Moreover, addition of increasing amounts (10100 ng/well) of pMT2-hGATA-2 expression vector elicited a concentration-dependent increase in promoter I.7 activity in HMEC-1 cells. Taken together, these results demonstrate the ability of GATA-2 to transactivate P450arom promoter I.7 in a concentration-dependent fashion in endothelial cells.
The GATA Binding Site at -196/-191 bp Is Protected by HMEC-1 Nuclear Protein(s) in Deoxyribonuclease I (DNase I) Footprinting
In vitro footprinting revealed that a nuclear extract from HMEC-1 cells protected a sequence flanking the -196/-191-bp GATA site from digestion by DNase I (Fig. 6
). Two additional footprints at an E47 and Ets binding sites positioned at -247/-243 bp and -166/-163 bp, respectively, were also detected (Fig. 6
).

View larger version (77K):
[in this window]
[in a new window]
|
Figure 6. DNase I Footprint Analysis of Promoter I.7
An end-labeled 196-bp SacIXbaI fragment comprising the -287/-92-bp promoter I.7 sequence was incubated with varying amounts of nuclear extract from HMEC-1 cells and subsequently digested with DNase I. The regions of the DNA probe protected from digestion are indicated by boxes. We added 0.5 µg of BSA to the sample that did not include nuclear extract (lane 3). M, End-labeled X174 DNA/HinfI markers (Promega Corp.); G+A, Maxam-Gilbert sequencing ladder of the probe DNA.
|
|
GATA-2 Binds to the -196/-191-bp GATA Site in Promoter I.7
To identity the protein(s) that binds at the -196/-191-bp GATA site, we performed EMSA using a 30-bp double-stranded oligonucleotide probe (representing the genomic region -209/-180 bp in promoter I.7) encompassing the putative GATA site at -196/-191 bp. As shown in Fig. 7
, incubation of nuclear extracts from HMEC-1 cells with the 32P-labeled oligonucleotide probe resulted in the formation of one specific band (lane 2), which was abolished by 50-fold molar excess of unlabeled probe (lane 3). In contrast, preincubation with an oligonucleotide, which also represents the sequence at -209/-180 bp but with a mutation in the -196/-191-bp GATA site (GATA
ctcg), did not compete with the wild-type probe (lane 4), demonstrating that the -196/-191-bp GATA sequence is a critical requirement for binding activity.

View larger version (68K):
[in this window]
[in a new window]
|
Figure 7. Specific Binding of the Transcription Factor GATA-2 to the -196/-191-bp GATA Element in Promoter I.7
Nuclear extracts (250 ng) derived from HMEC-1 cells were incubated with 32P- labeled double-stranded oligonucleotide identical to the -209/-180-bp region of promoter I.7 sequence containing the GATA consensus binding site at -196/-191 bp (lane 2). Competition reactions were performed using 50-fold molar excess of unlabeled wild-type sequence (lane 3) or with a mutated GATA binding site (lane 4). Supershift experiments were conducted by preincubating the nuclear extracts with antibodies raised against either GATA-1 (lane 5), GATA-2 (lane 6), or an unrelated protein (E47, lane 7). The specific complex (left arrow) could be supershifted only by the anti-GATA-2 antibody (right arrow).
|
|
To identify the transcription factor that binds to the -196/-191-bp motif, supershift experiments using antibodies directed against either GATA-1 or GATA-2 (Santa Cruz Biotechnology, Inc., Santa Cruz, CA) were performed. Preincubation of nuclear extracts with anti-GATA-2 antibody but not the anti- GATA-1 antibody supershifted the protein-DNA complex to a higher region in the gel (Fig. 6
, lanes 5 and 6, respectively). In contrast, the addition of a control antibody, anti-E47, had no effect (lane 7), confirming the specific binding of GATA-2 to the -196/-191-bp GATA site in promoter I.7.
 |
DISCUSSION
|
---|
We report herein a novel promoter that directs aromatase expression preferentially in breast cancer tissues. This promoter (I.7) is located approximately 36 kb upstream of the human CYP19 (P450arom) gene coding region and gives rise to a specific P450arom mRNA species, possibly by alternative splicing. Intriguingly, we also present evidence that promoter I.7 contains endothelial cell-specific cis-acting elements, and that significant levels of its activity could be demonstrated in an endothelial cell line but not in breast adipose fibroblasts or malignant epithelial cells.
The GATA transcription factor family consists of six homologous members, and, by virtue of its highly conserved zinc finger DNA binding domains, they all potentially recognize a (G/A) GAT (T/A)(A/G) motif (25, 26). Even though the GATA proteins display distinct but overlapping patterns of expression in a variety of tissues, gene disruption experiments revealed critical requirements for GATA-1 and GATA-2 in hematopoiesis and megakaryopoiesis (25, 26). In addition to a high level of expression in hematopoietic stem cells, mast cells, and megakaryocytes, abundant expression of GATA-2 has also been noted in vascular endothelial cells (26).
We demonstrated a positive concentration gradient of promoter I.7-specific P450arom mRNA levels in breast tissues bearing carcinomas, where mRNA levels were strikingly higher in cancer tissue and tissue proximal to carcinomas compared with adipose tissue biopsies distal to the tumor or cancer-free breast adipose tissues. We hypothesize that increased vascularity and, thus, numbers of endothelial cells with promoter I.7 activity in carcinoma tissue or breast tissues proximal to carcinomas form the basis for these observations. In fact, we and others consistently observed that P450arom protein is present in endothelial cells of the capillaries in breast cancer and surrounding adipose tissue (Sasano, H., personal communication, and our unpublished observations). A corollary to this observation is the body of literature showing P450arom expression in endothelial cells of human and murine vena cava, coronary artery, and aorta (22, 23, 24).
We do not know the functional implications of the high levels of promoter I.7-specific P450arom transcripts found in pregnant abdominal adipose tissue and normal male buttock adipose tissue. Increased vascularity in adipose tissue during pregnancy may be a potential explanation. Additionally, it is tempting to speculate that P450arom promoter I.7 activity levels exhibit a sexually dimorphic pattern in human adipose tissue. Experiments are under way to test these hypotheses.
Extensive angiogenesis, a hallmark of aggressive tumor growth, consists of endothelial cell proliferation, migration, and tube formation, yielding increased vascular density (27). We predict that endothelial cell-specific promoter I.7 activity may be the key to excessive expression of exon I.7 in breast tumors. Additionally, it is possible that paracrine factors from adjacent neoplastic or nonneoplastic cells may also contribute to the aberrant expression of exon I.7 in the vascular endothelium in breast cancer tissues. We propose that excessive levels of promoter I.7-specific P450arom mRNA in breast tumors contributing to the elevated in situ aromatase activity and local estrogen biosynthesis leading to the growth of breast cancer. In addition to its well recognized role as a potent mitogen to malignant breast epithelial cells, the ability of estrogen to stimulate angiogenesis, at least in part, by increasing vascular endothelial growth factor production leading to neovascularization, have also been demonstrated (28, 29, 30, 31). Estrogen was also shown to preserve the actin cytoarchitecture during metabolic stress, rescue endothelial cells from hypoxia-induced apoptosis, and induce the migration of endothelial cells leading to tube formation (32). Moreover, recent literature is suggestive of important physiological and pathological roles of estrogen in coronary artery disease and rheumatoid arthritis (33, 34, 35, 36). In summary, estrogen appears to play important roles in endothelial cell physiology and angiogenesis. In this regard, an endothelial cell-specific P450arom promoter giving rise to local estrogen formation with potential intracrine and paracrine effects may have significant impacts on human physiology and pathology.
 |
MATERIALS AND METHODS
|
---|
Materials
The MCF-7 (mammary gland adenocarcinoma) cell line was obtained from American Type Culture Collection (ATCC, Manassas, VA) and cultured as recommended by ATCC. The HMEC-1 (immortalized human microvascular endothelial cell line) was kindly provided by Dr. Thomas J. Lawely (37). Breast tissues were obtained from women undergoing mastectomy for breast cancer (n = 12) and reduction mammoplasty (n = 3) for macromastia. Male buttock adipose tissue sample was obtained by aspiration biopsy. A sample of abdominal adipose tissue was obtained from a pregnant woman at the time of cesarean section. Written or oral informed consent was obtained from all subjects before tissue samples were obtained. Consent forms and protocols were approved by the Institutional Review Boards of University of Texas Southwestern Medical Center at Dallas and University of Illinois at Chicago.
Fetal calf serum, media for cell cultures, LipofectAMINE PLUS reagent, Platinum Taq Polymerase, Superscript preamplification system, and oligonucleotides used as PCR primers or probes were obtained from Life Technologies, Inc. (Gaithersburg, MD). Sources of kits and reagents used in various experiments were indicated in appropriate contexts. Antibodies against GATA-1, GATA-2, and E47 used in supershift experiments were from Santa Cruz Biotechnology, Inc. (Santa Cruz, CA). Dual-Luciferase Reporter Assay System was purchased from Promega Corp. (Madison, WI). Other standard biochemicals and restriction enzymes used in the experiments were products of Sigma (St. Louis, MO), Life Technologies, Inc., or Roche Molecular Biochemicals (Indianapolis, IN). A human genomic library in EMBL3-SP6/T7 vector was procured from CLONTECH Laboratories, Inc. (Palo Alto, CA). A mammalian expression plasmid of human GATA-2 cDNA (pMT2 hGATA-2) was a gift from Dr. Stuart H. Orkin (18).
RNA Isolation, RT-PCR, and Southern Blot Analysis
Total RNA was isolated by the guanidinium thiocyanate-cesium chloride method (38). In the reverse transcription (RT) reaction 2 µg total RNA were reverse transcribed employing Superscript II enzyme (Life Technologies, Inc.) using random hexamers as the primer. For the amplification of exon I.7-specific P450arom mRNA species, the sense primer was designed as the +48/+67-bp sequence in exon I.7, whereas antisense primer was a 20-bp oligonucleotide (5'-ATTCCCATGCAGTAGCCAGG-3') complementary to the +272/+ 291-bp sequence in coding exon III of the P450arom cDNA. PCR was carried out employing the Ampli-Taq PCR system (Perkin-Elmer Cetus, Norwalk, CT) in a 20-µl reaction containing 8 pmol of each primer and 1 µl of the reverse transcription reaction as the PCR template. The PCR conditions were as the following: one cycle of 94 C for 1 min, followed by 30 cycles of 94 C for 30 sec, 55 C for 30 sec, and 72 C for 1 min with a final extension at 72 C for 5 min. The expected size of the PCR product was 259 bp. The PCR products were fractionated on a 2% agarose gel, transferred to a Zeta-Probe membrane (Bio-Rad Laboratories, Inc., Hercules, CA), and hybridized with a 18-bp radiolabeled oligonucleotide (5'-TCTGAGGTCAAGGAACAC-3') complementary to a sequence in the middle of the amplified PCR product. Additionally, we confirmed the fidelity of three amplified fragments by sequencing. As an internal control, the coding region of the housekeeping gene human glyceraldehyde-3-phosphate dehydrogenase was amplified from the reverse transcription mix and analyzed in an agarose gel containing ethidium bromide.
5'-RACE
5'-RACE was performed using the Smart RACE cDNA Amplification Kit of CLONTECH Laboratories, Inc. following the instructions provided by the manufacturer with minor modifications. Briefly, adaptor-ligated, double-strand cDNA was synthesized using ribonuclease H-reverse transcriptase, Superscript II (Life Technologies, Inc.), and a 33-bp antisense primer corresponding to the +485/+517-bp coding region (exon IV) of the human P450arom cDNA as the primer (5'-GCCTTTCTCATGCATACCGATGCACTGCAGGCC-3'). The first round of PCR was conducted using 1 µl of the adaptor-ligated 5'-RACE-ready cDNA as the template and universal primer, which is a forward primer complementary to the adaptor (component of the kit provided by CLONTECH Laboratories, Inc.) and a 21-bp antisense primer (5'-CAGGAATCTGCCGTGGGAGAT-3') complementary to the P450arom cDNA at +302/+323 bp (exon III). To increase the specificity of the reaction, a nested PCR amplification of the 1:100 diluted primary PCR reaction was performed using the nested universal primer complementary to the adaptor sequence and a 20-bp antisense primer (5'-CAGGCACGATGCTGGTGATG-3') representing the +160/+179-bp sequence (exon II) of the human P450arom mRNA. The resulting PCR products were fractionated using a 2% agarose gel. All bands within the range of 150 to 300 bp were excised and cloned into pGEM-T-Easy plasmid vector for sequence analysis. The novel exon I.7 was discovered as the 5'-end in certain portions of these clones. To map the transcription initiation site of exon I.7, exon I.7-specific PCR was used employing the RACE-ready cDNA as templates and the nested universal primer and a 32-bp antisense primer (5'-AGAGTCCCTTGCTGATTTCACCCCTTTCTCCG-3') complementary to the +76 to +101-bp sequence in exon I.7 plus the first 6 bp of the coding Exon II. The PlatinumTaq Polymerase enzyme system supplied by Life Technologies, Inc. was employed in all PCR experiments.
Colony Hybridization
Bacterial colonies (a total of 270) representing a library of 5'-RACE clones of P450arom transcripts expressed in breast cancer tissue specimens were replica-plated. Each replica plate was hybridized with an 18-bp radiolabeled oligonucleotide (5'-TCTGAGGTCAAGGAACAC-3') complementary to the sequence within the coding exon II of the P450arom cDNA or a 21-bp oligonucleotide complementary to the +27/+47-bp sequence (5'-GATGAAGTGTTGTATCTTAGG-3') of exon I.7.
Mapping of Exon I.7 to the Human Genome
Nucleotide sequence information of the exon I.7 (101 bp) and the human full-length aromatase cDNA (GenBank accession no.: M22246) were subjected to the BLAST homology search (39) against the Human Genome Resource database of the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/genome/seq/HsBlast.html).
Cloning of the 5'-Flanking Region of Exon I.7
The 5'-upstream region of the newly identified untranslated first exon I.7 in the human CYP19 gene was cloned from a EMBL3-SP6/T7 human genomic library constructed from human placental DNA partially digested with Sau3AI restriction enzyme (CLONTECH Laboratories, Inc.). For library screening employing oligonucleotides as probes, the protocol recommended by the manufacturer was followed (Lambda Library User Manual, PT1010-1, CLONTECH Laboratories, Inc.). Briefly, the human EMBL3 genomic library was sequentially screened by plaque lift hybridization using two oligonucleotide probes. For the first screening, we used a 21-bp oligonucleotide complementary to the +27/+47-bp sequence in exon I.7. The positive plaques were then rescreened with a second 20-bp oligonucleotide complementary to the +48/+67-bp sequence in exon I.7. The sequences of these two oligonucleotides were 5'-GATGAAGTGTTGTATCTTAGG-3' and 5'-GGCTCCATCTACAAGGATGA-3', respectively. Using purified DNA of a positive phage clone as the template, we amplified a region (
2.4 kb) by PCR using a primer pair complementary to +51/+81-bp sequence of exon I.7 and the T7 primer flanking the EMBL-SP6/T7 cloning site. For the PCR we employed the Expand Long Template PCR System (Roche Molecular Biochemicals) and 200 ng purified phage DNA (as template). The PCR conditions were as follows: one cycle of denaturation at 95 C for 3 min, followed by 25 cycles of 15 sec at 94 C, 1 min at 42 C, 12 min at 68 C, followed by a final extension at 68 C for 10 min. The gel-purified genomic fragment (
2.4 kb) was subsequently cloned into the pGEM-Teasy vector (Promega Corp.) employing the T-tailing method. The cloned fragment was sequenced to its entirety by the primer walking method using ABI PRISM BigDye Terminator Cycle sequencing kit and ABI PRISM 377 automated sequencer (PE Applied Biosystems, Foster City, CA).
Generation of Luciferase Reporter Constructs
Deletion of specific regions of the 2.4-kb cloned DNA fragment was accomplished through the use of available restriction sites and subcloning into compatible sites of the luciferase reporter plasmid pGL3-basic. In some cases, PCR was used to create flanking restriction sites for cloning purposes.
Mutations were performed using the QuikChange Site-Directed Mutagenesis Kit from Stratagene (La Jolla, CA), following the protocol provided by the manufacturer. A 380-bp fragment corresponding to the -299- to +81-bp region cloned into the pGL3-Basic vector was used as a template for PCR amplification. Mutations were introduced into putative GATA binding sites in this -299/+81-bp construct. For the -146/-141-bp GATA site, the sequence GAGATAAT was changed to GActcgAT, and for the -196/-191 GATA site, GTGATAGT was changed to GTctcgGT. The mutated bases are shown in lowercase. The mutations were confirmed by direct sequencing, and the 380-bp fragment (-299/+81 bp) containing the desired mutations was isolated from the original pGL3 construct and cloned into a separate MluI-BglII-digested empty pGL3-Basic vector.
Plasmids used in transfection experiments were purified using an EndoFree Plasmid Isolation Kit (QIAGEN, Valencia, CA), and purity was verified by A260/A280 absorbance and agarose gel electrophoresis. In all constructs, orientation of inserts was verified by direct sequencing.
Transfection of HMEC-1 and MCF-7 Cells with Luciferase Reporter Constructs
HMEC-1 cells were grown in endothelial basal medium MCDB-131 (Life Technologies, Inc.) containing appropriate supplements as described previously (16). MCF-7 cells were grown as recommended by ATCC. The day before transfection, cells were plated into six-well tissue culture dishes at a density such that the cells reached 7080% confluence by the time of transfection. Transfections were performed using LipofectAMINE PLUS reagent (Life Technologies, Inc.) following the protocol provided by the manufacturer. One microgram of the -299/+81-bp luciferase reporter construct or equimolar amounts of various deletion constructs were used. As an internal control, 25 ng of pRL (Renilla luciferase)-Null plasmid (Promega Corp.) was cotransfected. Three hours after transfection, the transfection medium was removed by aspiration, 3 ml of complete medium (containing serum and antibiotics) were added, and the plates were returned to the incubator. In the case of HMEC-1 cells, transfections were carried out using Lipofectin reagent (Life Technologies, Inc.) for 5 h. At 48 h post transfection, medium was removed and wells were rinsed with PBS to remove detached cells and residual growth medium. Next, 250 µl of 1x passive lysis buffer provided in the Dual-Luciferase Reporter Assay System (Promega Corp.) were added to each well. Cells were dispersed by scraping with a disposable plastic cell lifter and subjected to one cycle of freeze (-80 C)/thaw (room temperature) to ensure complete cell lysis. Ten microliters of the supernatant were used to assay luciferase activities. Firefly and Renilla luciferase activities were sequentially measured using the Dual-Luciferase Reporter Assay System (Promega Corp.) following the manufacturers instructions. Luciferase activities were determined using a Lumat LB 9507 Luminometer equipped with a variable auto dual injector (Berthold Technologies, Bad Wildbad, Germany) The protocol included a 1-sec predelay followed by a 20-sec measuring period. The Renilla luciferase activity, expressed from the pRL-Null vector, provided an internal control to monitor transfection efficiency. Firefly luciferase activities were normalized based on the Renilla luciferase activity in each well. Statistical analysis and graphic representation of the results were accomplished using the SigmaPlot 5 Program (SPSS, Inc., Chicago, IL).
Preparation of Nuclear Extracts from HMEC-1
Nuclear extracts were prepared from HMEC-1 cells using the NE-PER Nuclear and cytoplasmic extraction reagent kit supplied by Pierce Chemical Co. Protein content of the prepared nuclear extract was determined by the BCA-200 Protein Assay Kit supplied by Pierce Chemical Co.
DNase I Footprinting
For in vitro DNase I footprinting, a SacIXbaI fragment containing the -287/-92-bp region of promoter I.7 was gel purified and was labeled by filling in with Klenow fragment using
[32P]dCTP. Because SacI restriction digestion creates a 3'-protruding end, Klenow filling is restricted to the XbaI-digested end of the fragment. The labeled fragment was purified using a QIAquick PCR Purification Kit (QIAGEN), and radioactivity was quantitated by scintillation counting. The footprinting was performed using the SureTrack Footprinting Kit (Amersham Pharmacia Biotech, Arlington Heights, IL) and following the protocol provided by the manufacturer, except that the binding reaction was conducted in a 20-µl reaction volume, and the volume was adjusted to 50 µl by adding 30 µl of 1x binding buffer before proceeding with DNase I digestion. Samples were analyzed on a 6% sequencing gel. A sequencing ladder was obtained by loading a sample of the probe that had been treated with formic acid and piperidine to cleave at G and A residues.
EMSA
A 30-mer oligonucleotide (representing the -209- to -180-bp sequence), which harbors the -196/-191-bp GATA motif, was radiolabeled using
[32P]ATP and T4 polynucleotide kinase. The labeled oligonucleotide was purified using a QIAquick Nucleotide purification Kit (QIAGEN), and radioactivity was quantified by scintillation counting. Approximately 30,000 cpm of labeled probe and 250 ng of HMEC-1 nuclear extract were mixed with indicated amounts of oligonucleotide competitors or antibodies (supershift) in a gel shift binding buffer supplied by Promega Corp. The reaction mix had been incubated on ice for 1 h before adding the labeled probe (
30,000 cpm), after which the incubation was continued for 20 min at 30 C. Samples were then analyzed on a 5% nondenaturing acrylamide gel using 0.5x Tris-borate-EDTA as the running buffer. The dried gel was exposed to Biomax MR autoradiographic film (Eastman Kodak Co., Rochester, NY).
 |
ACKNOWLEDGMENTS
|
---|
We thank Dr. Stuart H. Orkin of Harvard Medical School (Boston, MA) for providing the expression plasmid pMT2-hGATA-2, Dr. David B. Wilson of Washington University School of Medicine (St. Louis, MO) for pMT2-empty vector, Dr. Thomas J. Lawley of Centers for Disease Control, (Atlanta, GA) for HMEC-1 cells, and Dr. C. Tiruppathi of the Department of Pharmacology, University of Illinois at Chicago (Chicago, IL) for timely help.
 |
FOOTNOTES
|
---|
This work was supported in part by United States Army Medical Research and Material Command Grant DAMD17-97-17025 and National Cancer Institute Grant CA67167 (to S.E.B.).
1 Present address: Department of Pathology, Duke University Medical Center 3712, Durham, North Carolina 27710. 
Abbreviations: BLAST, Basic local alignment search tool; DNase I, deoxyribonuclease I; HMEC-1, human microvascular endothelial cells; P450arom, aromatase cytochrome P450; RACE, rapid amplification of cDNA ends; UTR, untranslatory region.
2 The nucleotide sequences reported in this paper have been submitted to the GenBank/EBI data bank with accession nos. AF419337 and AF419338. 
Received for publication March 28, 2002.
Accepted for publication July 5, 2002.
 |
REFERENCES
|
---|
- ONeill JS, Elton RA, Miller WR 1988 Aromatase activity in adipose tissue from breast quadrants: a link with tumour site. Br Med J 296:741743[Medline]
- Yue W, Wang JP, Hamilton CJ, Demers, LM, Santen RJ 1998 In situ aromatization enhances breast tumor estradiol levels and cellular proliferation. Cancer Res 58:927932[Abstract]
- Bulun SE, Price TM, Aitken J, Mahendroo MS, Simpson ER 1993 A link between breast cancer and local estrogen biosynthesis suggested by quantification of breast adipose tissue aromatase cytochrome P450 transcripts using competitive polymerase chain reaction after reverse transcription. J Clin Endocrinol Metab 77:16221628[Abstract]
- Lu Q, Nakmura J, Savinov A, Yue W, Weisz J, Dabbs DJ, Wolz G, Brodie A 1996 Expression of aromatase protein and messenger ribonucleic acid in tumor epithelial cells and evidence of functional significance of locally produced estrogen in human breast cancers. Endocrinology 13:7306173068
- Brodie A M, Njar VC 2000 Aromatase inhibitors and their application in breast cancer treatment. Steroids 65:171179[CrossRef][Medline]
- Santen RJ, Harvey HA 1999 Use of aromatase inhibitors in breast carcinoma. Endocr Relat Cancer 6:7592[Abstract/Free Full Text]
- Bonneterre J, Buzdar A, Nabholtz JM, Robertson JF, Thurlimann B, von Euler M, Sahmoud T, Webster A, Steinberg M 2001 Anastrozole is superior to tamoxifen as first-line therapy in hormone receptor positive advanced breast carcinoma. Cancer 92:22472258[CrossRef][Medline]
- Sebastian S, Bulun SE 2001 A highly complex organization of the regulatory region of the human CYP19 (aromatase) gene revealed by the Human Genome Project. J Clin Endocrinol Metab 86:46004602[Free Full Text]
- Simpson ER, Michael MD, Agarwal VR, Hinshelwood MM, Bulun SE, Zhao Y 1997 Cytochromes P450 11: expression of the CYP19 (aromatase) gene: an unusual case of alternative promoter usage. FASEB J 11:2936[Abstract/Free Full Text]
- Simpson ER, Clyne C, Rubin G, Boon WC, Robertson K, Britt K, Speed C, Jones M 2002 Aromatasea brief overview. Annu Rev Physiol 64:93127[CrossRef][Medline]
- Bulun SE, Mahendroo MS, Simpson ER 1994 Aromatase gene expression in adipose tissue: relationship to breast cancer. Steroid Biochem Mol Biol 49:319326
- Agarwal VR, Bulun SE, Leitch M, Rohrich R, Simpson ER 1996 Use of alternative promoters to express the aromatase cytochrome P450 (CYP19) gene in breast adipose tissues of cancer-free and breast cancer patients. J Clin Endocrinol Metab 81:38433849[Abstract]
- Zhou C, Zhou D, Esteban J, Murai J, Siiteri PK, Wilczynski S, Chen S 1996 Aromatase gene expression and its exon I usage in human breast tumors. Detection of aromatase messenger RNA by reverse transcription-polymerase chain reaction. J Steroid Biochem Mol Biol 59:163171[CrossRef][Medline]
- Harada NJ 1997 Aberrant expression of aromatase in breast cancer tissues. Steroid Biochem Mol Biol 61:175184[CrossRef]
- Bucher P 1990 Weight matrix descriptions of four eukaryotic RNA polymerase II promoter elements derived from 502 unrelated promoter sequences. J Mol Biol 212:563578[Medline]
- Merika M, Orkin SH 1993 DNA-binding specificity of GATA family transcription factors. Mol Cell Biol 13:39994010[Abstract]
- Zhou J, Gurates B, Yang S, Sebastian S Bulun SE 2001 Malignant breast epithelial cells stimulate aromatase expression via promoter II in human adipose fibroblasts: an epithelial-stromal interaction in breast tumors mediated by CCAAT/enhancer binding protein ß. Cancer Res 61:23282334[Abstract/Free Full Text]
- Dorfman DM, Wilson DB, Bruns GA, Orkin SH 1992 Human transcription factor GATA-2. Evidence for regulation of preproendothelin-1 gene expression in endothelial cells. J Biol Chem 267:12791285[Abstract/Free Full Text]
- Lee ME, Temizer DH, Clifford JA, Quertermous T 1991 Cloning of the GATA-binding protein that regulates endothelin-1 gene expression in endothelial cells. J Biol Chem 66:1618816192
- Gumina RJ, Kirschbaum NE, Piotrowski K, Newman PJ 1997 Characterization of the human platelet/endothelial cell adhesion molecule-1 promoter: identification of a GATA-2 binding element required for optimal transcriptional activity. Blood 89:12601269[Abstract/Free Full Text]
- Zhang R, Min W Sessa, WC 1995 Functional analysis of the human endothelial nitric oxide synthase promoter. Sp1 and GATA factors are necessary for basal transcription in endothelial cells. J Biol Chem 270:1532015326[Abstract/Free Full Text]
- Sasano H, Murakami H, Shizawa S, Satomi S, Nagura H, Harada N 1999 Aromatase and sex steroid receptors in human vena cava. Endocr J 46:233242[Medline]
- Diano S, Horvath TL, Mor G, Register T, Adams M, Harada N, Naftolin F 1999 Aromatase and estrogen receptor immunoreactivity in the coronary arteries of monkeys and human subjects. Menopause 6:2128[Medline]
- Nathan L, Shi W, Dinh H, Mukherjee TK, Wang X, Lusis AJ, Chaudhuri G 2001 Testosterone inhibits early atherogenesis by conversion to estradiol: critical role of aromatase. Proc Natl Acad Sci USA 98:35893593[Abstract/Free Full Text]
- Weiss MJ, Orkin SH 1995 GATA transcription factors: key regulators of hematopoiesis. Exp Hematol 23:99107[Medline]
- Simon MC 1995 Gotta have GATA. Nat Genet 11:911[Medline]
- Weidner N, Semple JP, Welch WR, Folkman J 1991 Tumor angiogenesis and metastasis-correlation in invasive breast carcinoma. N Engl J Med 324:18[Abstract]
- Losordo DW, Isner JM 2001 Estrogen and angiogenesis. Arterioscler Thromb Vasc Biol 21:612[Abstract/Free Full Text]
- Shifren JL, Tseng JF, Zaloudek CJ, Ryan IP, Meng YG, Ferrara N, Jaffe RB, Taylor RN 1996 Ovarian steroid regulation of vascular endothelial growth factor in the human endometrium: implications for angiogenesis during the menstrual cycle and in the pathogenesis of endometriosis. J Clin Endocrinol Metab 81:31123118[Abstract]
- Mueller MD, Vigne JL, Minchenko A, Lebovic DI, Leitman DC, Taylor RN 2000 Regulation of vascular endothelial growth factor (VEGF) gene transcription by estrogen receptors
and ß. Proc Natl Acad Sci USA 97:1097210977[Abstract/Free Full Text]
- Concina P, Sordello S, Barbacanne MA, Elhage R, Pieraggi MT, Fournial G, Plouet J, Bayard F Arnal JF 2000 The mitogenic effect of 17ß-estradiol on in vitro endothelial cell proliferation and on in vivo reendothelialization are both dependent on vascular endothelial growth factor. J Vasc Res 37:202208[CrossRef][Medline]
- Razandi M, Pedram A, Levin ER 2000 Estrogen signals to the preservation of endothelial cell form and function. J Biol Chem 275:3854038546[Abstract/Free Full Text]
- Farhat MY, Lavigne MC, Ramwell PW 1996 The vascular protective effects of estrogen. FASEB J 10:615624[Abstract/Free Full Text]
- Russell KS, Haynes MP, Sinha, D, Clerisme E, Bender JR 2000 Human vascular endothelial cells contain membrane binding sites for estradiol, which mediate rapid intracellular signaling. Proc Natl Acad Sci USA 97:59305935[Abstract/Free Full Text]
- Kim-Schulze S, McGowan KA, Hubchak SC, Cid MC, Martin MB, Kleinman, HK, Greene GL, Schnaper HW 1996 Expression of an estrogen receptor by human coronary artery and umbilical vein endothelial cells. Circulation 94:14021407[Abstract/Free Full Text]
- Venkov CD, Rankin AB, Vaughan DE 1996 Identification of authentic estrogen receptor in cultured endothelial cells. A potential mechanism for steroid hormone regulation of endothelial function. Circulation 94:727733[Abstract/Free Full Text]
- Ades EW, Candal FJ, Swerlick RA, George VG, Summers S, Bosse DC, Lawley TJ 1992 HMEC-1: establishment of an immortalized human microvascular endothelial cell line. J Invest Dermatol 99:683690[Abstract]
- Chirgwin JM, Przybyla AE, MacDonald RJ, Rutter WJ 1979 Isolation of biologically active ribonucleic acid from sources enriched in ribonuclease. Biochemistry 18:52945299[Medline]
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ 1990 Basic local alignment search tool. J Mol Biol 215:403410[CrossRef][Medline]