| Nucleic Acids Research | Pages |
Transcription analysis of the telomeric repeat-specific retrotransposons TRAS1 and SART1 of the silkworm Bombyx mori
Introduction
Materials And Methods
Nucleotide sequence accession number
Northern hybridization
5[prime]-RACE analysis
Inverse PCR
Promoter assay
Results
TRAS1 and SART1 were transcribed as full-length units in various tissues and in the BmN4 cell line
Most genomic copies of TRAS1 and SART1 contained the transcription initiation sequences
TRAS1 possesses a promoter in its 5[prime]-UTR
Influence of the first 40 bp conserved in most insect retrotransposons on transcription of TRAS1
The neighboring telomeric repeats inhibit transcription of TRAS1
Discussion
Acknowledgements
References
Transcription analysis of the telomeric repeat-specific retrotransposons TRAS1 and SART1 of the silkworm Bombyx mori
ABSTRACT
INTRODUCTION
Most eukaryotic telomeres consist of direct repeats of oligonucleotides (telomeric repeats) which are elongated by a reverse trans-criptase named telomerase (1). Some insects, however, lack the telomeric repeats and offer interesting models for the study of the function of telomeres. These insects have alternative mechanisms to maintain their telomeres. The telomeres of Drosophila melanogaster are elongated by addition, approximately once every 100 generations per one chromosomal terminus, of the retro-transposable elements HeT-A and TART (2,3). The telomeres of the dipteran Chironomus papallidivittatus and Anopheles gambiae are likely to be maintained by recombinational events (4,5). These telomerase-independent telomere maintenance mechanisms may represent general back-up for telomerase in other eukaryotes, including yeast, mouse and human (6-10).
The telomere of the silkworm Bombyx mori has a structure intermediate between canonical repeat-type and Drosophila-type telomeres (Fig.
Figure 1. The telomeric repeat-specific retrotransposons in B.mori. (A) Schematic structure of the Bombyx chromosome. Black triangles represent the telomeric repeats, (TTAGG)n. Hatched boxes indicate the telomeric repeat-specific retrotransposons. The extreme ends (~6-8 kb) consist only of the telomeric repeats, while a large number of sequence-specific retrotransposons are mixed with the subtelomeric repeats. 50-100 copies of the retrotransposons may exist at one terminus. (B) Sequence specificity of the TRAS/SART families. The position at which the members of the TRAS or SART family insert into the telomeric repeats are shown with vertical arrows, based on the 5[prime]-RACE analysis in Figure 3. The ORFs and UTRs are depicted as open boxes and horizontal lines, respectively, above which nucleotide positions are indicated. Both large families comprise an endonuclease domain (EN), a reverse transcriptase domain (RT), zinc finger-like motifs (vertical lines near the C-termini of both ORFs) and a poly(A) tail. The TRAS families exist with their poly(A) tails facing toward the termini and the SART families are in the opposite orientation. Transcription is an important regulatory step for non-LTR retrotransposons because their RNA serves not only as the messenger for translation but also as the template for reverse transcription. Although TRAS1 and SART1 encode fully conserved ORFs, it is unknown if they are actually transcribed in vivo and if a few master copies or most of the genomic copies are transcribed. It is also unknown whether the sequence-specific, non-LTR elements have a unique transcription mechanism. Many non-sequence-specific, non-LTR elements such as Drosophila Jockey (18), I (19), F (20) and Doc (21) and human LINE-1 (22,23) have been shown to contain promoter activities in the 5[prime]-untranslated regions (UTRs). It seems rational for those non-LTR elements that transpose to many different genomic locations to have promoters downstream of the transcription initiation sites in order to retain their own promoter sequences after transcription and reverse transcription. TRAS1 and SART1, however, insert only at specific positions within the telomeric repeats and thus it is unclear if they have downstream promoters, as do other elements. There is also a possibility that their promoters are dependent on the existence of the telomeric repeats. It is interesting in this regard that transcription of the Drosophila HeT-A elements that exist as tandem repeats at telomeres is promoted by the neighboring copies (24). In this report, we focus on the transcription of TRAS1 and SART1, showing that TRAS1 and SART1 are transcribed in many tissues and in the BmN4 cell line. We demonstrate that the two elements are transcribed from their own 5[prime]-ends of the retrotransposon units without the telomeric repeats at the 5[prime]-ends. Using transient reporter assays, we identified promoter activity in the upstream region of the 5[prime]-UTR of TRAS1 in the BmN4 cell line. We discuss whether the elements are transcribed in vivo from most of the copies within the telomeric repeats or from a few master copies.
MATERIALS AND METHODS
Nucleotide sequence accession number
Nucleotide positions indicated in the text are referred to based on D38414 (TRAS1) and D85594 (SART1). In this paper, however, we redefine the transcription initiation sites as the +1 positions which were determined by 5[prime]-rapid amplification of cDNA ends (5[prime]-RACE).
Northern hybridization
Poly(A)+ RNAs were isolated from testes, ovaries, fat bodies, malpighian tubules and posterior silk glands of fifth instar larvae and from BmN4 cultured cells of B.mori with Micro-Fast Track[trade] (Invitrogen). Aliquots of 5 µg of poly(A)+ RNA per lane were electrophoresed at 5 V/cm on 18% formaldehyde, 20 mM MOPS (3-morpholinopropanesulfonic acid, pH 7.0), 5 mM sodium acetate, 1 mM EDTA, 0.9% agarose gels and blotted onto nylon membranes (Hybond-N; Amersham) in 10× SSC (1.5 M NaCl, 0.15 M sodium citrate). After prehybridization, the membranes were hybridized with ~106 c.p.m. of each probe at 42°C overnight in 40% formamide, 10× Denharts solution (0.2% each of BSA, Ficoll and polyvinylpyrrolidone), 5× SSC, 250 µg/ml salmon sperm DNA, 50 mM NaPO4 (pH 7.0), 10% dextran sulfate. The double-strand DNA probes shown in Figure
Table 1.
5[prime]-RACE analysis
5[prime]-RACE was conducted through self-ligation of synthesized cDNA (25) using the 5[prime]-Full Race Core Set (Takara). An aliquot of 20 ng of BmN4 poly(A)+ RNA was reverse transcribed with 5 U of AMV reverse transcriptase at 42°C for 60 min and then at 50°C for 30 min using 200 pmol of TRAS1 A307P primer or SART1 A287P primer, whose 5[prime]-ends were phosphorylated for later ligations. The cDNA was treated with 30 U of RNase H at 32°C and ethanol precipitated. The phosphorylated 5[prime]-end and the 3[prime]-end of the single-strand cDNA were self-ligated with 40 U of T4 RNA ligase. The ligated circular cDNA was amplified in two steps by nested PCR with LA Taq polymerase (Takara). The primer pairs for the nested PCR were as follows (Fig.
Inverse PCR
Inverse PCR was carried out as described by Triglia et al. (26). Genomic DNA was extracted from the silk glands of fifth instar larvae as described previously (13). Aliquots of 15 ng of the genomic DNA were digested with a four-base cutter, MspI or AluI, which cut restriction sites at the +323 or +544 positions in TRAS1 and SART1, respectively. The digested DNA was circularized through self-ligation using the DNA Ligation Kit V.2 (Takara). The DNA was amplified in two steps by PCR with Ex Taq polymerase (Takara). The primer sets for the inverse PCR were the same as those used for the 5[prime]-RACE PCR: TRAS1 first PCR, A142 and S203; TRAS1 second PCR, A79 and S245; SART1 first PCR, A130 and S141; SART1 second PCR, A80 and S177. The PCR products were cloned into the pGEM T Easy vectors (Promega) and sequenced as described above.
Promoter assay
The in vivo promoter assay was conducted by transient trans-fection of each experimental DNA cloned in the pGL3-Enhancer vector containing a firefly luciferase gene as reporter (the Dual-Luciferase[trade] Reporter system; Promega) into the BmN4 cell line. To normalize expression of the firefly luciferase, the pRL vector containing the Drosophila HSP70 promoter, which provides constitutive expression of Renilla luciferase in BmN4 cells, was co-transfected with each experimental DNA in the firefly luciferase vector. The Drosophila HSP70 promoter was a gift from Dr H. Maekawa and was subcloned into the HindIII site of the pRL-null vector (HSP70-pRL vector). In order to construct experimental firefly luciferase vectors, each DNA was amplified by PCR with LA Taq polymerase (Takara), with a pair of primers designed to make sequences for a NheI and a HindIII site at each end of the PCR product. The PCR products digested with NheI and HindIII were cloned into the pGL3-Enhancer vector. Each construct was sequenced and turned out to be free from PCR errors. The experimental DNAs were transfected into the BmN4 cell line by a liposome-mediated method using the Tfx[trade] 20 reagent (Promega). To ~1 × 105 BmN4 cells adherent to a well in a 96-well plate was added 40 µl of serum-free TC100 medium (Gibco-BRL) which contains 0.75 µg of each pGL3 construct, 0.01 µg of the HSP70-pRL vector and 2.3 µl of Tfx[trade] 20. Four hours later, 200 µl of TC100 medium with 10% fetal bovine serum was added. The luciferase assay was conducted 72 h later according to the manufacturers instructions using a Luminescence Reader BLR-201 (ALOKA).
Figure 2. Northern hybridization. Aliquots of 5 µg of poly(A)+ RNA from testis (T), ovary (O), fat body (F), malpighian tubule (M) and posterior silk gland (S) of fifth larva and from the BmN4 cultured cells (B) were hybridized with a TRAS1 probe (A), a SART1 probe (B) and an EF1-[alpha] probe (C). (D) The regions for the double-strand DNA probes used for TRAS1 and SART1 are shown by arrows (both ORFs and the 3[prime]-UTR). The band around 2 kb in (A) is probably rRNA and was never seen reproducibly. Note that hybridization of RNAs from the BmN4 cells with the TRAS1 probe was conducted in a separate experiment and thus the expression level cannot be compared to those in other tissues. To determine if TRAS1 and SART1 are transcribed in vivo, we extracted 5 µg of poly(A)+ RNAs from several tissues of fifth instar larvae and from BmN4 cells and carried out northern hybridization. We hybridized a single filter with a TRAS1 probe (Fig. Many non-sequence-specific, non-LTR elements are transcribed from their 5[prime]-ends by downstream promoters (18-23). It is unclear, however, if this is the case in sequence-specific non-LTR elements. In the R1 element, in fact, a low level of read-through transcript from 28S rDNA was found (27). To investigate the transcription initiation mechanism of the sequence-specific elements, we amplified the 5[prime]-ends of TRAS1 and SART1 transcripts (85 nt) from BmN4 cells by 5[prime]-RACE. Although we used RNAs from BmN4 cells as a template, the same result may be expected from other tissues because there seemed to be no apparent RNA size differences among tissues and BmN4 cells (Fig. Figure 3. Sequence comparison of the 5[prime]-ends of the transcripts of TRAS1 and SART1 with the counterparts in their genomic copies using 5[prime]-RACE (25) and inverse PCR (26). The position of the 5[prime]-phosphorylated primers used for reverse transcription in the 5[prime]-RACE and the primers for the nested PCRs in the 5[prime]-RACE and the inverse PCR are shown in Figure 1B. All the 5[prime]-RACE clones started at the same position (with one base uncertainty in the case of TRAS1) and contained almost the same sequences through 85 nt. Most of their genomic copies (9/10 and 13/18 for TRAS1 and SART1, respectively) which were amplified by inverse PCR contained the transcription initiation sites. Because the 5[prime]-ends of non-LTR retrotransposons are usually variable, probably due to incomplete reverse transcription and aberrant integration processes (28), only a few copies of TRAS1 and SART1 may contain their transcription start sequences. To examine what proportions of their genomic copies retain the transcription initiation sequence, we amplified DNA sequences around their 5[prime]-ends by inverse PCR (26). All the 5[prime]-end clones of TRAS1 and SART1 obtained by inverse PCR were adjacent to the telomeric repeats and retained their same sequence specificities as shown in the genomic clones (Fig. The discovery that TRAS1 and SART1 contain their own transcription initiation sites prompted us to examine whether their promoters lie downstream of the initiation sites, as do many non-sequence-specific, non-LTR elements. We subcloned their 5[prime]-UTRs into the pGL3-Enhancer vector, which contains the firefly luciferase gene as a reporter gene (Promega), and transiently transfected them into BmN4 cells. BmN4 cells are biologically relevant with regard to transcription of TRAS1 and SART1 as we proved by northern hybridization (Fig. The luciferase assays showed that TRAS1 contained its own promoter activity in its 5[prime]-UTR. The luciferase activities obtained from a series of 3[prime]-deleted constructs of the 5[prime]-UTR are shown in Figure Figure 4. Transient luciferase activity of 3[prime]-deletion derivatives of the TRAS1 5[prime]-UTR in BmN4 cells. (A) Structure of the -2/2432 construct which contains all the 5[prime]-UTR of TRAS1. The linker sequence and the firefly luciferase gene (Luc) are depicted as a gray box and a white box, respectively. This construct contains cytosine dinucleotides (CC) upstream of the transcription initiation sequence, GAGTT. The positions of ATG triplets that reduce translation of luciferase are indicated by vertical lines. (B) Luciferase activity of 3[prime]-deletion derivatives of the TRAS1 5[prime]-UTR. The activity for each experimental pGL3-Enhancer vector was normalized to those of the HSP70-containing pRL vector. The activity for -2/521 (maximal) is designated as 100% and values of the mean ± standard deviation for three independent experiments per construct are shown. HSP70 denotes the Drosophila HSP70 promoter that was used as a positive control. Activity for the pGL3-Enhancer vector alone is also shown below as a negative control. We tried to detect promoter activity in the SART1 5[prime]-UTR but the constructs, which contained all or parts of the SART1 5[prime]-UTR, did not show greater activity than the pGL3-Enhancer vector alone (data not shown). We then examined the possibility that transcription of SART1 is promoted by the 3[prime]-UTR of a neighboring SART1 copy, as is the case with the Drosophila HeT-A element (24), considering that at least half of all SART1 elements exist as tandem arrays interrupted by short (<100 bp) telomeric repeats (13). The constructs which contain the 3[prime]-UTR upstream of the telomeric repeats and the 5[prime]-UTR of SART1 did not, however, show significant activity (data not shown). Most non-LTR retrotransposons that have been well characterized in insects have downstream promoters in their 5[prime]-UTR. These elements have in common a TATA-less RNA polymerase II promoter, which have also been found in some LTR retrotrans-posons in insects (29), many Drosophila homeotic genes such as Antennapedia (30) and Engrailed (31) and the mammalian terminal deoxynucleotidyltransferase (TdT) gene (32), and so on. These retrotransposons share two conserved sites in the 5[prime]-terminal regions (29). The initiator [the (C/G)A(C/G/T)T motif] is located around the transcription initiation site and the CGT(G/T) motif is found at position +31 or +32 (Fig. Figure 5. Sequence comparison of the first 40 bp of non-LTR retrotransposons in insects for which transcription initiation sites have been determined. The initiator at the 5[prime]-end and the downstream CGT(G/T) motif are underlined. TRAS1 contains an initiator-like sequence while SART1 does not. To examine the importance of an initiator on transcription of TRAS1, we created a series of 5[prime]-deleted TRAS1 5[prime]-UTR constructs and performed luciferase assays (Fig. Figure 6. Luciferase activity of 5[prime]-variation derivatives of the TRAS1 5[prime]-UTR in BmN4 cells. The activity for +1/551 (maximal) is taken as 100% and values of the mean ± standard deviation for three independent experiments per construct are shown. The telomeric repeats are shown as (CCTAA)n. TRAS1 contained an internal promoter in its 5[prime]-UTR as described above. There is a possibility, however, that the neighboring telomeric repeats also have some influence on transcription of TRAS1. To examine this possibility, we constructed a series of plasmids that contained various short lengths of the telomeric repeats upstream of the first 551 bp of the TRAS1 element and performed luciferase assays. As shown in Figure Previous structural analyses, which showed the existence of undisrupted ORFs, a few truncated copies and sequence conser-vation among genomic copies, suggested that most copies of TRAS1 and SART1 in the genome have the ability to transpose via RNA (12,13). In this paper, we provide additional evidence for this idea by showing that TRAS1 and SART1 are actually expressed in many tissues and a cultured cell line as approximately full-length RNAs, which can serve not only as messengers for translation but as templates for target-primed reverse transcription. Their expression in testis and ovary tissues and in BmN4 cells suggests that they have been surviving by retrotransposition through vertical transmission and/or that they play some important role in telomere maintenance. Although most non-LTR retrotrans-posons show various lengths of transcripts due to 5[prime]-truncation and read-through transcription from neighboring genes, only single sizes of the TRAS1/SART1 RNAs were found. This is probably because only a few of their genomic copies are truncated and because the TRAS1 and SART1 copies are not transcribed from adjoining genes, being situated exclusively in the telomeric repeats. TRAS1 and SART1 are the first examples of sequence-specific retrotransposons that are actively transcribed in vivo. Although sequence-specific retrotransposons (R1 and R2) may be transcribed as read-through RNAs from neighboring rRNA genes (27,33), TRAS1 and SART1 were expressed as unit length RNAs. This may reflect differences in genomic localization between the rRNA genes and non-transcribed telomeric repeats. Major questions regarding transcription of TRAS1 and SART1 are concerned with whether the copies within the telomeric repeats are transcribed and whether they possess internal promoters. The data presented here suggest that most TRAS1 elements within the telomeric repeats are weakly transcribed in vivo by internal promoters. First, the TRAS1 RNA that was weakly detected on the northern blot corresponds in size (~8 kb) to the previously cloned TRAS1 element (TRAS1-[lambda]B1) that was adjacent to the telomeric repeats (12). Second, 5[prime]-RACE analysis showed that the 5[prime]-end of the TRAS1 RNA started precisely from the 5[prime]-end of TRAS1-[lambda]B1 and that the nucleotide sequence of the RNA was completely identical over 85 nt to that of a genomic TRAS1 element. Third, we have never found, by inverse PCR, a TRAS1 element that was inserted outside the telomeric repeats (Fig. Internal promoters have already been identified for some non-sequence-specific, non-LTR retrotransposons such as Jockey (18), I (19), F (20) and Doc (21) from D.melanogaster, the human LINE-1 element (22,23) and DRE (34), which integrates upstream of Dictyostelium discoideum tRNA genes with position specificity but not sequence specificity. This study showed that a sequence-specific, non-LTR element TRAS1 also comprised an internal promoter, as do many non-sequence-specific elements. Transcription of TRAS1 occurred independently of the telomeric repeats, which rather inhibited transcription. Non-LTR retrotrans-posons may be autonomous elements with regard to transcription and be independent of surrounding environments. This result does not necessarily mean, however, that all sequence-specific elements are transcribed by internal promoters. For example, R1 and R2 may be transcribed as one unit with rRNA (27,33). TRAS1 possesses an initiator-like sequence at its 5[prime]-end as do the retrotransposons I, Jockey, F and Doc. In Drosophila retrotransposons I, F and Doc, the first 40 bp segment contains all the promoter activity in cultured Drosophila S2 cells (19-21). Further, the sequences downstream of the first 40 bp regulate the tissue specificity in the case of I (35) and Fex (36). On the other hand, eliminating the first 43 bp of TRAS1 did not completely eliminate the promoter activity (Fig. Drosophila depends for maintenance of its telomeres upon retrotransposons (2,3), which have an intimate evolutionary association with telomerase (37). The fact that TRAS1 and SART1 are actively transcribed retrotransposons suggests the importance of retrotransposition of the TRAS/SART families on the evolution of insect telomeres. One can imagine that some insects lost the ability to produce telomerase at some point during their evolution, but maintained telomere length by insertion of TRAS/SART-like elements into the telomeric repeats. Various mechanisms such as addition of retrotransposons onto the very ends of telomeres and recombinational events may have replaced telomerase activity. Since so far we have not succeeded in detecting telomerase activity (Y.Sasaki and H.Fujiwara, unpublished data), telomerase might have been lost and an alternative mechanism might have arisen in B.mori, although some other insects may have telomerase. Functional analysis of the TRAS/SART families may provide insight into the mechanism of telomere maintenance in the silkworm. We are grateful to Dr H. Maekawa for providing a plasmid including a Drosophila HSP70 promoter and to Dr T. Shimada for helping to use the luminometer. This work was supported by grants from the Ministry of Education, Science and Culture of Japan (MESCJ) and Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists of MESCJ.
RESULTS
TRAS1 and SART1 were transcribed as full-length units in various tissues and in the BmN4 cell line
Most genomic copies of TRAS1 and SART1 contained the transcription initiation sequences
TRAS1 possesses a promoter in its 5[prime]-UTR
Influence of the first 40 bp conserved in most insect retrotransposons on transcription of TRAS1
The neighboring telomeric repeats inhibit transcription of TRAS1
DISCUSSION
ACKNOWLEDGEMENTS
REFERENCES
This article has been cited by other articles:
This page is run by Oxford University Press, Great Clarendon Street, Oxford OX2 6DP, as part of the OUP Journals
Comments and feedback: jnl.info{at}oup.co.uk
Last modification: 10 Apr 1999
Copyright©Oxford University Press, 1999.
![]()
CiteULike
Connotea
Del.icio.us What's this?
![]()
![]()

![]()
![]()
![]()
P. H. Maxwell, J. M. Belote, and R. W. Levis
Identification of multiple transcription initiation, polyadenylation, and splice sites in the Drosophila melanogaster TART family of telomeric retrotransposons
Nucleic Acids Res.,
November 14, 2006;
34(19):
5498 - 5507.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
T. Matsumoto, M. Hamada, M. Osanai, and H. Fujiwara
Essential Domains for Ribonucleoprotein Complex Formation Required for Retrotransposition of Telomere-Specific Non-Long Terminal Repeat Retrotransposon SART1.
Mol. Cell. Biol.,
July 1, 2006;
26(13):
5168 - 5179.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
K. K. Kojima, T. Matsumoto, and H. Fujiwara
Eukaryotic Translational Coupling in UAAUG Stop-Start Codons for the Bicistronic RNA Translation of the Non-Long Terminal Repeat Retrotransposon SART1
Mol. Cell. Biol.,
September 1, 2005;
25(17):
7675 - 7686.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
N. Maita, T. Anzai, H. Aoyagi, H. Mizuno, and H. Fujiwara
Crystal Structure of the Endonuclease Domain Encoded by the Telomere-specific Long Interspersed Nuclear Element, TRAS1
J. Biol. Chem.,
September 24, 2004;
279(39):
41067 - 41076.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
M. Osanai, H. Takahashi, K. K. Kojima, M. Hamada, and H. Fujiwara
Essential Motifs in the 3' Untranslated Region Required for Retrotransposition and the Precise Start of Reverse Transcription in Non-Long-Terminal-Repeat Retrotransposon SART1
Mol. Cell. Biol.,
September 15, 2004;
24(18):
7902 - 7913.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
Y. Yamamoto, Y. Fujimoto, R. Arai, M. Fujie, S. Usami, and T. Yamada
Retrotransposon-mediated restoration of Chlorella telomeres: accumulation of Zepp retrotransposons at termini of newly formed minichromosomes
Nucleic Acids Res.,
August 1, 2003;
31(15):
4646 - 4653.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
T. Anzai, H. Takahashi, and H. Fujiwara
Sequence-Specific Recognition and Cleavage of Telomeric Repeat (TTAGG)n by Endonuclease of Non-Long Terminal Repeat Retrotransposon TRAS1
Mol. Cell. Biol.,
January 1, 2001;
21(1):
100 - 108.
[Abstract]
[Full Text]
![]()
This Article ![]()
![]()
Abstract
![]()
Print PDF (142K)
![]()
Alert me when this article is cited
![]()
Alert me if a correction is posted
![]()
Services ![]()
![]()
Email this article to a friend
![]()
Similar articles in this journal
![]()
Similar articles in ISI Web of Science
![]()
Similar articles in PubMed
![]()
Alert me to new issues of the journal
![]()
Add to My Personal Archive
![]()
Download to citation manager
![]()
Search for citing articles in:
ISI Web of Science (21)
![]()
Request Permissions ![]()
Commercial Re-use Guidelines
for Open Access NAR Content
![]()
Google Scholar ![]()
![]()
Articles by Takahashi, H.
![]()
Articles by Fujiwara, H.
![]()
Search for Related Content
![]()
PubMed ![]()
![]()
PubMed Citation
![]()
Articles by Takahashi, H.
![]()
Articles by Fujiwara, H.
![]()
Social Bookmarking ![]()
![]()
What's this?