|
|
|
|
Vol. 14, Issue 6, 2314-2326, June 2003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




* Department of Molecular and Cellular Biology, University of California
Berkeley, Berkeley, California 94720;
Department of Microbiology and Immunology, University of California, San
Francisco, San Francisco, California 94143-0414
Submitted January 20, 2002;
Revised February 13, 2003;
Accepted February 19, 2003
Monitoring Editor: Tim Stearns
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
H. capsulatum exists in two morphological forms: a mycelial (or filamentous) form in soil and a yeast form in the host. The mycelial form produces vegetative spores, or conidia.
Conidia or mycelial fragments are inhaled by the host and then taken up by
macrophages and other phagocytic cells
(Eissenberg and Goldman,
1994
). Once inside the host, conversion of the mycelial form to
the budding yeast form is triggered within hours. Yeast cells evade killing
and multiply within macrophages (Bullock,
1993
). Subsequently, yeast cells spread to multiple organs of the
reticuloendothelial system such as the spleen, liver, lymph nodes, and bone
marrow. In patients with disseminated disease, a variety of additional organs
can be colonized (Eissenberg and Goldman,
1991
).
The ability of H. capsulatum to grow in the mycelial form in soil
and shift to the yeast form in the host is important for infection and
disease. On disturbance of the soil, mycelial fragments and associated conidia
are aerosolized and inhaled by the host, unlike yeast cells. Once introduced
into the host, transformation to the yeast form must occur for the fungus to
survive and proliferate (Maresca et
al., 1977
; Medoff et
al., 1986
). Despite the fundamental roles that these two
forms play in infection and disease, little is known about their molecular
differences or what regulates the transition between them.
The transformation between mycelial and yeast phases, or vice versa, can be
recapitulated in culture by shifting the growth temperature from 25 to
37°C, or the reverse (Maresca and
Kobayashi, 1989
; Maresca
et al., 1994
). This characteristic has made it possible
to identify phase-specific genes. Five genes have been identified that are
expressed predominantly in the yeast phase
(Keath et al., 1989
;
Di Lallo et al.,
1994
; Keath and Abidi,
1994
; Gargano et al.,
1995
; Patel et al.,
1998
). One of these genes, CBP1 (calcium binding
protein), is a virulence factor in the host
(Sebghati et al.,
2000
). Several genes specific to the mycelial phase have also been
identified (Harris et al.,
1989a
,b
;
Tian and Shearer, 2001
;
Johnson et al., 2002
;
Tian and Shearer, 2002
).
To develop a more complete description of the two morphological phases, we
performed a large-scale analysis of gene expression in H. capsulatum.
Because the sequence of the H. capsulatum genome has not yet been
completed, we constructed a 10,000-element array containing random genomic
fragments. Using this array, we identified
500 clones whose expression
was differentially induced in either the yeast or mycelial forms, including
several potential regulatory genes. This work sets the stage for uncovering
the function of these genes in the growth phases of H. capsulatum as
well as for applying genomic approaches to other questions in this fungal
pathogen.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Library
Genomic DNA was isolated from 100 ml of the virulent G217B strain by using
genomic tips and genomic DNA buffer set (both from QIAGEN, Valencia, CA). For
construction of the mini-array, DNA was partially digested with
Sau3AI and size fractionated (0.52 kb) on a 1% low melt
agarose gel. This procedure was repeated, size fractionating from 1 to 2 kb,
to create a library for the large array. The resulting fragments were ligated
into pBluescript KS+ precut with BamHI. The library was
transformed into SUREII cells (Stratagene, La Jolla, CA).
Colony Polymerase Chain Reaction (PCR) Amplification and Microarray
Construction
A 96-well format was used to grow 1536 (mini-array) or 9600 (large array)
individual bacterial cultures, each harboring an independent library
transformant (see library construction above). Individual colonies from the
library transformation were inoculated into 100 µl of Luria Broth plus 100
µg/ml carbenicillin in individual wells of 96-well plates and incubated for
1618 h at 37°C on an orbital shaker at 150 rpm. Two microliters
from each well were used for PCR amplification of the inserts by using
M13-forward (5'-GTTTTCCCAGTCACGAC-3') and M13-reverse
(5'-GCGGATAACAATTTCACACAGG-3') primers; these primers were
complementary to the vector. Glycerol was added to the remaining bacterial
cultures to 25%, and the 96-well dishes frozen at -80°C. The PCR
products were analyzed on 1% agarose gels and then precipitated, washed, and
printed on glass slides as described previously
(DeRisi et al.,
1997
). Later iterations of the array included
400 clones from
a yeast-phase cDNA library.
Sequencing
Copies of the 96-well bacterial cultures were sent to Incyte Genomics (Palo
Alto, CA) and the Genome Sequencing Center (Washington University, St. Louis,
MO). Incyte Genomics (Palo Alto, CA) sequenced one side of each clone by using
the M13-forward primer. As part of the ongoing genome project, the Genome
Sequencing Center sequenced each clone by using the M13-forward and
M13-reverse primers. The Incyte sequence is available at
http://gregor.berkeley.edu,
and the Genome Sequencing Center sequence is available at
http://www.genome.wustl.edu/projects/hcapsulatum/.
RNA Preparation
Cultures of yeast and mycelial phase H. capsulatum were harvested
by filtration. Cells were disrupted in RNA extraction buffer (4 M guanidine
thiocyanate, 25 mM sodium citrate, 0.5% sarkosyl
[N-lauroyl-sarcosine], 0.1 M
-mercaptoethanol) by vortexing in
the presence of glass beads. RNA was extracted once with acid phenol,
chloroform, and 0.1 M NaOAc, pH 4.0, and then extracted twice with
equilibrated (pH 8) phenol/chloroform. RNA was then precipitated with
isopropanol, washed with 70% ethanol and resuspended in double distilled
H2O.
cDNA Synthesis, Labeling, and Analysis
Fluorescently labeled cDNA was made by incorporating amino-allyl dUTP
during reverse transcription of poly-adenylated (poly-A)selected RNA.
Cy3 or Cy5 dyes (Amersham Biosciences, Piscataway, NJ) were coupled to the
amino-allyl group as described previously
(DeRisi et al.,
1997
). For the yeast stationary-phase experiments, an equal mass
of each time point was pooled to generate a reference sample, which was
labeled with Cy3. Each time point was individually labeled with Cy5 and
competitively hybridized against the reference sample. Yeast and mycelial
cDNAs were labeled with Cy3 and Cy5, respectively. Dyes were reversed for the
reverse fluor control.
Northern Analysis
Total RNA (510 µg) was separated on a 1.5% denaturing
agaroseformaldehyde gel and transferred to a GeneScreen Plus membrane
(PerkinElmer Life Sciences, Boston, MA). To generate probes, either the entire
insert of a library clone was isolated through restriction digest, or
gene-specific primers were used to amplify regions specific to the gene of
interest. Oligonucleotide sequences can be found in supplemental materials.
Probes were labeled using rediprime (Amersham Biosciences) and
[
-32P]dCTP. The membrane was probed in hybridization buffer
(1 M NaCl, 50% formamide, 1% SDS, 10% dextran sulfate, 33 µg/ml salmon
sperm DNA) at 42°C overnight, and then washed twice in 2x SSC, 1%
SDS at 65°C for 1 h each before exposure to film and PhosphorImager screen
(Molecular Dynamics, Amersham Biosciences, Piscataway, NJ).
Data Analysis
Arrays were scanned on a GenePix 4000B scanner (Axon Instruments, Foster
City, CA) and analyzed using GENEPIX PRO version 3.0, NOMAD
(http://derisilab5.ucsf.edu/NOMAD),
CLUSTER, and TREEVIEW (Eisen et
al., 1998
). For yeast stationary-phase experiments, because
the reference was a pooled sample, ratio measurements from the time-course
data were normalized relative to the first time point (day 2). The expression
ratios for each clone on a given array were divided by the corresponding
ratios measured from the day 2 array. CLUSTER analysis was performed on two
independent stationary-phase experiments and three yeast-versus-mycelial
experiments. Only one stationary-phase experiment is shown in
Figure 3. Data from all
experiments are available in supplemental materials.
|
BLAST Analysis and E-Values for Homologs
For initial annotation of partial array clone sequences, DNA sequences from
Incyte Genomics were compared against available protein databases by using
BLASTx (Gish and States,
1993
). BLASTx hits with an E-value less than or equal to 1 x
10-6 were evaluated. For further annotation of yeast and mycelial
induced clones, array clone end-sequences from Incyte Genomics and the Genome
Sequencing Center were compared against the Genome Sequencing Center H.
capsulatum genome sequencing project contigs by using BLASTn. Contigs
corresponding to microarray clones were compared against the National Center
for Biotechnology Information nr protein database by using BLASTx. Array
clones and potential BLASTx hits were mapped to contig sequences and evaluated
for overlap. Array clones that clearly contained a single open reading frame
(ORF) BLASTx hit with an E-value ≤1 x 10-12 were
annotated. Those that contained more than one ORF were not annotated.
5' and 3' Rapid Amplification of cDNA Ends (RACE) and
Reverse Transcription-Polymerase Chain Reaction (RT-PCR)
Poly-A RNA was purified from total RNA isolated from yeast or mycelia by
using Oligotex mRNA kit (QIAGEN). FirstChoice RNA-ligase mediated (RLM)-RACE
kit (Ambion, Austin, TX) was used to identify 5' and 3' ends of
cDNAs as specified in the kit. Oligonucleotide sequences used for these
analyses can be found in supplemental materials. The coding sequence for open
reading frames was amplified by PCR with gene-specific primers from cDNA
synthesized from poly-A RNA. PCR products were cloned using TOPO-TA
(Stratagene) and sequenced using M13-forward, M13-reverse, and gene-specific
primers as needed.
| RESULTS |
|---|
|
|
|---|
25 Mb)
(Carr and Shearer, 1998
30 H.
capsulatum genes with available genomic and cDNA sequences predicted that
the average intron size is fairly small (100 nucleotides) with 06
introns per gene. Thus, intron sequences were unlikely to interfere with the
ability of a cDNA probe to bind its cognate genomic fragment. These
expectations were tested by construction of a 1500-element mini-array
containing random genomic inserts of
1 kb (see MATERIALS AND METHODS). To
determine the fraction of elements in this array that contained coding
sequence, we carried out a competitive hybridization with two samples: 1)
genomic DNA labeled with the green fluor Cy3 and 2) cDNA made from yeast-phase
cells labeled with the red fluor Cy5. As expected, the genomic DNA gave a
signal for all of the array clones. The cDNA hybridized to
75% of the
array clones (our unpublished data). Because the cDNA represented genes
expressed only under one growth condition, these data indicate that a minimum
of 75% of the array elements contained coding sequence. We therefore proceeded
to construct a large-scale genomic shotgun array (see MATERIALS AND METHODS).
Previously identified H. capsulatum genes were also spotted at known
locations on the array.
|
Array Content
Based on the estimated genome size, the number of clones on the array, and
the average size of the clones, the array covers approximately one-third of
the genome. For purposes of clone identification, we obtained partial sequence
information (one or two sequence reads per clone) for the majority of clones
on the array (see MATERIALS AND METHODS). BLASTx analysis (using an E-value
cut-off of 10-6) of single sequence reads revealed that a diverse
set of gene families was represented on the array
(Figure 1B). Although only 10%
of these reads identified a clear homolog by BLAST, mapping and annotation of
the remainder of the clones is ongoing as the genome sequencing project
develops. Only 2% of the clones contained ribosomal DNA (rDNA); thus, a large
percentage of clones likely contain protein-coding sequence.
Identification of Yeast and Mycelial Phase-regulated Clones
To compare the gene expression profiles of yeast cells with the mycelial
form of H. capsulatum, yeast cells were grown to mid-logarithmic
phase at 37°C. Mycelial cells were grown by shifting a dilution of yeast
cells from 37°C to room temperature; mycelial form cells were allowed to
accumulate for
34 wk without shaking. Microscopic observation of
these cells confirmed the expected morphology and revealed that the mycelial
sample also contained conidia, or vegetative spores, which arise from the
mycelial form. Polyadenylated RNA was isolated from both yeast and mycelial
cells. Lower yields of poly-adenylated RNA (but not total RNA) were obtained
from mycelia than from yeast for unknown reasons. cDNAs generated from yeast
and mycelial RNA were labeled differentially (Cy3 for the yeast sample and Cy5
for the mycelial sample) and subjected to competitive hybridization on the
microarray. As expected from the pilot mini-array, the majority of the array
elements contained protein coding sequence, as indicated by the ability of the
cDNA probe to bind to most array elements
(Figure 2A). A histogram
showing the distribution of the signal for each fluorophore over the entire
array indicated that at least 500 genes were expressed at significantly higher
levels (≥5-fold) in one phase compared with the other
(Figure 2B).
|
CBP1 and yps-3, two previously identified yeast-specific
genes (Keath et al.,
1989
; Patel et al.,
1998
), were used as control spots on the array. As expected, both
were highly expressed in yeast cells. In addition, CBP1 and
yps-3 were also represented on the array in the set of random genomic
array clones; these spots also showed the same differential expression. To
provide an independent test of whether other clones identified as
differentially expressed in the array analysis were phase-regulated, Northern
analysis was performed on total RNA from yeast and mycelial cells for six of
these clones (Figure 2C). These
clones were recovered from the bacterial archives of the array library, and
the insert corresponding to the genomic DNA fragment on the array was purified
from these clones and used as a probe. Four yeast-specific clones, including
the previously identified CBP1, were confirmed as being primarily
expressed in the yeast form. The enhanced expression of two mycelial-specific
clones was also confirmed. It should be noted that we use "yeast
specific" and "mycelial specific" to refer to quantitatively
different expression levels in the two growth phases.
Annotation of Phase-regulated Genes
Because the Northern analysis confirmed the validity of the array data, the
most highly phase-regulated array clones were annotated. The single-sequence
reads that marked the boundaries of each array clone were mapped onto contigs
from the ongoing genome-sequencing project
(http://genome.wustl.edu/projects/hcapsulatum)
and then a BLASTx analysis of each contig was performed against the National
Center for Biotechnology Information nr database. Clones were annotated if 1)
the entire sequence of the clone encoded a single BLASTx hit, or 2) the
majority of the clone's sequence corresponded to a single BLASTx hit and
flanking sequence was unlikely to correspond to a second gene
(Table 1). Each H.
capsulatum gene was given a three-letter name based on the putative
function of its ortholog. The three-letter code, annotation of each gene,
accession number of its ortholog, and ratio of expression in the two
morphological forms are displayed in Table
1. Differentially expressed clones that contained sequence from
transposable elements or rDNA are excluded from this table.
|
|
Identification of Mycelial-specific Genes
The data revealed several categories of genes showing significantly higher
expression in mycelia compared with yeast. Several orthologs of genes involved
in conidiation (a process confined to mycelia) in other fungi were observed.
For example, an ortholog of fluG, an Aspergillus nidulans
gene that regulates conidial differentiation and secondary metabolite
production was mycelial specific. The H. capsulatum ortholog
FLU1 may be of particular interest because A. nidulans fluG
is required for the accumulation of a presumed extracellular, diffusible
factor that stimulates the differentiation of mycelia into conidia
(Lee and Adams, 1994
). In
A. nidulans, fluG functions as an upstream member of a conidiation
pathway that includes the transcription factor flbC and the regulator
wetA (Marshall and Timberlake,
1991
; Wieser et al.,
1994
). H. capsulatum orthologs of both of these genes
(FBC1 and WET1) were among the set of genes expressed in
mycelia. A formate dehydrogenase homolog (FDH1) was also expressed at
higher levels in mycelia than in yeast. Expression of the N. crassa
formate dehydrogenase is induced under conditions that favor the formation of
conidia (Chow and RajBhandary,
1993
). Because the mycelial sample contained both mycelia and
conidia, it is unknown whether these H. capsulatum orthologs were
expressed in mycelia because they differentiate into conidia, in the conidia
themselves, or both.
Genes that affect mycelial or polar growth in other organisms were also
identified. TYR1 is an ortholog of a tyrosinase gene (MelC2)
from the bacterium Streptomyces griseus. Expression of MelC2
is both necessary and sufficient to produce melanin and stimulate the
production of aerial hyphae (Endo et
al., 2001
). The connection between the production of melanin
and morphology is unclear, but under normal culture conditions, the H.
capsulatum mycelial form produces melanin, whereas yeast cells do not
(Nosanchuk et al.,
2002
). TYR1 is an intriguing candidate for a regulator of
melanin production and establishment or maintenance of the mycelial phase in
H. capsulatum. We also found an ortholog of phenylalanyl-tRNA
synthetase (TRF1), which is required for polar hyphal growth in
A. nidulans (Osherov et
al., 2000
).
Genes that might affect the accumulation of other phase-regulated
transcripts, such as orthologs of the splicing factor PRP8
(Strauss and Guthrie, 1991
)
and the silencing factor HST4, were also induced in mycelia. The
latter is a presumptive NAD+-dependent protein deacetylase that has
roles in silencing in Schizosaccharomyces pombe
(Freeman-Cook et al.,
1999
) as well as other cellular processes in S.
cerevisiae (Brachmann et al.,
1995
). PRP8 and HST4 might influence the overall
mycelial gene expression program by affecting the splicing and expression of
mycelial genes.
Finally, we identified a set of genes that might reflect the unique
physiology of the mycelial form. One of these, NIR1, encodes a
nitrite reductase most similar to the Neisseria gonorrhea aniA gene,
which is required for growth under oxygen-limiting conditions
(Hoehn and Clark, 1992
;
Mellies et al.,
1997
). In other organisms, nitrite reductase functions in
denitrification (Zumft, 1997
),
converting nitrogen oxides into molecular nitrogen. Denitrification is
performed by many soil organisms. Because soil is the natural environment of
the mycelial form of H. capsulatum, the expression of enzymes
required for denitrification may be coupled with growth in the mycelial form.
Similarly, we identified genes that encode several transporters and enzymes
(such as DPP1, which encodes a zinc-regulated diacylglycerol
pyrophosphate phosphatase, and OXO1, which encodes a 3-oxoacyl
reductase) whose differential expression might reflect the different growth
requirements of mycelia and yeast. In addition, we observed two genes
previously identified as mycelial-specific in H. capsulatum: the
MS95 gene (GenBank, Glen Shearer), as well as tubulin, which has been
reported to be expressed fivefold higher in the mycelial phase
(Harris et al.,
1989b
).
Identification of Yeast-specific Genes
Genes with higher expression in the yeast phase than the mycelial phase
were also identified (Table 1),
including several involved in sulfur metabolism. Sulfur metabolism influences
the morphological state of H. capsulatum and other dimorphic fungi
(Medoff et al., 1987
;
Maresca and Kobayashi, 1989
).
In H. capsulatum, at least some strains exhibit a requirement for the
presence of cystine or cysteine in the culture medium to establish the yeast
phase (Salvin, 1949
;
Scherr, 1957
;
Maresca et al.,
1977
). Additionally, a cysteine dioxygenase activity was
previously identified as specific to the yeast phase of H. capsulatum
(Kumar et al., 1983
),
but the gene encoding this protein had not been identified. We identified a
yeast-specific cysteine dioxygenase gene that was expressed 11-fold higher in
yeast than mycelia. In addition, we identified several yeast-expressed genes,
which share sequence similarity to genes involved in sulfur metabolism in
other organisms: choline sulfatase; ATP sulfurylase (the first enzyme in the
sulfate-assimilation arm of the methionine/cysteine biosynthetic pathway);
glutamate-cysteine ligase (which affects glutathione and glutamate
metabolism); and methionine permease (which can mediate both methionine and
cysteine uptake in S. cerevisiae;
Kosugi et al.,
2001
).
The yeast-expressed genes included those that might affect other aspects of metabolism and nutrient availability, such as a lysine permease, an oxoprolinase (which affects L-glutamate production), a 4-hydroxyphenylpyruvate dioxygenase (which affects tyrosine metabolism), a dihydrolipoamide dehydrogenase (which affects carbon and amino acid metabolism), and several ATP-binding cassette and major facilitator superfamily efflux transporters. Expression of these genes predominantly in the yeast form over the mycelial form may reflect different metabolic requirements for the two forms independent of the growth medium. Because the yeast phase represents the form of H. capsulatum found in the host, expression of these genes could reflect the ability of the yeast form to adapt to growth-limiting environments during infection.
Other interesting yeast-regulated genes included those that may affect
aspects of the cell cycle, such as BUB1, which encodes a spindle
assembly checkpoint kinase (Roberts et
al., 1994
), and SMC1, which is involved in sister
chromatid cohesion (Strunnikov et
al., 1993
). The enhanced expression of these genes in the
yeast phase may reflect differences between cell cycle processes in a
unicellular yeast versus a multicellular mycelium. Additionally, there was
markedly increased expression of an ortholog (TRI11) of trichothecene
C-15 hydroxylase in yeast compared with mycelia. The ortholog of this gene
from the plant pathogen Fusarium sporotrichioides encodes a
cytochrome P450 monooxygenase required for production of the toxin
trichothecene (Alexander et al.,
1998
). Trichothecene is thought to inhibit protein synthesis in
many eukaryotes. The role of TRI11 in the virulence of H.
capsulatum yeast is unknown.
Comparison of Yeast and Mycelial Phase-regulated Genes with Genes
Induced during Stationary Phase
The 500 yeast- and mycelial-phaseregulated clones were further
characterized by examining their expression under other conditions. Analysis
of the expression profile of these clones during stationary phase of yeast
cells was useful for the following reasons: Because the mycelial cultures must
be grown for extended periods, a subset of the putative mycelial-expressed
genes might be genes that are induced during stationary phase independent of
morphology. Additionally, because yeast-phase cells transform into mycelia if
kept in stationary phase for extended periods
(Maresca and Kobayashi, 1989
),
genes required for establishment of the mycelial phase might be induced in
stationary-phase yeast cells. In contrast, genes that are required only for
the maintenance of the mycelial phase might show no induction in
stationary-phase yeast cells. A large culture of yeast-phase cells was
inoculated on day 0 and the culture sampled daily from days 2 to 10. The
culture grew exponentially until day 3 and then remained stationary.
Morphologically, the cells remained in the yeast phase throughout the
experiment. Poly-adenylated RNA was used to produce a labeled cDNA probe from
each sample, and the gene expression profile of each time point was determined
relative to the first logarithmic phase time point (day 2).
Cluster analysis, which groups genes that show similar expression profiles
(Eisen et al., 1998
),
was used to determine which of the 500 phase-regulated genes were also
transcriptionally regulated during stationary phase. This analysis revealed a
variety of clusters, several of which are highlighted in
Figure 3. Cluster 1
(Figure 3A) represents clones
that were induced only in mycelial cells and not in stationary-phase yeast
cells. In contrast, cluster 2 represents clones that were induced in both
mycelial cells and stationary-phase yeast cells. By graphing the average fold
change of all clones in each of the highlighted clusters for each time point,
it is apparent that only the second cluster showed significant expression in
stationary phase (Figure
3B).
Cluster 1 includes multiple clones encoding FBC1, NIR1, and TYR1, as well as the MFS4 transporter. The expression of WET1, which falls outside cluster 1, was also restricted to mycelia. Because these genes were not induced in stationary phase, their expression was unlikely to reflect simple nutrient deprivation due to extended growth time. Instead, these genes may be expressed only once cells have switched to mycelial growth.
In contrast, cluster 2 includes clones that were induced in stationary phase of yeast cells as well as in the mycelial form. This cluster includes FLU1. Because the stationary-phase cultures contained only yeast cells and no mycelia or conidia, the expression of FLU1 may precede mycelial conversion and the expression of WET1 and FBC1, which are likely to be involved in production of conidia. The cluster also contains SDH1 and ABC4, which encode sorbitol dehydrogenase and an ATP-binding cassette transporter, respectively.
Some of the yeast-specific clones also showed informative expression patterns during stationary phase. Cluster 3 represents yeast-specific clones that become even more highly expressed in stationary-phase yeast cells compared with mid-logarithmic yeast cells. This cluster includes CHO1, LYP1, TIF3, TRI11, SIP1, and ABC3. Because these genes were induced in stationary-phase yeast cells but were not highly expressed in mycelia, they may reflect nutrient requirements that are specific to yeast-phase cells. Genes such as TRI11 could be induced as the density of the yeast culture increases if quorum sensing regulates production of a toxin.
Because coregulated genes cluster together over a variety of conditions, we
were able to draw some conclusions about the regulation of phase-specific
genes by examining the content of different clusters. Most notably, the
previously identified yeast-specific genes yps-3 and CBP1
fall into two different clusters (Figure 3,
A and B, cluster 4 and 5). Although each of these genes was
significantly more expressed in yeast cells than mycelia, yps-3 was
repressed as yeast cells went into stationary phase, whereas the expression of
CBP1 remained constant as yeast cells aged. This observation was
consistent with the prior suggestion that yps-3 and CBP1 are
subject to different regulation (Rooney
et al., 2001
).
Finally, two other clusters were identified. First, clones that contain rDNA sequence were expressed more strongly in yeast than in mycelia and were repressed in stationary-phase yeast cells. These results likely reflect a difference in metabolic activity among mycelial cells, mid-logarithmic yeast cells, and stationary-phase yeast cells. Second, a group of array clones contains regions of DNA homologous to retrotransposon sequences, which are frequent in the H. capsulatum genome (Mardis, unpublished data). These sequences exhibit increased expression in the yeast form over the mycelial form. It is unclear whether these transposons are active and whether expression of these genes correlates with transposition.
Confirming Differential Expression with Northern Analysis
Northern analysis was used to confirm the phase regulation of a subset of
the genes described above. Although most of the annotated clones clearly
contained a single ORF, a subset of array clones (FLU1, NIR1, WET1, FDH1,
TYR1, MPS1, DPP1, OXO1, GST1, LYP1, ABC1, CHO1, MS95, and ASY1)
contained significant amounts of flanking sequence in addition to the complete
or partial homologous ORF that specified the annotation. For these clones
(with the exception of MS95), the annotation was confirmed by
designing ORF-specific probes for each of the previously mentioned genes.
Northern analysis of yeast and mycelial total RNA confirmed that these genes
were differentially expressed (Figure
4). FLU1, NIR1, WET1, FDH1, TYR1, MPS1, DPP1, OXO1, and
GST1 were more highly expressed in mycelia than in yeast, whereas
LYP1, ABC1, CHO1, and ASY1 were more highly expressed in
yeast than in mycelia.
|
Molecular Basis of Difference in Transcript Size
Surprisingly, Northern analysis revealed that, in addition to their
differential expression, six of the nine mycelial-regulated genes (FLU1,
NIR1, WET1, FDH1, DPP1, and OXO1) exhibited differences in
transcript size in yeast and mycelia. To determine whether this difference was
due to differential transcription initiation sites, poly-adenylation sites, or
splicing, we used 5' RLM-RACE to map the transcription initiation sites
of capped transcripts and 3' RACE to map the sites of poly-adenylation
for FLU1, NIR1, WET1, and FDH1. In addition, RT-PCR was used
to clone the full-length cDNAs from the yeast and mycelial forms. The sequence
of the full-length cDNAs as well as the 5' and 3' RACE products
was compared with the genomic sequence to determine the location of introns.
For all of these genes, we observed major differences in the length of the
5' untranslated region (UTR), likely due to differences in the start
site of transcription (Figure
5B); these differences are reflected by the size of the 5'
RACE products from yeast and mycelia
(Figure 5A). The start site of
the longest yeast transcripts ranged from between 0.8 and 1.5 kb upstream of
the mycelial transcripts. In contrast, TYR1, a mycelial-specific gene
that did not show an altered form in yeast cells by Northern, gave no evidence
of an altered transcript by 5' RACE analysis (Gebhart and Sil,
unpublished data). In addition to altered transcription initiation sites, the
yeast-form transcripts of WET1, FDH1, and NIR1 have introns
in their long 5'UTRs. FLU1, FDH1, and NIR1 had short
5' UTR transcripts in the yeast form in addition to the long 5'
UTR transcripts. In the case of NIR1, this short yeast-form
transcript may encode a cytosolic nitrite reductase rather than the larger
predicted membrane-bound form encoded in mycelia.
|
Although minor differences at the 3' ends of the transcript were revealed by sequencing the 3' RACE clones, these differences contributed to, but could not account for, the large difference in transcript size between yeast and mycelia. Similarly, differential splicing was observed (Figure 5B) but did not account for the large difference in transcript size. Thus, the major mechanism underlying differential transcript size between yeast and mycelia was the use of different transcript initiation sites. These observations suggest that H. capsulatum uses both differential expression and differential usage of transcript initiation site to regulate the transcriptional profile of the two morphological forms.
| DISCUSSION |
|---|
|
|
|---|
We sorted the phase-regulated genes identified herein into broad functional categories (Figure 6). Yeast and mycelial cells differ in morphology, sulfur metabolism, growth rate, growth environment (host vs. soil), melanin production, mating competence, and conidiation. Based on the function of their orthologs (Table 1), we have identified phase-regulated genes that are implicated in several of these processes (Figure 6). Because the ability of yeast cells to survive in the phagosome of host cells may be dependent on nutrient acquisition, we hypothesize that genes that may affect growth rate might also affect host survival (Figure 6). A deeper analysis of the genes identified in this study is necessary before their function can be determined.
|
It may be particularly informative to examine the function of the
yeast-specific genes involved in sulfur metabolism. Although little is
understood about the regulation of the morphological phases, for H.
capsulatum and other systemic dimorphic fungi, the addition of exogenous
sulfhydryl reducing agents (dithiothreitol) to the media traps cells in the
yeast form independent of temperature, whereas the addition of sulfhydryl
oxidizing agents (p-chloromercurphenisulphonic acid) traps cells in
the mycelial form independent of temperature
(Maresca et al.,
1977
; Medoff et al.,
1986
,
1987
)
(Figure 6). Presumably, the
effectors and downstream targets of these regulatory pathways will emerge from
examining gene expression in response to dithiothreitol and
p-chloromercurphenisulfonic acid. Disrupting the genes involved in
sulfur metabolism identified by this work will test the contributions of these
genes to the two growth forms.
We expect that these findings on H. capsulatum will be relevant to
other fungi. There have been only a few published examples of gene expression
profiling by microarray in fungal pathogens (mainly Candida albicans)
(De Backer et al.,
2001
; Kahmann and Basse,
2001
; Murad et al.,
2001
; Lan et al.,
2002
; Lorenz,
2002
; Rogers and Barker,
2002
), and no examples of large-scale microarray studies in the
systemic dimorphic fungi (H. capsulatum, Coccidioides immitis, Blastomyces
dermitiditis, Paracoccidioides brasiliensis, and Penicillium
marneffii). Because alterations in physiology and morphology play key
roles in the pathogenesis of many fungi, the genes and regulatory circuits we
have identified will be relevant to understanding similar switches in other
species. For example, we have identified orthologs of three
Aspergillus regulatory factors required for conidiation:
fluG, which initiates conidiophore development; wetA, whose
expression is both necessary and sufficient for conidia-specific gene
expression; and flbC, which encodes a zinc-finger transcription
factor that is thought to affect the timing of conidiation
(Marshall and Timberlake,
1991
; Lee and Adams,
1994
; Wieser et al.,
1994
). The expression patterns of the FLU1/fluG,
WET1/wetA, and FBC1/flbC orthologs from Histoplasma and
Aspergillus suggest the evolutionary conservation of conidial
development from mycelial cells. This hypothesis is strengthened by the
presence of orthologs of other genes required for conidiation in
Aspergillus, such as brlA and flbA
(Wieser et al., 1994
)
in the H. capsulatum genome. Further annotation of the genes
identified in these experiments will deepen these insights.
In addition to differential gene expression, we observed further regulation
of a subset of phase-regulated genes. Several genes predominantly expressed in
mycelia unexpectedly displayed modest accumulation of transcripts of
significantly altered size in the yeast phase. This phenomenon arose from
different sites of transcript initiation in yeast and mycelia. The underlying
rationale for the production of longer yeast transcripts, some of which
initiated >1 kb upstream of the mycelial transcript, is unclear. Perhaps
the long 5' UTR might be refractory to translation at 37°C but
permissive for translation at 25°C. Such a mechanism might allow the
initial production of mycelial-specific proteins by yeast cells that are
shifted to 25°C until the normal mycelial transcript is produced.
Temperature-dependent regulation of translation of a 5' UTR has been
observed recently in the bacterial pathogen Listeria monocytogenes
(Johansson et al.,
2002
). Translation of transcripts with particular 5' UTRs
can also be influenced by the abundance of translation initiation factors
(Browning et al.,
1988
; Calkhoven et
al., 2000
). Interestingly, we observed two translation
initiation factors, eIF2 and eIF3, that were significantly more expressed in
the yeast form over the mycelial form.
At present, the mechanism of generating the longer 5' UTR transcripts
is unknown. Our observations suggest that an H. capsulatum ortholog
of the SPT6 gene might be significantly expressed in mycelia compared
with yeast. In S. cerevisiae, SPT6 can influence the site of
transcription initiation (Clark-Adams and
Winston, 1987
). Future experiments will be necessary to determine
the prevalence, implications, and mechanism of this regulation of transcript
initiation in H. capsulatum.
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
L.H. was supported by a Damon Runyon Walter Winchell postdoctoral fellowship and a Leukemia and Lymphoma senior fellowship (310703). J.R. is supported by NIHGM35827. A.S. is supported by a Burroughs-Wellcome Molecular Mycology Award (993243), an American Cancer Society Research Scholar Award (RSG-01039-01-MBC), NIAID U01 (UO1-AI-5093401), and the UCSF Fellows Program.
| Footnotes |
|---|
Present address: Department of Organismic and Evolutionary Biology, Harvard
University, Cambridge, MA 02138. ![]()
Corresponding author. E-mail address:
sil{at}cgl.ucsf.edu.
| REFERENCES |
|---|
|
|
|---|
Brachmann, C.B., Sherman, J.M., Devine, S.E., Cameron, E.E.,
Pillus, L., and Boeke, J.D. (1995). The SIR2 gene
family, conserved from bacteria to humans, functions in silencing, cell cycle
progression, and chromosome stability. Genes Dev.
9,
2888-2902.
Browning, K.S., Lax, S.R., Humphreys, J., Ravel, J.M., Jobling,
S.A., and Gehrke, L. (1988). Evidence that the
5'-untranslated leader of mRNA affects the requirement for wheat germ
initiation factors 4A, 4F, and 4G. J. Biol. Chem.
263,
9630-9634.
Bullock, W.E. (1993). Interactions between human phagocytic cells and Histoplasma capsulatum. Arch. Med. Res. 24, 219-223.[Medline]
Calkhoven, C.F., Muller, C., and Leutz, A. (2000).
Translational control of C/EBP
and C/EBP
isoform expression.
Genes Dev. 14,
1920-1932.
Carr, J., and Shearer, G., Jr. (1998). Genome size,
complexity, and ploidy of the pathogenic fungus Histoplasma
capsulatum. J. Bacteriol.
180,
6697-6703.
Chow, C.M., and RajBhandary, U.L. (1993).
Developmental regulation of the gene for formate dehydrogenase in
Neurospora crassa. J. Bacteriol.
175,
3703-3709.
Clark-Adams, C.D., and Winston, F. (1987). The
SPT6 gene is essential for growth and is required for
-mediated transcription in Saccharomyces cerevisiae.
Mol. Cell. Biol. 7,
679-686.
De Backer, M.D., Ilyina, T., Ma, X.J., Vandoninck, S., Luyten,
W.H., and Vanden Bossche, H. (2001). Genomic profiling of the
response of Candida albicans to itraconazole treatment using a DNA
microarray. Antimicrob. Agents Chemother.
45,
1660-1670.
De Marie, S. (2000). New developments in the diagnosis
and management of invasive fungal infections. Haematologica
85,
88-93.
De Pauw, B.E. (2001). Treatment of documented and suspected neutropenia-associated invasive fungal infections. J. Chemother. 13, 181-192.
DeRisi, J.L., Iyer, V.R., and Brown, P.O. (1997).
Exploring the metabolic and genetic control of gene expression on a genomic
scale. Science 278,
680-686.
Di Lallo, G., Gargano, S., and Maresca, B. (1994). The Histoplasma capsulatum cdc2 gene is transcriptionally regulated during the morphologic transition. Gene 140, 51-57.[CrossRef][Medline]
Dixon, D.M., McNeil, M.M., Cohen, M.L., Gellin, B.G., and La Montagne, J.R. (1996). Fungal infections: a growing threat. Public Health Rep. 111, 226-235.[Medline]
Eisen, M.B., Spellman, P.T., Brown, P.O., and Botstein, D.
(1998). Cluster analysis and display of genome-wide expression
patterns. Proc. Natl. Acad. Sci. USA
95,
14863-14868.
Eissenberg, L.G., and Goldman, W.E. (1991).
Histoplasma variation and adaptive strategies for parasitism: new
perspectives on histoplasmosis. Clin. Microbiol. Rev.
4,
411-421.
Eissenberg, L.G., and Goldman, W.E. (1994). The interplay between Histoplasma capsulatum and its host cells. Balliere's Clin. Infect. Dis. 1, 265-283.
Endo, K., Kamo, K., Hosono, K., Beppu, T., and Ueda, K. (2001). Characterization of mutants defective in melanogenesis and a gene for tyrosinase of Streptomyces griseus. J. Antibiot. 54, 789-796.[Medline]
Freeman-Cook, L.L., Sherman, J.M., Brachmann, C.B., Allshire, R.C.,
Boeke, J.D., and Pillus, L. (1999). The Schizosaccharomyces
pombe hst4(+) gene is a SIR2 homologue with silencing and
centromeric functions. Mol. Biol. Cell
10,
3171-3186.
Gargano, S., Di Lallo, G., Kobayashi, G.S., and Maresca, B.
(1995). A temperature-sensitive strain of Histoplasma
capsulatum has an altered
9-fatty acid desaturase gene.
Lipids 30,
899-906.[Medline]
Gish, W., and States, D.J. (1993). Identification of protein coding regions by database similarity search. Nat. Genet. 3, 266-272.[CrossRef][Medline]
Harris, G.S., Keath, E.J., and Medoff, J. (1989a). Characterization of alpha and beta tubulin genes in the dimorphic fungus Histoplasma capsulatum. J. Gen. Microbiol. 135, 1817-1832.[Medline]
Harris, G.S., Keath, E.J., and Medoff, J. (1989b).
Expression of
- and
-tubulin genes during dimorphic-phase
transitions of Histoplasma capsulatum. Mol. Cell.
Biol., 9,
2042-2049.
Hoehn, G.T., and Clark, V.L. (1992). Isolation and
nucleotide sequence of the gene (aniA) encoding the major
anaerobically induced outer membrane protein of Neisseria
gonorrhoeae. Infect Immun.
60,
4695-4703.
Johansson, J., Mandin, P., Renzoni, A., Chiaruttini, C., Springer, M., and Cossart, P. (2002). An RNA thermosensor controls expression of virulence genes in Listeria monocytogenes. Cell 110, 551-561.[CrossRef][Medline]
Johnson, C.H., Klotz, M.G., York, J.L., Kruft, V., and McEwen, J.E.
(2002). Redundancy, phylogeny and differential expression of
Histoplasma capsulatum catalases. Microbiology
148,
1129-1142.
Kahmann, R., and Basse, C. (2001). Fungal gene expression during pathogenesis-related development and host plant colonization. Curr. Opin. Microbiol. 4, 374-380.[CrossRef][Medline]
Keath, E.J., and Abidi, F.E. (1994). Molecular cloning and sequence analysis of yps-3, a yeast-phase-specific gene in the dimorphic fungal pathogen Histoplasma capsulatum. Microbiology 140, 759-767.[Abstract]
Keath, E.J., Painter, A.A., Kobayashi, G.S., and Medoff, G.
(1989). Variable expression of a yeast-phase-specific gene in
Histoplasma capsulatum strains differing in thermotolerance and virulence.
Infect. Immun. 57,
1384-1390.
Kosugi, A., Koizumi, Y., Yanagida, F., and Udaka, S. (2001). MUP1, high affinity methionine permease, is involved in cysteine uptake by Saccharomyces cerevisiae. Biosci. Biotechnol. Biochem. 65, 728-731.[CrossRef][Medline]
Kumar, V., Maresca, B., Sacco, M., Goewert, R., Kobayashi, G.S., and Medoff, G. (1983). Purification and characterization of a cysteine dioxygenase from the yeast phase of Histoplasma capsulatum. Biochemistry 22, 762-768.[CrossRef][Medline]
Lan, C.Y., Newport, G., Murillo, L.A., Jones, T., Scherer, S.,
Davis, R.W., and Agabian, N. (2002). Metabolic specialization
associated with phenotypic switching in Candida albicans.
Proc. Natl. Acad. Sci. USA 99,
14907-14912.
Lee, B.N., and Adams, T.H. (1994). The Aspergillus
nidulans fluG gene is required for production of an extracellular
developmental signal and is related to prokaryotic glutamine synthetase I.
Genes Dev. 8,
641-651.
Lorenz, M.C. (2002). Genomic approaches to fungal pathogenicity. Curr. Opin. Microbiol. 5, 372-378.[CrossRef][Medline]
Maresca, B., Carratu, L., and Kobayashi, G.S. (1994). Morphological transition in the human fungal pathogen Histoplasma capsulatum. Trends Microbiol. 2, 110-114.[CrossRef][Medline]
Maresca, B., and Kobayashi, G.S. (1989). Dimorphism in
Histoplasma capsulatum: a model for the study of cell differentiation
in pathogenic fungi. Microbiol. Rev.
53,
186-209.
Maresca, B., Medoff, G., Schlessinger, D., and Kobayashi, G.S. (1977). Regulation of dimorphism in the pathogenic fungus Histoplasma capsulatum. Nature 266, 447-448.[CrossRef][Medline]
Marshall, M.A., and Timberlake, W.E. (1991).
Aspergillus nidulans wetA activates spore-specific gene expression.
Mol. Cell. Biol. 11,
55-62.
Medoff, G., Painter, A., and Kobayashi, G.S. (1987).
Mycelial- to yeast-phase transitions of the dimorphic fungi Blastomyces
dermatitidis and Paracoccidioides brasiliensis. J.
Bacteriol. 169,
4055-4060.
Medoff, G., Sacco, M., Maresca, B., Schlessinger, D., Painter, A.,
Kobayashi, G.S., and Carratu, L. (1986). Irreversible block of
the mycelial-to-yeast phase transition of Histoplasma capsulatum.
Science 231,
476-479.
Mellies, J., Jose, J., and Meyer, T.F. (1997). The Neisseria gonorrhoeae gene aniA encodes an inducible nitrite reductase. Mol. Gen. Genet. 256, 525-532.[CrossRef][Medline]
Momin, F., and Chandrasekar, P.H. (1995).
Antimicrobial prophylaxis in bone marrow transplantation. Ann. Intern.
Med. 123,
205-215.
Murad, A.M., d'Enfert, C., Gaillardin, C., Tournu, H., Tekaia, F., Talibi, D., Marechal, D., Marchais, V., Cottin, J., and Brown, A.J. (2001). Transcript profiling in Candida albicans reveals new cellular functions for the transcriptional repressors CaTup1, CaMig1 and CaNrg1. Mol. Microbiol. 42, 981-993.[CrossRef][Medline]
Nosanchuk, J.D., Gomez, B.L., Youngchim, S., Diez, S., Aisen, P.,
Zancope-Oliveira, R.M., Restrepo, A., Casadevall, A., and Hamilton, A.J.
(2002). Histoplasma capsulatum synthesizes melanin-like
pigments in vitro and during mammalian infection. Infect.
Immun. 70,
5124-5131.
Osherov, N., Mathew, J., and May, G.S. (2000). Polarity-defective mutants of Aspergillus nidulans. Fungal Genet. Biol. 31, 181-188.[CrossRef][Medline]
Patel, J.B., Batanghari, J.W., and Goldman, W.E.
(1998). Probing the yeast phase-specific expression of the
CBP1 gene in Histoplasma capsulatum. J.
Bacteriol. 180,
1786-1792.
Roberts, B.T., Farr, K.A., and Hoyt, M.A. (1994). The
Saccharomyces cerevisiae checkpoint gene BUB1 encodes a
novel protein kinase. Mol. Cell. Biol.
14,
8282-8291.
Rogers, P.D., and Barker, K.S. (2002). Evaluation of
differential gene expression in fluconazole-susceptible and -resistant
isolates of Candida albicans by cDNA microarray analysis.
Antimicrob. Agents Chemother.
46,
3412-3417.
Rooney, P.J., Sullivan, T.D., and Klein, B.S. (2001). Selective expression of the virulence factor BAD1 upon morphogenesis to the pathogenic yeast form of Blastomyces dermatitidis: evidence for transcriptional regulation by a conserved mechanism. Mol. Microbiol. 39, 875-889.[CrossRef][Medline]
Salvin, S.B. (1949). Cysteine and related compounds in the growth of the yeast-like phase of H. capsulatum. J. Infect. Dis., 84, 275-283.[Medline]
Samonis, G., and Bafaloukos, D. (1992). Fungal infections in cancer patients: an escalating problem. In Vivo 6, 183-193.[Medline]
Scherr, G.H. (1957). The roles of SH group and incubation temperature. Exp. Cell Res. 12, 92-102.[CrossRef][Medline]
Sebghati, T.S., Engle, J.T., and Goldman, W.E. (2000).
Intracellular parasitism by Histoplasma capsulatum: fungal virulence
and calcium dependence. Science
290,
1368-1372.
Strauss, E.J., and Guthrie, C. (1991). A
cold-sensitive mRNA splicing mutant is a member of the RNA helicase gene
family. Genes Dev. 5,
629-641.
Strunnikov, A.V., Larionov, V.L., and Koshland, D.
(1993). SMC1:an essential yeast gene encoding a putative
head-rod-tail protein is required for nuclear division and defines a new
ubiquitous protein family. J. Cell Biol.
123,
1635-1648.