Previous Article | Next Article ![]()
Molecular and Cellular Biology, April 2000, p. 2488-2497, Vol. 20, No. 7
Department of Microbiology and
Immunology,1 Department of Chemistry and
Biochemistry,2 and Molecular Biology
Institute,3 University of California, Los
Angeles, California 90095, and Howard Hughes Medical Institute,
Department of Medicine and Cardiovascular Research Institute,
University of California, San Francisco, California
941434
Received 25 October 1999/Returned for modification 1 December
1999/Accepted 5 January 2000
A number of microaerophilic eukaryotes lack mitochondria but
possess another organelle involved in energy metabolism, the hydrogenosome. Limited phylogenetic analyses of nuclear genes support a
common origin for these two organelles. We have identified a protein of
the mitochondrial carrier family in the hydrogenosome of
Trichomonas vaginalis and have shown that this protein,
Hmp31, is phylogenetically related to the mitochondrial ADP-ATP carrier (AAC). We demonstrate that the hydrogenosomal AAC can be targeted to
the inner membrane of mitochondria isolated from Saccharomyces cerevisiae through the Tim9-Tim10 import pathway used for the assembly of mitochondrial carrier proteins. Conversely, yeast mitochondrial AAC can be targeted into the membranes of hydrogenosomes. The hydrogenosomal AAC contains a cleavable, N-terminal presequence; however, this sequence is not necessary for targeting the protein to
the organelle. These data indicate that the membrane-targeting signal(s) for hydrogenosomal AAC is internal, similar to that found for
mitochondrial carrier proteins. Our findings indicate that the membrane
carriers and membrane protein-targeting machinery of hydrogenosomes and
mitochondria have a common evolutionary origin. Together, they provide
strong evidence that a single endosymbiont evolved into a progenitor
organelle in early eukaryotic cells that ultimately give rise to these
two distinct organelles and support the hydrogen hypothesis for the
origin of the eukaryotic cell.
A variety of phylogenetically
diverse eukaryotes, including ciliates, fungi, and amoeboflagellates,
lack typical eukaryotic organelles such as the mitochondrion and the
peroxisome. Interestingly, these organisms often contain a
double-membrane-bounded organelle called the hydrogenosome (32,
40). Organisms that contain this organelle live in
microaerophilic habitats and rely on the hydrogenosome for fermentative
carbohydrate metabolism. Similar to mitochondria, the hydrogenosome
produces ATP, requiring its exchange with ADP from the cytosol. The
hydrogenosome participates in carbohydrate metabolism, producing
ATP, carbon dioxide, acetate, and molecular hydrogen from the
fermentation of pyruvate. In this respect, hydrogenosomes can be
regarded as the anaerobic equivalents of mitochondria as ATP generators
in these cells.
Although similar in many aspects, mitochondria and hydrogenosomes
differ significantly in structure and function. The absence of cristae,
DNA, F1F0 ATPase, respiratory-chain components,
and cardiolipin and the presence of the enzymes hydrogenase and
pyruvate-ferredoxin oxidoreductase in the hydrogenosome set it apart
from the mitochondrion (17, 32). Based on the detection of
enzymes typically present in anaerobic bacteria, it was proposed that
the hydrogenosome originated from an endosymbiont related to the strict
anaerobe Clostridium (31). Later, emerging
similarities between mitochondria and hydrogenosomes led to the
proposal that hydrogenosomes were converted mitochondria that lost
their respiratory function as a result of movement into anaerobic
habitats (6). Unfortunately, the lack of hydrogenosomal DNA
has precluded a direct analysis of the origin of hydrogenosomes similar
to that carried out using mitochondrial DNA which demonstrated that
mitochondria evolved from an endosymbiont of the Recently, molecular analyses of hydrogenosomal heat shock proteins from
Trichomonas vaginalis (5, 11, 16, 41) have demonstrated a close phylogenetic relationship between the nuclear genes encoding these proteins and their mitochondrial counterparts, suggesting a common symbiotic origin for the two organelles. The latest
theory is one based on the metabolic force behind the symbiotic event A critical step in the evolution of the ancestral endosymbiont to an
organelle would be the evolution of membrane proteins to allow
communication between the organelle and its host cell. Hydrogenosomes
and mitochondria undergo biogenesis by binary fission, followed by
translocation of the nuclear-encoded proteins required for their
functions (37, 40). Hence, membrane proteins would be
essential not only for intracellular transport of substrates and
products but also for translocation of host-encoded proteins during
organelle biogenesis. These membrane proteins would have evolved at the
time of or shortly prior to DNA transfer from symbiont to host nucleus.
One of the first such membrane proteins that would have evolved in the
case of ATP-producing organelles would be an ADP-ATP exchanger that
would provide ATP to the cytosol (6). The evolution of this
translocator would not be complete without developing a translocation
machinery and organelle-targeting signals, all of which would involve a
series of rare mutations. Therefore, the presence of
phylogenetically-related membrane proteins and the use of similar
translocation pathways for mitochondria and hydrogenosomes would reveal
their evolution from a common progenitor if indeed a single ancestor
gave rise to them as proposed (4, 29). Previous studies have
indicated that common import signals are used for targeting proteins to
the matrix of mitochondria and hydrogenosomes (4, 15, 50);
however, studies comparing membrane translocation pathways have not
been reported.
Here, we report the characterization of the first membrane protein
isolated from T. vaginalis hydrogenosomes and demonstrate that it is a member of the mitochondrial carrier family (MCF) (23). Within this group of MCF proteins, phylogenetic
analysis indicates that the hydrogenosomal protein has a common origin with ADP-ATP carrier (AAC) proteins. In vivo and in vitro translocation analyses with hydrogenosomes and mitochondria demonstrate that AACs
from the two organelles utilize similar translocation pathways and rely
on internal signals for specific membrane targeting. These data reveal
the presence of conserved membrane carriers in hydrogenosomes and
mitochondria and indicate the coevolution of membrane protein-targeting
pathways in the two organelles.
Isolation and sequencing of Hmp31.
Hydrogenosomes were
purified from T. vaginalis C1 (ATCC 30001) cultures as
described previously (4). The organelles were alkaline
extracted by incubation for 30 min at 4°C at a concentration of 0.03 mg/ml in 0.1 M sodium carbonate, pH 11.5 (10). The insoluble fraction was separated from the soluble proteins by centrifugation at
200,000 × g for 1 h at 4°C. The insoluble
proteins were resolved by sodium dodecyl sulfate (SDS)-12%
polyacrylamide gel electrophoresis (PAGE) and visualized by Coomassie
brilliant blue staining. A major band of 31 kDa was excised from the
gel and digested with endoproteinase Lys-C (Boehringer Mannheim), and
the resulting peptides were subjected to microsequencing. For
N-terminal sequencing, the protein was transferred to polyvinylidene
difluoride membranes prior to Edman degradation.
Construction and screening of T. vaginalis cDNA and
genomic DNA libraries.
A unidirectional cDNA library was
constructed in the
0270-7306/00/$04.00+0
Copyright © 2000, American Society for Microbiology. All rights reserved.
Presence of a Member of the Mitochondrial Carrier
Family in Hydrogenosomes: Conservation of Membrane-Targeting Pathways
between Hydrogenosomes and Mitochondria
![]()
ABSTRACT
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
![]()
INTRODUCTION
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
-proteobacterial
family (14).
the hydrogen hypothesis (29). This hypothesis
propounds that an
-proteobacterium which produced molecular hydrogen
and carbon dioxide established a symbiosis with a methanogenic archaeon that utilized these as sources of energy. As these gases became depleted with changes in the early Earth's atmosphere, the host (the
archaeon) would have become dependent on its symbiotic partner for its
needs. This association further led to the transfer of genes from the
symbiont to the host and to its successful establishment as an
organelle. This ancient cell is proposed to be ancestral to eukaryotes
that have diverged into a respiratory (mitochondria-containing) or a
fermentative (hydrogenosome-containing) fate, depending on their habitats.
![]()
MATERIALS AND METHODS
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
ZAPII vector (Stratagene) using T. vaginalis poly(A)+ RNA. A fragment of 257 bp, FP257,
was amplified from a cDNA pool of 5 × 107 phage using
the T7 promoter primer on the vector as a reverse primer and LA1 (Table
1), a degenerate forward primer designed from the peptide PIYSGMMQAF. A genomic library constructed in the
FIXII vector (Stratagene) from 9- to 23-kb T. vaginalis
BamHI genomic DNA fragments was screened with FP257, yielding
positive clones with 9.4-kb inserts. A 1.7-kb HindIII
fragment from one of these was further subcloned into pBluescript KS
(Stratagene) and was found by sequencing to contain the complete open
reading frame, 948 bp long. DNA and protein sequences were analyzed
using the MacVector program (Oxford Molecular Group).
TABLE 1.
Oligonucleotide primers used in this study
Constructs for selectable transformation of T. vaginalis.
The
-SCSB-CAT (
-succinyl coenzyme A synthetase
B-chloramphenicol transferase) construct described previously
(8) was modified by introducing an NdeI
restriction site at the initiation codon of the
-SCSB gene by
site-directed mutagenesis. All primers used are listed in Table 1. The
construct was digested with SacI and NdeI to
replace the 1.66-kb 5' untranslated region of
-SCSB with a shorter
SacI-NdeI-digested 340-bp 5' untranslated region fragment generated by PCR using the primers 5'
-SCSBFor and 5'
-SCSBRev. To generate the dihemagglutinin [(HA)2]
epitope tag, primers with KpnI (+2HAKpnBam) and
BamHI (
2HAKpnBam) overhangs, respectively, were hybridized
to each other and used to replace the CAT gene in the modified
-SCSB-CAT template. The Fd- and Hmp31-coding regions were amplified
by PCR to introduce an NdeI site at the initiation codon and
a KpnI after the ultimate codon for restriction digestion,
followed by ligation into the (HA)2 template. For the
Hmp31-(HA)2 and
Hmp31-(HA)2 constructs, the forward PCR primers used were AACHAN and DAACHAN, respectively, with
the reverse primer AACHAK. For the Fd-(HA)2 and
Fd-(HA)2 constructs, forward PCR primers were FdHAN and
DFdHAN, respectively, with the reverse primer FdHAK. For presequence
insertion into
Fd-(HA)2 and
Hmp31-(HA)2,
two sets of primers with NdeI overhangs were designed. For
the Hmp31 presequence, these were AACL+ and AACL
; for the Fd
presequence, these were FdL+ and FdL
. In each case, the primers
corresponded to the presequence and the first three amino acids of the
mature protein. These primer sets were annealed and ligated into the
NdeI site of the
Hmp31-(HA)2 or
Fd-(HA)2 construct to create
FdL
Hmp31-(HA)2 and
Hmp31L
Fd-(HA)2. For the Saccharomyces
cerevisiae AAC1-(HA)2 construct, primers AAC1F and
AAC1R were used. A ClaI restriction fragment containing the
neomycin phosphotransferase (neo) cassette was transferred into each of
the above constructs from the
-TUB-neo construct previously
described (27).
Selectable transformation of T. vaginalis. Electroporation of T. vaginalis strain C1 was carried out as described previously (8) with 30 to 50 µg of circular plasmid DNA. The transformants were selected with G418 (100 µg/ml) prior to crude fractionation and organelle purification.
Crude fractionation of T. vaginalis cells.
T.
vaginalis transformant cultures were harvested and resuspended in
SMDI (250 mM sucrose, 10 mM morpholinepropanesulfonic acid [MOPS]
[pH 7.2], 10 mM dithiothreitol,
N
-p-tosyl-L-lysine chloromethyl ketone
[TLCK] [50 µg/ml], and leupeptin [10 µg/ml]). The cells were
broken in a cell disrupter (Energy Service Co.) and 100 µl of broken
cells was removed as a whole-cell aliquot. One milliliter of the
resuspended broken cells was centrifuged at 12,000 × g
after which an aliquot of 100 µl of supernatant was removed. The
remaining pellet, consisting of the organelles was resuspended in 1 ml
of SMDI, from which an aliquot of 100 µl was removed. Samples were
analyzed by SDS-PAGE.
Phylogenetic analysis. Protein sequences of 25 representative MCF members and Hmp31 were aligned using CLUSTAL W (49) and edited manually with LINEUP (Wisconsin package version 10.0; Genetics Computer Group). Phylogenetic analysis was performed on regions corresponding to residues 24 to 280 of the hydrogenosomal protein using maximum-parsimony analysis from PAUP (48) by 50 rounds of random sequence addition heuristic searches with tree bisection reconnection branch swapping on 100 bootstrap replicates. Only groups that occurred in more than 50% of the bootstrap replicates were included in the consensus tree. Protein distance matrices were calculated using the PROTDIST program of the PHYLIP (9) package on 100 bootstrap replicates of the alignment. A majority consensus tree was generated using the FITCH program (9) with five multiple jumbles per replicate and with global rearrangement. Both trees were rooted using the midpoint method.
Import of radiolabelled proteins into isolated mitochondria. Mitochondria were purified from lactate-grown yeast cells of the tim10-1 mutant and parental strains (13, 21). The protein-coding regions from T. vaginalis Hmp31 and the S. cerevisiae AAC1 genes were subcloned into pSP65 (Promega), and SP6 polymerase was used for in vitro transcription. Proteins were then synthesized in a rabbit reticulocyte lysate in the presence of [35S]methionine. The reticulocyte lysate containing radiolabelled precursor was incubated at 25°C with isolated mitochondria in import buffer (1 mg of bovine serum albumin per ml, 0.6 M sorbitol, 150 mM KCl, 10 mM MgCl2, 2.5 mM EDTA, 2 mM ATP, 2 mM NADH, 20 mM HEPES-KOH [pH 7.4]). Where indicated, the potential across the mitochondrial inner membrane was dissipated with 1 µM valinomycin and 25 µM carbonyl cyanide p-(trifluoromethoxy)phenylhydrazone. Nonimported radiolabelled proteins were removed by treatment with 100 µg of trypsin per ml; trypsin was then inhibited with 200 µg of soybean trypsin inhibitor per ml. Following reisolation, mitochondria were subjected to sodium carbonate extraction, and membrane proteins were separated on SDS-12% PAGE gels, followed by fluorography. Import was quantitated with a laser-scanning densitometer and expressed as a percentage of total import with the four-minute time point of import into wild-type mitochondria set at 100%.
Miscellaneous. Polyclonal antibodies against endogenous Hmp31 were raised in rabbits and used for immunodecoration of Hmp31 in Western blot analysis. Western blot analyses of the epitope-tagged transformants were performed using an anti-hemagglutinin (HA) mouse immunoglobulin G. Detection in both cases was carried out using the respective secondary horseradish peroxidase-linked immunoglobulin G, followed by enhanced chemiluminescence with Amersham substrates. Quantitation of the scanned images was performed using the IQUANT (Molecular Dynamics) program.
Nucleotide sequence accession number. The sequence for Hmp31 has been deposited with GenBank with the accession number AF216971.| |
RESULTS |
|---|
|
|
|---|
Isolation and characterization of Hmp31.
To identify
hydrogenosomal translocases, we isolated integral membrane proteins
from purified T. vaginalis hydrogenosomes by alkaline
extraction (10). One of the most abundant proteins (Fig.
1A), is an approximately 31-kDa protein
designated Hmp31 (for hydrogenosomal membrane protein 31), which was
microsequenced. A fragment of the gene encoding Hmp31 was amplified
from a cDNA library by PCR with primers designed from the peptide
sequences. The complete gene was isolated from a genomic DNA library
and sequenced, revealing a protein-coding region of 316 amino acids of
a calculated molecular mass of 33 kDa. The presence of all six
proteolytic peptide sequences generated by analysis of the purified
hydrogenosomal membrane protein in the translated sequence confirmed
that the correct gene had been isolated and characterized (Fig. 1B).
Interestingly, internal protein sequence comparison of Hmp31 revealed a
tripartite repeat structure with a frequency of about 100 amino acids
(Fig. 1B).
|
2 or
3
relative to the cleavage site, unlike the Hmp31 presequence. Moreover,
the matrix presequences (12 of 13) have either a phenylalanine or an
asparagine at position
1 relative to the cleavage site that is not
found in the Hmp31 presequence. These differences could reflect the
presence of two distinct presequence-processing peptidases in the hydrogenosome.
The Hmp31 N-terminal presequence is not necessary for translocation. To study the role of the Hmp31 presequence in translocation, we have expressed mutant and chimeric hydrogenosomal proteins in T. vaginalis (8, 27) using the influenza HA epitope to detect the proteins. Constructs carrying the HA-tagged genes also carry a neo gene cassette (24) for selection of transformants. The selectable transformants were lysed and separated into soluble cellular and organellar fractions. Hydrogenosomes were then purified from the organellar fraction (4) and further subjected to alkaline sodium carbonate extraction to separate the hydrogenosomal membrane proteins from the soluble proteins. Western analyses using an anti-HA epitope antibody revealed the location of the protein of interest.
To determine whether the Hmp31 presequence is necessary for translocation, we created two selectable transformants (Fig. 2A) that express Hmp31-(HA)2, a full-length Hmp31 construct with the presequence and
Hmp31-(HA)2, a mutant Hmp31 construct lacking amino acid
residues 2 to 12, corresponding to a presequence-minus version of the
protein. Western analysis of a crude fractionation of
Hmp31-(HA)2-transformed cells showed that 48% of the
expressed protein was present as a precursor-sized species in the
soluble cellular fraction and 52% was targeted to the organellar
fraction. The targeted protein corresponded in size to the
presequence-minus protein, indicating that it had been cleaved (Fig.
2B). Within the organelles, Western blot analysis of Hmp31 revealed its
location exclusively in the membrane fraction (Fig. 2B), whereas a
known matrix protein,
-succinyl coenzyme A synthetase
(24) localized to the supernatant fraction (data not shown).
The presence of precursor-sized Hmp31-(HA)2 protein
in the soluble cellular fraction could be due to misfolding or
overabundance. Within the hydrogenosomes from this transformant, the
tagged protein is targeted exclusively to the membrane (Fig. 2B).
Similar analyses on the transformant expressing the presequence-minus
version revealed that removal of the presequence did not affect
translocation to the hydrogenosomes (36% of the expressed protein in
the soluble cellular fraction and 64% in the organellar fraction) and
also resulted in targeting exclusively to the membrane fraction (Fig.
2C). These results show that the Hmp31 presequence is not necessary for
the translocation of Hmp31 to the hydrogenosomal membrane, indicating
that the protein bears internal membrane-targeting signals.
However, we cannot exclude the possibility that the presequence
participates in the translocation process.
|
Hmp31-(HA)2.
Interestingly, this fusion protein was exclusively targeted to the
organellar fraction, with no detectable trace in the soluble cellular
fraction (Fig. 2D). It appears that replacing the Hmp31 presequence
with the Fd presequence results in more efficient targeting. Within the
hydrogenosomes, the fusion protein was found entirely in the membrane
fraction. These data show that the matrix presequence cannot override
the internal membrane-targeting signal in the Hmp31 protein and cause
mistargeting to the soluble fraction of the hydrogenosome. The Fd
presequence was not cleaved upon translocation to the organelle
membrane. This is likely to reflect a requirement for a
matrix-processing peptidase for cleavage of this presequence. The N
terminus of the fusion protein may not be facing the matrix, or the
conformation of the presequence may be altered in the fusion protein
such that it is no longer a substrate for the processing peptidase.
The Hmp31 N-terminal presequence can replace a hydrogenosomal
matrix presequence.
We have previously shown that a
presequence-minus version of the Fd protein,
Fd, was unable to bind
to and hence translocate into isolated hydrogenosomes in vitro
(4). To test whether the Hmp31 presequence can restore
translocation of
Fd, we transformed T. vaginalis cells
with the chimeric construct Hmp31L
Fd-(HA)2, consisting of the Hmp31 presequence fused to presequence-minus Fd (Fig.
3A). As controls, we also generated the
transformants Fd-(HA)2 and
Fd-(HA)2, which
express full-length Fd and presequence-minus Fd, respectively (Fig.
3A). Full-length Fd was targeted to the soluble fraction (91%) of
purified hydrogenosomes as expected from its matrix location (Fig. 3B).
A small amount (9%) was found to be associated with the membranes,
perhaps due to saturation of the organelle with this protein. The
presequence-minus Fd protein was 100% in the soluble cellular fraction
(Fig. 3C), confirming previous in vitro translocation results
(4) and showing that the Fd presequence is necessary for
translocation into the hydrogenosome in vivo. The addition of the Hmp31
presequence to the presequence-minus Fd restored translocation to the
hydrogenosome, with 47% of the protein going to the organellar
fraction and 53% in the soluble cellular fraction. As with the fusion
protein FdL
Hmp31-(HA)2, no cleavage of
the Hmp31 presequence was observed. Within the purified hydrogenosomes,
this fusion protein was found only in the soluble fraction. These data
show that the Hmp31 presequence, although not necessary for
translocation, can substitute for a targeting signal that has been
shown to be necessary to direct proteins to the hydrogenosomal soluble
fraction (4). Furthermore, the results show that the Hmp31
presequence does not appear to bear any membrane-targeting signal.
Thus, it appears that there are at least two types of targeting signals
within the precursor Hmp31 protein.
|
Hmp31 appears to be an inner-membrane hydrogenosomal
protein.
To determine whether Hmp31 is on the outer or the
inner membrane of the hydrogenosome, we subjected intact purified
wild-type T. vaginalis hydrogenosomes to protease
treatment with increasing concentrations of trypsin. The effect on
endogenous Hmp31 was monitored by Western analysis of digested
hydrogenosomes using a polyclonal anti-Hmp31 antibody (Fig.
4). To ascertain whether the Hmp31
protein is inherently resistant to the protease, we also included
Triton X-100-solubilized hydrogenosomes that were exposed to 20 and 1 mg of trypsin per ml, the lowest and the highest trypsin
concentrations used in the assay. The Hmp31 protein appeared to be
unaffected by the trypsin treatment in intact hydrogenosomes, but the
levels were greatly reduced even at concentrations of 20 µg of
trypsin per ml and completely abolished at 1 mg of trypsin per ml in
the solubilized samples (Fig. 4). Control experiments using an antibody
against a matrix protein showed that this protein was also intact in
trypsin-treated samples and was digested upon solubilization of the
organelles with detergent (data not shown). This assay indicates that
the Hmp31 membrane protein is in a protease-protected location,
suggesting that it is an inner-membrane hydrogenosomal protein.
|
Hmp31 is a member of the MCF. Database homology searches with BLAST showed that Hmp31 displays 16 to 30% identity and 34 to 56% similarity with members of the MCF. The carrier family consists of mostly mitochondrial and some peroxisomal membrane proteins of approximately 30 kDa that translocate anions across the organelle membranes. Protein comparisons and phylogenetic analyses have revealed a common ancestral origin for these carriers (23). Hmp31 is structurally similar to these carriers which are also characterized by three repeats, each about 100 amino acids in length (44). The repeats of MCF proteins consist of two transmembrane domains, resulting in proteins that traverse the membrane six times (23). The mitochondrial MCF proteins are inner-membrane proteins, as appears to be the case for Hmp31 (Fig. 4).
A further characteristic is the presence at three equivalent locations of the degenerate signature sequence P (Hy) (D, E) X X (K, R) where Hy designates a hydrophobic residue and X designates any residue (33). This sequence is present in triplicate in all characterized mitochondrial MCF proteins but is found only in the first two repeat units of peroxisomal MCF members. The alignment of approximately 100 amino-acid segments of Hmp31 revealed such a tripartite structure with three highly conserved MCF signature sequences (Fig. 1B). The highest database hits for Hmp31 are to the mitochondrial Graves' disease carriers of unknown function (30% identity and 48% similarity), followed by the peroxisomal-mitochondrial PECA protein (51), a Ca2+-dependent carrier of unknown function (29% identity and 47% similarity), and finally the mitochondrial AACs (21 to 26% identity and 36 to 44% similarity). The latter group of carriers exchange ATP produced in the mitochondrial matrix with ADP from the cytosol. Using maximum-parsimony and distance methods, phylogenetic analyses of Hmp31 with representative members of each subgroup of the MCF proteins showed that this protein groups with high bootstrap values (78% for the maximum-parsimony tree and 100% for the distance tree) with the mitochondrial AACs, the peroxisomal-mitochondrial PECA protein, and the mitochondrial Graves' disease carriers, indicating a common ancestry with these carriers (Fig. 5). The grouping of the hydrogenosomal protein Hmp31 with the PECA protein might suggest a peroxisomal nature for this protein. However, the common ancestral node with the exclusively mitochondrial AAC and Graves' disease carriers, together with the presence of the carrier signature sequence in all three repeats of the hydrogenosomal protein, points more towards a mitochondrial nature (Fig. 1B). Furthermore, the PECA protein is synthesized with an N-terminal extension of 193 amino acids, which is not present in either the T. vaginalis protein or the mitochondrial AACs and Graves' disease carriers.
|
Conservation of charge-pair network residues indicates that Hmp31
is a homologue of the mitochondrial AAC.
The Hmp31 protein shows
26% identity and 44% similarity with the most extensively studied MCF
member, AAC2 of S. cerevisiae. In addition to the conserved
motifs discussed above, Hmp31 contains charge-pair networks (Fig.
6) which have been found to be essential for AAC function in yeast (33, 34). Of the 12 residues in the charge-pair networks, 7 are identical in the hydrogenosomal MCF-like protein, 4 have conserved charge substitutions, and 1 has an
opposite charge. Six of the conserved residues appear to be AAC
specific (33), including the positive charges on the arginine cluster at positions 252 to 254 in the yeast sequence which
are thought to interact with the negative charges on the phosphate
group of ATP on the matrix-facing side of the membrane. The motif
RRRMMM, found in all the AACs characterized to date, is highly
conserved in the hydrogenosomal Hmp31 sequence RKRMML except for a
positive charge substitution and a nonpolar substitution. This level of
conservation drops in the PECA protein and in Graves' disease
carriers, where the sequences at this position are RTRMQA and RRRMQL,
respectively. Similarities in these key motifs provide strong evidence
that the Hmp31 protein is a functional homologue of the mitochondrial
AAC. Moreover, the AAC is the most abundant protein in the
mitochondrial membrane (20), paralleling the high levels of
Hmp31 in hydrogenosomal membranes (Fig. 1A).
|
Import of the yeast AAC into hydrogenosomes in vivo.
Given the
similarity between the hydrogenosomal and mitochondrial AAC membrane
proteins, we examined whether the yeast mitochondrial AAC could be
targeted to hydrogenosomes in vivo. We transformed T. vaginalis cells with a ScAAC1-(HA)2 construct for
expression of C-terminal HA-tagged yeast AAC1 (Fig.
7). The mitochondrial AAC was very
efficiently targeted to the organelle fraction (91% in the organelle
pellet and 9% in the soluble cellular fraction). Within the
hydrogenosomes, 87% of the translocated fusion protein was present in
the membrane fraction and 13% was present in the soluble fraction. We
have obtained similar results using in vitro translocation of
radiolabelled AAC1 into isolated hydrogenosomes, followed by protease
treatment (data not shown). These studies, the first to examine
translocation of proteins to hydrogenosomal membranes, show that a
hydrogenosomal and a mitochondrial membrane protein are targeted with
similar efficiencies. Moreover, they strongly indicate that the
hydrogenosomal translocation machinery can recognize a
membrane-targeting signal on yeast AAC and can correctly direct this
heterologous protein to the membrane.
|
Import of hydrogenosomal AAC into isolated mitochondria in
vitro.
Conversely, we tested whether the hydrogenosomal AAC
protein, Hmp31, can be correctly targeted to yeast mitochondria in
vitro. We found that radiolabelled Hmp31 was translocated into isolated yeast mitochondria in a membrane potential-dependent manner, as was the
mitochondrial AAC (39). As described above for translocation of hydrogenosomal and yeast AAC into hydrogenosomes, both proteins were
specifically inserted into the mitochondrial membrane (Fig. 8A). Depletion of the membrane potential
resulted in a reduction in the amount of protease-protected,
membrane-inserted T. vaginalis Hmp31, in accordance with the
pathway taken by the mitochondrial AAC (39).
|
| |
DISCUSSION |
|---|
|
|
|---|
We have identified and characterized the first membrane protein from hydrogenosomes and found that it is phylogenetically related to the mitochondrial carrier family of proteins. Specifically, this protein is a homologue of the mitochondrial ADP-ATP carrier, an inner-membrane component that exchanges ADP and ATP between the cytosol and the organelle. Like the mitochondrial AAC, Hmp31 is one of the most abundant membrane proteins of hydrogenosomes. We have established a system to study targeting of proteins to hydrogenosomes in vivo and have shown that the hydrogenosomal AAC contains an N-terminal, cleavable presequence; however, this sequence is not necessary for translocation of the protein into hydrogenosomes. These data indicate that internal sequences are used to target hydrogenosomal membrane carriers, similar to that known for targeting mitochondrial carriers (38, 46). Our results also demonstrate that the translocation machinery of hydrogenosomes can effectively recognize and target a mitochondrial AAC. Conversely, we have shown that the hydrogenosomal AAC homologue can be targeted to mitochondrial membranes in vitro utilizing the same pathway used for assembly of mitochondrial inner-membrane carrier proteins (19, 21, 26, 39, 43, 45). This compatibility in the translocation machinery in the two organelles is quite remarkable given the biochemical differences between hydrogenosomes and mitochondria.
We have shown that the mature Hmp31 protein, lacking the presequence, is translocation competent. Unlike Hmp31, however, the majority of MCF proteins do not possess an N-terminal cleavable targeting sequence. Two notable exceptions are the plant AACs and mammalian phosphate carriers (35). In plants, the AACs are synthesized with long N-terminal presequences (12) which have been found not to be necessary for translocation (30, 52). Although the Hmp31 presequence is shorter than those on plant AACs, it has an amino acid composition similar to that of the N-terminal extremity of plant AAC presequences. The latter generally start with a Met-Ala dipeptide, are enriched in the amino acids glutamine and alanine, and likewise have a negatively charged amino acid residue (12).
The other MCF proteins that bear a cleavable N-terminal presequence, the mammalian phosphate carriers, are synthesized with an extension of 44 to 49 amino acids. In this case, deletion of the presequence still results in an import competent protein; however, efficiency is reduced to 40 to 50% of normal. Thus, although the presequence is not necessary for translocation, it probably has a role in enhancing import efficiency. It was further observed that the presequence alone was able to target a passenger protein into mitochondria, albeit at a very low efficiency. The mitochondrial phosphate carriers therefore seem to have multiple targeting signals (53). These comparisons reveal a similarity between import signals on membrane proteins of mitochondria and hydrogenosomes. Both appear to rely upon internal membrane-targeting signals; however, when present, a cleavable N-terminal presequence also appears to be capable of restoring translocation of proteins to the organelles.
We have demonstrated that the addition of a matrix-targeting signal to the N terminus of Hmp31 that lacks its own presequence does not cause mistargeting to the matrix, indicating that the matrix-targeting signal is not capable of overriding the internal membrane signal. A similar experiment was done in mitochondria using the presequence of preornithine carbamyltransferase, a matrix protein, fused to the mitochondrial uncoupling protein, an inner-membrane protein of the MCF which normally lacks a presequence. In this case, the chimeric protein was translocated across the inner membrane into the matrix, where the matrix presequence was cleaved and the protein was not inserted into the membrane (28). These results indicate that the membrane-targeting signal on the uncoupling protein can be overridden by the matrix-targeting signal, contrary to our observations using a chimeric construct consisting of a matrix presequence fused to Hmp31.
Translocation of mitochondrial inner-membrane proteins, such as the AAC and the Tim components, follows a different import pathway than that of proteins destined for the mitochondrial matrix (19, 21, 26, 39, 43, 45). After passing through the Tom complex, these precursor proteins are assisted by a family of small Tim proteins, Tim8, Tim9, Tim10, and Tim13, which guide the precursors through the intermembrane space to inner-membrane translocase complexes that mediate membrane insertion. We have shown that hydrogenosomal AAC is imported efficiently into yeast mitochondria and that its import requires the Tim9-Tim10 complex as does mitochondrial AAC. These data indicate that the import pathway used for inner-membrane proteins is conserved between hydrogenosomes and mitochondria. We have also demonstrated that yeast mitochondrial AAC is translocated efficiently into hydrogenosomes, providing additional evidence for a conservation of the import pathway for inner-membrane proteins. One would therefore predict that homologues of the Tim proteins are present in hydrogenosomes. Our data also indicate that a common inner-membrane-targeting signal for MCF proteins is used by the two organelles. Specific targeting sequences have not been identified for mitochondrial inner-membrane proteins; however, deletion studies of yeast AAC indicate that the targeting information is internal and occurs in two distinct parts of the protein (38, 46).
Previous studies have shown that the matrix-targeting pathway is conserved between hydrogenosomes and mitochondria (4, 15, 50). Taken together with our current studies on the translocation of proteins into the membranes of hydrogenosomes, it can be concluded that both matrix and membrane-targeting pathways are far more similar than is reasonable to predict if these two organelles were to have had independent origins.
Current phylogenetic analyses support the theory that the mitochondrion
originated from an endosymbiotic ancestor related to the contemporary
members of the
-proteobacteria (14). The debate over the
origins of the hydrogenosome and of the mitochondrion basically
involves two possibilities: either the two organelles evolved from two
independent symbioses of related
-proteobacteria, or there was one
common bacterial ancestor that subsequently diverged into the
specific organelles.
Phylogenetic analyses of heat shock proteins in hydrogenosome-containing T. vaginalis (5, 11, 16, 41) as well as in Giardia lamblia (42) and Entamoeba histolytica (7), two eukaryotes that harbor neither organelle but appear to have lost mitochondria secondarily, indicate that these proteins are closely related to mitochondrial heat shock proteins. These findings, however, are still compatible with independent symbioses of two related proteobacteria (47) bearing similar heat shock proteins, which are indeed still conserved in their modern-day relatives. On the other hand, the MCF proteins have no known proteobacterial or archaebacterial homologues, and statistical analyses of the MCF repeat units have indicated that the tripartite proteins have evolved from the triplication of an ancestral gene around the time mitochondria first appeared in eukaryotes (23) before diverging into the various specific solute carriers. Hence, the origin of these MCF proteins appears to be closely associated with the symbiotic event that gave rise to the mitochondrion. It seems unlikely that such a triplication event occurred independently in different progenitors of mitochondria and hydrogenosomes, to give rise to the mitochondrial carrier proteins in the former case and to Hmp31 in the latter case. Moreover, mitochondrial AAC proteins are phylogenetically unrelated to the ADP-ATP translocators of Rickettsia, the microbe which is most closely related to mitochondria (3). This implies that the event which gave rise to mitochondrial AAC is different from that in which the mitochondria as a whole originated (2).
Therefore, the results presented here, showing the presence of a mitochondrial-type membrane carrier protein in hydrogenosomes and a compatible translocation mechanism for targeting membrane proteins into both organelles, provide unequivocal evidence that the two organelles originated from a common progenitor organelle. This does not, however, preclude the possibility that later symbioses delivered specific proteins to the two organelles after their divergence.
| |
ACKNOWLEDGMENTS |
|---|
We thank Miriam Makabi for technical assistance, Harry Hahn for help with the phylogeny programs, and Gottfried Schatz, James Lake, Alexander van der Bliek, and members of our laboratory for helpful advice and critical comments on the manuscript.
This work was supported by National Institutes of Health grant AI27857 to P.J.J., postdoctoral fellowships from the Damon Runyon-Walter Winchell Cancer Research Foundation and the National Science Foundation to C.M.K. and from the Deutsche Forschungsgemeinschaft to E.P., and a predoctoral training grant (NIH AI07323) to P.J.B. P.J.J. is a recipient of a Burroughs-Wellcome Scholar in Molecular Parasitology award.
| |
FOOTNOTES |
|---|
* Corresponding author. Mailing address: UCLA Department of Microbiology and Immunology, 1602 Molecular Sciences Building, 405 Hilgard Ave., Los Angeles, CA 90095-1489. Phone: (310) 825-4870. Fax: (310) 206-5231. E-mail: johnsonp{at}ucla.edu.
| |
REFERENCES |
|---|
|
|
|---|
| 1. | Adam, A., M. Endres, C. Sirrenberg, F. Lottspeich, W. Neupert, and M. Brunner. 1999. Tim9, a new component of the TIM22.54 translocase in mitochondria. EMBO J. 18:313-319[CrossRef][Medline]. |
| 2. | Andersson, S. G. E., and C. G. Kurland. 1999. Origins of mitochondria and hydrogenosomes. Curr. Opin. Microbiol. 2:535-541[CrossRef][Medline]. |
| 3. | Andersson, S. G. E., A. Zomorodipour, J. P. Andersson, T. Sicheritz-Ponten, U. C. M. Alsmark, R. M. Podowski, A. K. Naslund, A.-S. Eriksson, H. H. Winkler, and C. G. Kurland. 1998. The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature 396:133-140[CrossRef][Medline]. |
| 4. | Bradley, P. J., C. J. Lahti, E. Plümper, and P. J. Johnson. 1997. Targeting and translocation of proteins into the hydrogenosome of the protist Trichomonas: similarities with mitochondrial protein import. EMBO J. 16:3484-3493[CrossRef][Medline]. |
| 5. |
Bui, E. T. N.,
P. J. Bradley, and P. J. Johnson.
1996.
A common evolutionary origin for mitochondria and hydrogenosomes.
Proc. Natl. Acad. Sci. USA
93:9651-9656 |
| 6. | Cavalier-Smith, T. 1987. The simultaneous symbiotic origin of mitochondria, chloroplasts, and microbodies. Ann. N. Y. Acad. Sci. 503:55-71[Medline]. |
| 7. |
Clark, C. G., and A. J. Roger.
1995.
Direct evidence for secondary loss of mitochondria in Entamoeba histolytica.
Proc. Natl. Acad. Sci. USA
92:6518-6521 |
| 8. |
Delgadillo, M. G.,
D. R. Liston,
K. Niazi, and P. J. Johnson.
1997.
Transient and selectable transformation of the parasitic protist Trichomonas vaginalis.
Proc. Natl. Acad. Sci. USA
94:4716-4720 |
| 9. |
Felsenstein, J.
1993.
PHYLIP phylogeny inference package, version 3.57c.
University of Washington, Seattle.
|
| 10. |
Fujiki, Y.,
A. L. Hubbard,
S. Fowler, and P. B. Lazarow.
1982.
Isolation of intracellular membranes by means of sodium carbonate treatment: application to endoplasmic reticulum.
J. Cell Biol.
93:97-102 |
| 11. |
Germot, A.,
H. Philippe, and H. Le Guyader.
1996.
Presence of a mitochondrial-type 70-kDa heat shock protein in Trichomonas vaginalis suggests a very early mitochondrial endosymbiosis in eukaryotes.
Proc. Natl. Acad. Sci. USA
93:14614-14617 |
| 12. | Glaser, E., S. Sjöling, M. Tanudji, and J. Whelan. 1998. Mitochondrial protein import in plants. Signals, sorting, targeting, processing and regulation. Plant Mol. Biol. 38:311-338[CrossRef][Medline]. |
| 13. | Glick, B. S., and L. A. Pon. 1995. Isolation of highly purified mitochondria from Saccharomyces cerevisiae. Methods Enzymol. 260:213-223[Medline]. |
| 14. |
Gray, M. W.,
G. Burger, and B. F. Lang.
1999.
Mitochondrial evolution.
Science
283:1476-1481 |
| 15. | Hausler, T., Y. D. Stierhof, J. Blattner, and C. Clayton. 1997. Conservation of mitochondrial targeting sequence function in mitochondrial and hydrogenosomal proteins from the early-branching eukaryotes Crithidia, Trypanosoma and Trichomonas. Eur. J. Cell Biol. 73:240-251[Medline]. |
| 16. | Horner, D. S., R. P. Hirt, S. Kilvington, D. Lloyd, and T. M. Embley. 1996. Molecular data suggest an early acquisition of the mitochondrion endosymbiont. Proc. R. Soc. Lond. Series B. Biol. Sci. 263:1053-1059[Medline]. |
| 17. | Johnson, P. J., P. J. Bradley, and C. J. Lahti. 1995. Cell biology of trichomonads: protein targeting to the hydrogenosome, p. 399-411. In J. C. Boothroyd, and R. Komuniecki (ed.), Molecular approaches to parasitology. Wiley-Liss, Inc., New York, N.Y. |
| 18. |
Johnson, P. J.,
C. E. d'Oliveira,
T. E. Gorrell, and M. Müller.
1990.
Molecular analysis of the hydrogenosomal ferredoxin of the anaerobic protist Trichomonas vaginalis.
Proc. Natl. Acad. Sci. USA
87:6097-6101 |
| 19. |
Kerscher, O.,
J. Holder,
M. Srinivasan,
R. S. Leung, and R. E. Jensen.
1997.
The Tim54p-Tim22p complex mediates insertion of proteins into the mitochondrial inner membrane.
J. Cell Biol.
139:1663-1675 |
| 20. | Klingenberg, M. 1976. The enzymes of biological membranes, vol. 3. Membrane transport. , p. 383-438. Plenum Publishing Corp., New York, N.Y. |
| 21. |
Koehler, C. M.,
E. Jarosch,
K. Tokatlidis,
K. Schmid,
R. J. Schweyen, and G. Schatz.
1998.
Import of mitochondrial carriers mediated by essential proteins of the intermembrane space.
Science
279:369-373 |
| 22. | Koehler, C. M., S. Merchant, W. Oppliger, K. Schmid, E. Jarosch, L. Dolfini, T. Junne, G. Schatz, and K. Tokatlidis. 1998. Tim9p, an essential partner subunit of Tim10p for the import of mitochondrial carrier proteins. EMBO J. 17:6477-6486[CrossRef][Medline]. |
| 23. | Kuan, J., and M. H. Saier, Jr. 1993. The mitochondrial carrier family of transport proteins: structural, functional, and evolutionary relationships. Crit. Rev. Biochem. Mol. Biol. 28:209-233[Medline]. |
| 24. |
Lahti, C. J.,
C. E. d'Oliveira, and P. J. Johnson.
1992.
-succinyl CoA synthetase from Trichomonas vaginalis is a soluble hydrogenosomal protein with an amino-terminal sequence that resembles mitochondrial presequences.
J. Bacteriol.
174:6822-6830 |
| 25. | Lahti, C. J., and P. J. Johnson. 1991. Trichomonas vaginalis hydrogenosomal proteins are synthesized on free polyribosomes and may undergo processing upon maturation. Mol. Biochem. Parasitol. 46:307-310[CrossRef][Medline]. |
| 26. | Leuenberger, D., N. A. Bally, G. Schatz, and C. M. Koehler. 1999. Different import pathways through the mitochondrial intermembrane space for inner membrane proteins. EMBO J. 18:4816-4822[CrossRef][Medline]. |
| 27. |
Liston, D. R., and P. J. Johnson.
1999.
Analysis of a ubiquitous promoter element in a primitive eukaryote: early evolution of the initiator element.
Mol. Cell. Biol.
19:2380-2388 |
| 28. |
Liu, X. Q.,
K. B. Freeman, and G. C. Shore.
1990.
An amino-terminal signal sequence abrogates the intrinsic membrane-targeting information of mitochondrial uncoupling protein.
J. Biol. Chem.
265:9-12 |
| 29. | Martin, W., and M. Müller. 1998. The hydrogen hypothesis for the first eukaryote. Nature 392:37-41[CrossRef][Medline]. |
| 30. | Mozo, T., K. Fischer, U. I. Flugge, and U. K. Schmitz. 1995. The N-terminal extension of the ADP/ATP translocator is not involved in targeting to plant mitochondria in vivo. Plant J. 7:1015-1020[CrossRef][Medline]. |
| 31. | Müller, M. 1980. The hydrogenosome. In G. W. Gooday, D. Lloyd, and A. P. J. Trinci (ed.), The eukaryotic microbial cell. Cambridge University Press, Cambridge, England. |
| 32. | Müller, M. 1993. The hydrogenosome. J. Gen. Microbiol. 139:2879-2889[Medline]. |
| 33. | Nelson, D. R., C. M. Felix, and J. M. Swanson. 1998. Highly conserved charge-pair networks in the mitochondrial carrier family. J. Mol. Biol. 277:285-308[CrossRef][Medline]. |
| 34. | Nelson, D. R., J. E. Lawson, M. Klingenberg, and M. G. Douglas. 1993. Site-directed mutagenesis of the yeast mitochondrial ADP/ATP translocator. Six arginines and one lysine are essential. J. Mol. Biol. 230:1159-1170[CrossRef][Medline]. |
| 35. | Neupert, W. 1997. Protein import into mitochondria. Annu. Rev. Biochem. 66:863-917[CrossRef][Medline]. |
| 36. | Pfanner, N. 1998. Mitochondrial import: crossing the aqueous intermembrane space. Curr. Biol. 8:R262-R265[CrossRef][Medline]. |
| 37. | Pfanner, N., E. A. Craig, and A. Hönlinger. 1997. Mitochondrial preprotein translocase. Annu. Rev. Cell Dev. Biol. 13:25-51[CrossRef][Medline]. |
| 38. |
Pfanner, N.,
P. Hoeben,
M. Tropschug, and W. Neupert.
1987.
The carboxyl-terminal two-thirds of the ADP/ATP carrier polypeptide contains sufficient information to direct translocation into mitochondria.
J. Biol. Chem.
262:14851-14854 |
| 39. |
Pfanner, N., and W. Neupert.
1987.
Distinct steps in the import of ADP/ATP carrier into mitochondria.
J. Biol. Chem.
262:7528-7536 |
| 40. | Plümper, E., P. J. Bradley, and P. J. Johnson. 1998. Implications of protein import on the origin of hydrogenosomes. Protist 149:303-311. |
| 41. |
Roger, A. J.,
C. G. Clark, and W. F. Doolittle.
1996.
A possible mitochondrial gene in the early-branching amitochondriate protist Trichomonas vaginalis.
Proc. Natl. Acad. Sci. USA
93:14618-14622 |
| 42. |
Roger, A. J.,
S. G. Svärd,
J. Tovar,
C. G. Clark,
M. W. Smith,
F. D. Gillin, and M. L. Sogin.
1998.
A mitochondrial-like chaperonin 60 gene in Giardia lamblia: evidence that diplomonads once harbored an endosymbiont related to the progenitor of mitochondria.
Proc. Natl. Acad. Sci. USA
95:229-234 |
| 43. | Ryan, M. T., and N. Pfanner. 1998. The preprotein translocase of the mitochondrial outer membrane. Biol. Chem. 379:289-294[Medline]. |
| 44. | Saraste, M., and J. E. Walker. 1982. Internal sequence repeats and the path of polypeptide in mitochondrial ADP/ATP translocase. FEBS Lett. 144:250-254[CrossRef][Medline]. |
| 45. | Sirrenberg, C., M. Endres, H. Fölsch, R. A. Stuart, W. Neupert, and M. Brunner. 1998. Carrier protein import into mitochondria mediated by the intermembrane proteins Tim10/Mrs11 and Tim12/Mrs5. Nature 391:912-915[CrossRef][Medline]. |
| 46. | Smagula, C. S., and M. G. Douglas. 1988. ADP-ATP carrier of Saccharomyces cerevisiae contains a mitochondrial import signal between amino acids 72 and 111. J. Cell. Biochem. 36:323-327[CrossRef][Medline]. |
| 47. | Sogin, M. 1997. History assignment: when was the mitochondrion founded? Curr. Opin. Genet. Dev. 7:792-799[CrossRef][Medline]. |
| 48. | Swofford, D. L. 1997. PAUP: phylogenetic analysis using parsimony, version 4.0. Sinauer Associates, Inc., Sunderland, Mass. |
| 49. |
Thompson, J. D.,
D. G. Higgins, and T. J. Gibson.
1994.
CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.
Nucleic Acids Res.
22:4673-4680 |
| 50. | van der Giezen, M., J. A. Kiel, K. A. Sjollema, and R. A. Prins. 1998. The hydrogenosomal malic enzyme from the anaerobic fungus Neocallimastix frontalis is targeted to mitochondria of the methylotrophic yeast Hansenula polymorpha. Curr. Genet. 33:131-135[CrossRef][Medline]. |
| 51. |
Weber, F. E.,
G. Minestrini,
J. H. Dyer,
M. Werder,
D. Boffelli,
S. Compassi,
E. Wehrli,
R. M. Thomas,
G. Schulthess, and H. Hauser.
1997.
Molecular cloning of a peroxisomal Ca2+-dependent member of the mitochondrial carrier superfamily.
Proc. Natl. Acad. Sci. USA
94:8509-8514 |
| 52. | Winning, B. M., C. J. Sarah, P. E. Purdue, C. D. Day, and C. J. Leaver. 1992. The adenine nucleotide translocator of higher plants is synthesized as a large precursor that is processed upon import into mitochondria. Plant J. 2:763-773[CrossRef][Medline]. |
| 53. |
Zara, V.,
F. Palmieri,
K. Mahlke, and N. Pfanner.
1992.
The cleavable presequence is not essential for import and assembly of the phosphate carrier of mammalian mitochondria but enhances the specificity and efficiency of import.
J. Biol. Chem.
267:12077-12081 |
This article has been cited by other articles: