- Research article
- Open Access
Differential expression of lipoprotein genes in Mycoplasma pneumoniae after contact with human lung epithelial cells, and under oxidative and acidic stress
© Hallamaa et al; licensee BioMed Central Ltd. 2008
- Received: 26 February 2008
- Accepted: 23 July 2008
- Published: 23 July 2008
Mycoplasma pneumoniae is a human pathogen that is a common cause of community-acquired pneumonia. It harbours a large number of lipoprotein genes, most of which are of unknown function. Because of their location on the cell surface, these proteins are likely to be involved in the bacterial response to environmental changes, or in the initial stages of infection. The aim of this study was to determine if genes encoding surface lipoproteins are differentially expressed after contact with a human cell line, or after exposure to oxidative or acidic stress.
Using qRT-PCR assays, we observed that the expression of a number of lipoprotein genes was up-regulated when M. pneumoniae was placed in contact with human cells. In contrast, lipoprotein expression was generally down-regulated or unchanged when exposed to either hydrogen peroxide or low pH (5.5). When exposed to low pH, the mRNA levels of four polycistronically transcribed genes in Lipoprotein Multigene Family 6 formed a gradient of decreasing quantity with increasing distance from a predicted promoter.
The demonstrated transcriptional changes provide evidence for the functionality of these mostly unassigned genes and indicate that they are regulated in response to changes in environmental conditions. In addition we have shown that the members of Lipoprotein Gene Family 6 may be expressed polycistronically.
Mycoplasma pneumoniae is one of the most common causes of community acquired pneumonia . It is thought to be an extracellular pathogen that colonizes the surface of epithelial cells, although there are reports suggesting that it is capable of invading host cells [2, 3]. Mycoplasmas lack a cell wall, and, hence, their membrane-bound proteins are exposed to the external environment and to their host's immune response. In spite of their apparent vulnerability to immune responses, they establish chronic infections.
The M. pneumoniae genome has a high proportion of genes predicted to encode lipoproteins (6.68% of all genes) and 67 of these are members of six lipoprotein gene families [4, 5]. A full description of these gene families, including their genomic locations, has been published previously . The functionality of these genes is unknown, but because their products are cell surface proteins they could potentially play a role in adhesion and/or the initial stages of infection. The members of Lipoprotein Gene Family 2 are homologs of an M. hyopneumoniae lipoprotein located immediately upstream of a conserved ABC transport operon containing an exonuclease . Little is known about general stress responses in M. pneumoniae, with only one study exploring the effects of heat shock , and another reporting changes in the global protein synthesis pattern in response to the presence of glycerol . It has been suggested that the genome lacks many of the transcriptional regulators found in other bacteria. For example, no two-component regulatory system and only one sigma factor has been identified [4, 9]. Hence, the extent of transcriptional regulation and the mechanisms involved remain largely unknown in this organism.
Oligonucleotides used for qPCR
Gene annotation (NCBI)
Bacterial strain and culture conditions
M. pneumoniae strain M129 was cultured in 150 cm2 tissue culture flasks (Falcon, Becton Dickinson) in Modified Hayflick Medium  supplemented with 20% horse serum (Sigma H-1138) at 37°C until a colour change was seen in the medium (1–2 days).
The human cell line and culture conditions
The human lung cell line A549 was initially grown in flasks at 37°C in Minimum Essential Medium (MEM) supplemented with 5% fetal calf serum (FCS), 5 mM sodium hydrogen carbonate, 10 mM HEPES (pH 7.4) and ampicillin (50 μg/ml). The cells were subsequently cultured in 6-well tissue culture plates at 37°C in Dulbecco's Modified Eagle's Medium (DMEM) supplemented with 5% FCS, 44 mM sodium hydrogen carbonate and ampicillin (50 μg/ml) in 5% CO2/air, or in Leibovitz's L15 medium supplemented with 5% FCS, 10 mM HEPES (pH 7.4) and ampicillin (50 μg/ml) in air.
M. pneumoniae M129 cultured in a 150 cm2 cell culture flask were scraped from the plastic into MEM supplemented with 100 U penicillin/ml. The suspension was passed through a 27-gauge needle to separate the cells.
Two millilitres of fresh MEM containing 5% FCS was added to the wells containing confluent cultures of A549 cells and 250 μl of the mycoplasma suspension was added to each well (an approximate multiplicity of infection of 100). As a control, the mycoplasma suspension was added to a well containing MEM supplemented with 5% FCS but not cells. The plates were incubated at 37°C for 15 min, 30 min, 1 h or 2 h. All time points were chosen to be less than the doubling time for M. pneumoniae to ensure that replication did not affect interpretation of the assay. Following incubation the medium was removed and the wells were washed twice with 2 ml of PBS at 37°C to remove non-adherent mycoplasmas. A 600 μl volume of RLT buffer (RNeasy kit, Qiagen) was added to each well and the cells were collected using a cell scraper. The cell lysate was passed through a 20-gauge needle at least 5 times and RNA was isolated using a Qiagen RNeasy kit as per the manufacturer's instructions. RNA was eluted in 50 μl of DEPC-treated water and stored at -70°C. Three independent experiments were carried out, each containing three biological replicates.
M. pneumoniae M129 was cultured in 25 cm2 cell culture flasks. The cells were then incubated for 10 min in Modified Hayflick Medium at 37°C supplemented with 1 mM or 20 mM H2O2. No H2O2 was added to the negative control. The medium was discarded and RNA was isolated as described above. This experiment was repeated four times.
Low pH stress test
M. pneumoniae M129 was cultured in 25 cm2 cell culture flasks. The cells were then incubated for 10 min in Modified Hayflick Medium, with the pH adjusted to 5.5 or 7.4 (control) at 37°C. The medium was discarded and RNA was isolated as described above. This experiment was repeated three times.
DNase treatment of RNA
RNA was treated with DNase prior to RT-PCR. A 15 μl aliquot of RNA was treated with 2 U of RNase-free DNase I (Novagen) in a final volume of 25 μl at 37°C for 30 minutes. The reaction was stopped by the addition of 50 mM EDTA to a final concentration of 4.5 mM, followed by incubation at 70°C for 5 minutes.
Reverse transcription was carried out by incubating 0.5 μl random primers (3 μg/μl, Invitrogen), 1 μl dNTPs (10 nM each) and 10 μl DNase-treated RNA at 65°C for 10 minutes, and then adding 2.5 μl 10× reaction buffer (F-570B, Finnzymes), 0.5 μl RNase inhibitor (RNAguard, Amersham Biosciences), 0.5 μl AMV-RT (F-570S, Finnzymes) and DEPC-treated water to a final volume of 25 μl. The reaction was incubated for 10 minutes at 25°C, 30 min at 42°C, and then 2 min at 94°C. The cDNA was diluted 1:10 with distilled water and used as template for qPCR. As a control for detecting possible DNA contamination in samples, a duplicate reverse transcription reaction for each sample was performed without the addition of AMV-RT. The control reactions were treated identically to samples in the subsequent dilution step and in qPCR.
The qPCR reactions were carried out using a Stratagene MX3000P® System and Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen) chemistry. Expression ratios were calculated using REST 2005 . Reactions consisted of 1× Platinum SYBR Green qPCR SuperMix-UDG containing the ROX reference dye at a final concentration of 50 nM, 0.75 μM forward primer, 0.75 μM reverse primer, 5 μl of template and distilled water to a final volume of 20 μl. The reactions were incubated at 50°C for 2 min, then 95°C for 5 min, then 40 cycles of 95°C for 30 s and 60°C for 30 s. Melting curve analysis was performed at the end of the amplification.
The qPCR efficiency for each primer pair was determined by running a tenfold dilution series in triplicate using M129 genomic DNA as template, at concentrations ranging from 0.5 pg to 5 ng per reaction, covering a cycle threshold (CT) range from approximately 15 to 30. The specificity of each primer pair was confirmed by melting curve analysis and subsequent agarose gel electrophoresis of qPCR products. All primer pairs produced a PCR product with a single peak in melting curve analyses and a single band after gel electrophoresis. Oligonucleotide primer sequences used in this study are provided in Table 1.
Three housekeeping genes, for glyceraldehyde-3-phosphate dehydrogenase (gap, also known as MPN430), elongation factor Tu (tuf, also known as MPN665) and DNA gyrase subunit B (gyrB, also known as MPN003), were validated in an adhesion assay by ensuring that the relative expression of each gene remained constant before and after contact of M. pneumoniae M129 cells with A549 cells. Two housekeeping genes were used for normalisation when analysing H2O2-treated samples, glyceraldehyde-3-phosphate dehydrogenase (gap or MPN430) and elongation factor Tu (tuf or MPN665). The expression of gyrB was not stable after exposure to H2O2. The expression of all three housekeeping genes was stable after exposure to low pH, so all three genes were used for normalisation of low pH stress samples.
Statistical analysis of qPCR data was carried out using the randomisation test included in REST 2005 . A P value of less than 0.05 was considered significant.
Several lipoprotein genes were up-regulated after mycoplasma contact with human A549 cells
Changes in M. pneumoniae lipoprotein gene expression after contact with human A549 cells for 2 hours, or after incubation in 20 mM H2O2 or at pH 5.5 for 10 min
Gene annotation (NCBI)
Change in expression (fold) after:
Contact with A549 cells#
20 mM H2O2 for 10 min†
pH 5.5 for 10 min§
Two lipoprotein genes were down-regulated in response to hydrogen peroxide
Two lipoprotein genes, MPN489 in Family 3 and MPN639 in Family 6, were down-regulated when mycoplasma cells were incubated in growth medium supplemented with 20 mM H2O2 (Table 2). In addition, several genes in Family 2, MPN408, MPN200 and MPN152, MPN436 in Family 3, MPN012 and MPN505 in Family 5, and MPN646, MPN647 and MPN641 in Family 6 appeared to be down-regulated, although the results were not statistically significant. The lower concentration of H2O2 (1 mM) had no significant effect on lipoprotein gene expression (see Additional file 2: Change in lipoprotein gene expression in M. pneumoniae after exposure to 1 mM H2O2 for 10 minutes).
Three lipoprotein genes were down-regulated after acid stress
Three genes in Family 2, MPN408, MPN200 and MPN152, were significantly down-regulated after the cells were exposed to pH 5.5 for ten minutes (Table 2). Interestingly, after low pH stress, the relative mRNA concentrations of four consecutive genes in Family 6, MPN647, MPN646, MPN645 and MPN644, formed a decreasing gradient (R2 = 0.997). By employing oligonucleotides bridging adjacent ORFs in RT-PCR, these genes were shown to be polycistronically expressed (see Additional file 3: Four consecutive genes in Family 6, MPN647, MPN646, MPN645 and MPN644 are polycistronically expressed as shown by RT-PCR).
Family 2 had the most diverse expression profiles under different environmental conditions
MPN200 in Family 2 was the only gene affected by more than one of the environmental conditions tested. This gene was up-regulated after contact with A549 cells and down-regulated after exposure to low pH.
In this study we report transcriptional changes in the expression of M. pneumoniae lipoprotein genes in response to contact with human epithelial cells in vitro, and provide further evidence of transcriptional regulation in mycoplasmas in response to environmental conditions [7, 17–19]. We chose qRT-PCR as the method to study changes in gene expression. Although microarray and qRT-PCR results generally correlate well, both oligonucleotide and cDNA arrays have a tendency to underestimate the fold-change ratios of the underlying mRNAs compared to qRT-PCR . Thus, qRT-PCR is more sensitive in detecting smaller changes in gene expression .
The hydrogen peroxide concentrations used in this study were chosen to have a minimal impact on cell viability, based on unpublished studies by Zimmermann and Herrmann . A pH of 5.5 was used in the acidic stress test and was chosen to be below the pH of exhaled breath concentrate (pH 6.4)  but above the pH of 4.8 found in lysosomes. In this study mycoplasma cells were exposed to low pH or hydrogen peroxide for 10 minutes before changes in gene expression were measured. Similar hydrogen peroxide treatment times have been used previously in gene expression studies of M. pneumoniae [7, 20] and Haemophilus influenzae . In our cell culture infection model we allowed the mycoplasma cells to be in contact with eukaryotic cells for 15, 30, 60 and 120 minutes. This timeframe was chosen to investigate changes in gene expression during the initial stages of infection, but prior to any replication.
Global transcriptional changes in M. pneumoniae after heat shock have been studied using microarray technology . Forty-seven genes were up-regulated after heat shock at 43°C, while 30 genes were down-regulated in the same conditions. Of the genes analysed in our study, only three genes showed differential expression in response to heat shock; MPN199 and MPN200 in Family 2 were up-regulated, whereas MPN591 in Family 1 was down-regulated. This is in agreement with our observation that Family 2 genes appear to be regulated in response to a variety of environmental conditions. Recently, in a microarray study of transcriptional regulation in M. gallisepticum associated with eukaryotic cells, 58 genes were found to be either up- or down-regulated . The expression of M. gallisepticum genes homologous to M. pneumoniae lipoproteins in multigene Families 1, 2, 3, and 4 was not affected by contact with eukaryotic cells as detected by microarray.
All genes, apart from MPN456, that were differentially expressed in this study are hypothetical genes with no assigned function. The results from this study imply that these genes could be functional and that their gene products may play a role in adaptation to different environmental conditions. Two genes that belong to a gene family unique to M. pneumoniae, MPN011 and MPN411, were up-regulated when the cells were in contact with A549 cells, suggesting that the genes in this family could be involved in adhesion or in the initial stages of infection. MPN456, a homolog of a predicted component of a predicted oligopeptide ABC transport system (COG4166), was also up-regulated after adhesion to A549 cells. This could indicate regulation in response to a requirement for transport of substrates from the host cell or involvement in adhesion, as in M. hominis, in which a substrate-binding domain of a peptide transport system (OppA) mediates cytadherence [23, 24]. However, further studies are needed to understand the role of these genes in vivo and whether the changes observed in gene expression of M. pneumoniae during association with human cells are due to adhesion, invasion or other interactions with the cell, although some recent reports of the ability of M. pneumoniae to invade human cells have been contradictory [2, 3, 25].
After acidic stress the mRNA of four genes in Family 6 formed a decreasing concentration gradient, in correlation with their distance from a predicted promoter . By employing RT-PCR and oligonucleotides bridging the adjacent coding sequences, these genes were shown to be expressed polycistronically, contrary to our previous findings . A positive correlation between mRNA levels and proximity to a promoter was also observed by Benders et al.  when they measured mRNA levels of the ftsZ gene cluster in M. pneumoniae. Groups of functionally related genes in M. pneumoniae and M. genitalium are often preceded by promoters. However no transcriptional terminators are located downstream of these genes to prevent follow-through of transcription to subsequent ORFs . It has been proposed that there are only a few transcriptional terminators in the mycoplasma genomes [28, 29] and, therefore, that the transcription of many genes may be initiated from upstream promoters located in front of other genes, not all of which are necessarily functionally related to those downstream.
In summary, we have assessed transcriptional responses of M. pneumoniae after association with host cells and after exposure to hydrogen peroxide or low pH. A 20 mM concentration of H2O2 induced a trend towards general down-regulation of lipoprotein genes, whereas low pH had a more diverse effect on lipoprotein gene expression. Adhesion of mycoplasma cells to A549 cells resulted in up-regulation of selected genes, but the expression of most genes remained unaffected. These observations of transcriptional changes provide evidence for the functionality of these mostly unassigned genes and indicate that they may have a role in responses by M. pneumoniae to environmental changes.
This work was supported by an NHMRC project grant.
- Foy HM: Infections caused by Mycoplasma pneumoniae and possible carrier state in different populations of patients. Clin Infect Dis. 1993, 17 Suppl 1: S37-46.View ArticlePubMedGoogle Scholar
- Meseguer MA, Alvarez A, Rejas MT, Sanchez C, Perez-Diaz JC, Baquero F: Mycoplasma pneumoniae: a reduced-genome intracellular bacterial pathogen. Infect Genet Evol. 2003, 3: 47-55. 10.1016/S1567-1348(02)00151-X.View ArticlePubMedGoogle Scholar
- Yavlovich A, Tarshis M, Rottem S: Internalization and intracellular survival of Mycoplasma pneumoniae by non-phagocytic cells. FEMS Microbiol Lett. 2004, 233: 241-246. 10.1111/j.1574-6968.2004.tb09488.x.View ArticlePubMedGoogle Scholar
- Dandekar T, Huynen M, Regula JT, Ueberle B, Zimmermann CU, Andrade MA, Doerks T, Sánchez-Pulido L, Snel B, Suyama M, Yuan YP, Herrmann R, Bork P: Re-annotating the Mycoplasma pneumoniae genome sequence: adding value, function and reading frames. Nucleic Acids Res. 2000, 28: 3278-3288. 10.1093/nar/28.17.3278.PubMed CentralView ArticlePubMedGoogle Scholar
- Hallamaa KM, Browning GF, Tang SL: Lipoprotein multigene families in Mycoplasma pneumoniae. J Bacteriol. 2006, 188: 5393-5399. 10.1128/JB.01819-05.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmidt JA, Browning GF, Markham PF: Mycoplasma hyopneumoniae mhp379 is a Ca2+-dependent, sugar-nonspecific exonuclease exposed on the cell surface. J Bacteriol. 2007, 189: 3414-3424. 10.1128/JB.01835-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Weiner J, Zimmerman CU, Göhlmann HW, Herrmann R: Transcription profiles of the bacterium Mycoplasma pneumoniae grown at different temperatures. Nucleic Acids Res. 2003, 31: 6306-6320. 10.1093/nar/gkg841.PubMed CentralView ArticlePubMedGoogle Scholar
- Halbedel S, Eilers H, Jonas B, Busse J, Hecker M, Engelmann S, Stülke J: Transcription in Mycoplasma pneumoniae: analysis of the promoters of the ackA and ldh Genes. J Mol Biol. 2007, 371: 596-607. 10.1016/j.jmb.2007.05.098.View ArticlePubMedGoogle Scholar
- Himmelreich R, Hilbert H, Plagens H, Pirkl E, Li BC, Herrmann R: Complete sequence analysis of the genome of the bacterium Mycoplasma pneumoniae. Nucleic Acids Res. 1996, 24: 4420-4449. 10.1093/nar/24.22.4420.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayakawa M, Taguchi H, Kamiya S, Fujioka Y, Watanabe H, Kawai S, Kobayashi H: Animal model of Mycoplasma pneumoniae infection using germfree mice. Clin Diagn Lab Immunol. 2002, 9: 669-676. 10.1128/CDLI.9.3.669-676.2002.PubMed CentralPubMedGoogle Scholar
- Svenstrup HF, Nielsen PK, Drasbek M, Birkelund S, Christiansen G: Adhesion and inhibition assay of Mycoplasma genitalium and M. pneumoniae by immunofluorescence microscopy. J Med Microbiol. 2002, 51: 361-373. 10.1016/S0167-7012(02)00113-6.View ArticlePubMedGoogle Scholar
- Yang J, Hooper WC, Phillips DJ, Talkington DF: Regulation of proinflammatory cytokines in human lung epithelial cells infected with Mycoplasma pneumoniae. Infect Immun. 2002, 70: 3649-3655. 10.1128/IAI.70.7.3649-3655.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Dallas PB, Gottardo NG, Firth MJ, Beesley AH, Hoffmann K, Terry PA, Freitas JR, Boag JM, Cummings AJ, Kees UR: Gene expression levels assessed by oligonucleotide microarray analysis and quantitative real-time RT-PCR – how well do they correlate?. BMC Genomics. 2005, 6: 59-10.1186/1471-2164-6-59.PubMed CentralView ArticlePubMedGoogle Scholar
- Yuen T, Wurmbach E, Pfeffer RL, Ebersole BJ, Sealfon SC: Accuracy and calibration of commercial oligonucleotide and custom cDNA microarrays. Nucleic Acids Res. 2002, 30: e48-10.1093/nar/30.10.e48.PubMed CentralView ArticlePubMedGoogle Scholar
- Freundt EA: A collaborative assay of mycoplasma reference antisera. J Biol Stand. 1983, 11: 227-240. 10.1016/S0092-1157(83)80010-9.View ArticlePubMedGoogle Scholar
- Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.PubMed CentralView ArticlePubMedGoogle Scholar
- Cecchini KR, Gorton TS, Geary SJ: Transcriptional responses of Mycoplasma gallisepticum strain R in association with eukaryotic cells. J Bacteriol. 2007, 189: 5803-5807. 10.1128/JB.00667-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Madsen ML, Nettleton D, Thacker EL, Edwards R, Minion FC: Transcriptional profiling of Mycoplasma hyopneumoniae during heat shock using microarrays. Infect Immun. 2006, 74: 160-166. 10.1128/IAI.74.1.160-166.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Madsen ML, Nettleton D, Thacker EL, Minion FC: Transcriptional profiling of Mycoplasma hyopneumoniae during iron depletion using microarrays. Microbiology. 2006, 152: 937-944. 10.1099/mic.0.28674-0.View ArticlePubMedGoogle Scholar
- Zimmerman C-UR, Herrmann R: Transcriptome analysis of Mycoplasma pneumoniae Subtypes M129 and FH During Hydrogen Peroxide Induced Oxidative Stress. International Organisation for Mycoplasmology 16th International Congress. Edited by: Ayling RD, Citti C, Nicholas RAJ. 2006, Cambridge, UK, 99-Google Scholar
- Czebe K, Kullmann T, Csiszer E, Barat E, Horvath I, Antus B: Variability of exhaled breath condensate pH in lung transplant recipients. Respiration. 2007, 75 (3): 322-7.Google Scholar
- Harrison A, Ray WC, Baker BD, Armbruster DW, Bakaletz LO, Munson RS: The OxyR regulon in nontypeable Haemophilus influenzae. J Bacteriol. 2007, 189: 1004-1012. 10.1128/JB.01040-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Henrich B, Hopfe M, Kitzerow A, Hadding U: The adherence-associated lipoprotein P100, encoded by an opp operon structure, functions as the oligopeptide-binding domain OppA of a putative oligopeptide transport system in Mycoplasma hominis. J Bacteriol. 1999, 181: 4873-4878.PubMed CentralPubMedGoogle Scholar
- Hopfe M, Henrich B: OppA, the substrate-binding subunit of the oligopeptide permease, is the major Ecto-ATPase of Mycoplasma hominis. J Bacteriol. 2004, 186: 1021-1928. 10.1128/JB.186.4.1021-1028.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Krunkosky TM, Jordan JL, Chambers E, Krause DC: Mycoplasma pneumoniae host-pathogen studies in an air-liquid culture of differentiated human airway epithelial cells. Microb Pathog. 2007, 42: 98-103. 10.1016/j.micpath.2006.11.003.View ArticlePubMedGoogle Scholar
- Weiner J, Herrmann R, Browning GF: Transcription in Mycoplasma pneumoniae. Nucleic Acids Res. 2000, 28: 4488-4496. 10.1093/nar/28.22.4488.PubMed CentralView ArticlePubMedGoogle Scholar
- Benders GA, Powell BC, Hutchison CA: Transcriptional analysis of the conserved ftsZ gene cluster in Mycoplasma genitalium and Mycoplasma pneumoniae. J Bacteriol. 2005, 187: 4542-4551. 10.1128/JB.187.13.4542-4551.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Ermolaeva MD, Khalak HG, White O, Smith HO, Salzberg SL: Prediction of transcription terminators in bacterial genomes. J Mol Biol. 2000, 301: 27-33. 10.1006/jmbi.2000.3836.View ArticlePubMedGoogle Scholar
- Washio T, Sasayama J, Tomita M: Analysis of complete genomes suggests that many prokaryotes do not rely on hairpin formation in transcription termination. Nucleic Acids Res. 1998, 26: 5456-5463. 10.1093/nar/26.23.5456.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.