- Research article
- Open Access
- Open Peer Review
SOD2polymorphisms: unmasking the effect of polymorphism on splicing
© Shao et al; licensee BioMed Central Ltd. 2007
- Received: 23 August 2006
- Accepted: 01 March 2007
- Published: 01 March 2007
The SOD2 gene encodes an antioxidant enzyme, mitochondrial superoxide dismutase. SOD2 polymorphisms are of interest because of their potential roles in the modulation of free radical-mediated macromolecular damage during aging.
We identified a new splice variant of SOD2 in human lymphoblastoid cell lines (LCLs). The alternatively spliced product was originally detected by exon trapping of a minigene in order to examine the consequences of an intronic polymorphism found upstream of exon 4 (nucleotide 8136, 10T vs 9T). Examination of the transcripts derived from the endogenous loci in five LCLs with or without the intron 3 polymorphism revealed low levels of an in-frame deletion of exon 4 that were different from those detected by the exon trap assay. This suggested that exon trapping of the minigene unmasked the effect of the 10T vs 9T polymorphism on the splicing of the adjacent exon.
We also determined the frequencies of single nucleotide polymorphisms in a sample of US African-Americans and non-African-Americans ages 65 years and older who participated in the 1999 wave of the National Long Term Care Survey (NLTCS). Particularly striking differences between African-Americans and non-African-Americans were found for the frequencies of genotypes at the 10T/9T intron 3 polymorphism.
Exon trapping can unmask in vitro splicing differences caused by a 10T/9T intron 3 polymorphism. Given the recent evidence that SOD2 is in a region on chromosome 6 linked to susceptibility to hypertension, it will be of interest to investigate possible associations of this polymorphism with cardiovascular disorders.
- Reverse Transcriptase Polymerase Chain Reaction
- Cystic Fibrosis Transmembrane Conductance Regulator
- Lymphoblastoid Cell Line
- Diabetic Polyneuropathy
- SOD2 mRNA
The free radical theory of aging is currently the leading hypothesis for a fundamental mechanism of aging. It potentially unifies aging mechanisms across an array of species [1–5]. The theory postulates that reactive oxygen species (ROS) such as superoxide and hydroxyl radicals, mainly generated as toxic by-products of normal cellular metabolism during aerobic respiration in mitochondria, give rise to damage to wide variety of macromolecules, including DNA. Animals possess enzymes that scavenge ROS. Among the most important of such antioxidant defenses are the superoxide dismutases (SOD), which facilitate the dismutation of superoxide free radical (O2 -.) . There are three types of SODs in higher eukaryotes, a copper/zinc containing cytoplasmic enzyme, SOD1; a manganese containing mitochondrial enzyme, SOD2; and a copper/zinc containing SOD3 that is prevalent in extracellular spaces [6, 7].
SOD2 is constitutively expressed in most cells . SOD2 gene dosage correlates with longevity in fruit flies and in mice . Mice that are heterozygous for an SOD2 knockout accumulate excessive damage to mitochondrial DNA as they age . In these mouse lines, there were no compensatory increases of other antioxidant enzymes such as SOD1, catalase or glutathione peroxidase . There is little understanding of the role of SODs in the modulation of longevity or health status in aging humans. Human disorders due to the SOD2 mutations have yet to be identified. However, a number of studies describe associations of SOD2 polymorphisms (particularly Ala16Val) with age-related disorders. One such study of elderly people with an Ashkenazi ethnic background suggested a possible association between Val at codon 16 of SOD2 and longevity . Individuals who were homozygous for the valine allele at position 16 (also described as the -9 position of the mitochondrial targeting sequence) were shown to have an increased risk of non-small cell lung cancer [13, 14]. The Val allele may also be a predisposing factor for diabetic polyneuropathy (DPN) in Russians with type 1 diabetes mellitus . Finally, homozygosity for the Val allele was also associated with non-familial idiopathic dilated cardiomyopathy among Japanese . This was thought to be due to a decreased efficiency of SOD2 transport into mitochondria .
In contrast to the above reports, the alanine allele has been associated with high risk of breast cancer among U.S. and Finnish Caucasians [17, 18] and with an increased risk of developing psoriatic arthritis in the Taiwanese population .
As is typical of such association studies, most have not been confirmed in independent populations. Linkage studies are therefore of interest. A particularly cogent example is a report of a genome scan of microsatelite markers among hypertensive African-Americans; linkage was described at 6q24 and 21q21 , with Z scores, respectively, of 4.14 and 4.34. Markers associated with hypertension in African-Americans are between 112 cM and 166 cM of 6q. SOD2 is located at 160 cM.
The National Long-Term Care Survey (NLTCS) reports on the functional status of cross-sectional and longitudinal population-based samples from U.S. Medicare recipients ≥ age 65 . During the 1999 wave of that survey, biological samples were obtained. This provided an opportunity to examine SOD2 polymorphisms among different ethnic groups.
The 1999 NLTCS populations
DNA was obtained from a total of 644 blood samples. Each DNA isolation gave an average of 40 ug of DNA, sufficient for 800 PCR reactions. We also received 2,102 buccal wash samples of which 2,078 provided DNA suitable for banking. There was sufficient DNA from buccal wash samples to permit the performance of up to ~400 PCR reactions.
Screening studies were performed using genomic DNA from peripheral blood. Because a striking difference in the frequencies for the polymorphism at the nucleotide position 8,136 was observed among African-Americans (see below), studies of that polymorphism were extended to the larger sample of genomic DNA from buccal samples.
SOD2 primer sequences and PCR conditions used for this study
Primers (Top-forward primer, Middle-reverse primer, Bottom-sequencing primer)
Size of amplicon
P-GTCCACTGTCGCCATTG (for exon 1)
Q-GTGGTACGCTGACTGAC (for exon 2)
Frequencies of investigated SOD2 polymorphisms in the 1996 NLTCS blood samples (N = 613–639).
Is the intron 3 T9/T10 polymorphism associated with splicing variations?
Dramatically altered 8136 allele frequencies within the African-American population
The 644 blood samples were derived from 596 Caucasians, 43 African-Americans, 2 Asian/Pacific islanders, and 3 "others". Statistically significant differences between African-Americans and Caucasians were observed for three out of five polymorphisms. At the Ala16Val polymorphism in exon 2 (which has been reported to alter the risk of several disorders), the Val (GTT) was more frequent than Ala (GCT) in the African-American population (P = 0.044). The C8039G polymorphism in intron 3 also showed a statistically significant difference between African-Americans and Caucasians (P = 0.031). A statistically significant difference was not observed for the case of the adjacent G8116T polymorphism (P = 0.093).
Among the five polymorphisms studied, the most striking finding was a difference in the frequencies of the polymorphism at nucleotide 8,136 (Table 2). A single base pair deletion was designated as T9 and the wild type designated as T10. Among 40 African-Americans, 18% carried one T9 allele, significantly higher than the average for the total population (1.7%) and strikingly higher when compared to the population of non-African-Americans (0.7%). We then examined this T9/T10 polymorphism in buccal samples that were obtained from African-Americans. Out of 67 African-Americans who provided buccal wash samples, 54 (81%) were T10/T10 and 13 (10%) were T10/T9. No individual has been found to be homozygous for the T9 allele. The allele frequencies from blood and buccal samples were comparable. We concluded that the T9 allele is far more frequent in elderly African-Americans than in non-African-Americans (p < 0.0001).
The T9/T10 SOD2 polymorphism was newly identified in our laboratory in the course of screening of the 1999 NLTCS population for SOD2 polymorphisms. Our analysis indicated that the T9 allele is highly prevalent in elderly African-Americans as compared to elderly Caucasian-Americans (p < 0.0001). The location of T9/T10 at an intron/exon boundary region suggested a role in pre-mRNA processing. An exon trap analysis was quite consistent with this interpretation. There are a number of examples in the human genetics literature in which intron polymorphisms appear to have altered the ratio of wild type and alternative splicing patterns [22, 23]. In some cases, altered splicing was shown be responsible for the pathogenetic mechanism leading to disease. Interesting examples are the sequence variations at the 3' end of intron 8 of cystic fibrosis transmembrane conductance regulator (CFTR). These are thought to affect splicing of exon 9, resulting in the reduction of functional CFTR . For the case of our SOD2 T9/T10 polymorphism, however, it was not possible to confirm the splicing bias using RNA from intact cells (lymphoblastoid cell lines). Similarly, we were unable to determine what fractions of the transcript lacking exon 4 were derived from T10 and T9 alleles in the heterozygote.
Previous work from our laboratory showed that intronic polymorphisms generated cryptic splicing products by exon trapping or by assays utilizing a mini gene and RT-PCR, but there was little evidence that this occurred in vivo . The most likely interpretation of this discordance in the results of test tube-based versus cell-based methods is that the former exclude important cis-acting domains residing in more distal domains of DNA. In addition, trans-acting factors may differ among different cell types and species. These elements are likely to be capable of modulating splicing.
Although we have so far failed to provide molecular evidence for a functional distinction between the polymorphic allele of interest, given the striking difference in prevalence among African-Americans as compared to Caucasian-Americans, it will be important to carry out appropriate epidemiological studies of the impact of the two alleles (and, eventually, contrasting haplotypes), upon health and survival in the two ethnic groups. The NLTCS should provide valuable material for such an investigation, as there is linkage to Medicare data and data on functional status and survival. A definitive study, however, must await the availability of biological material from future cohorts, as the 1999 wave of the study could not provide biological specimens from many institutionalized subjects.
Perhaps the most immediate question would be the potential impact of this allelic variation of SOD2 upon susceptibility to hypertension among African-Americans. First, hypertension is pervasive in this population. Secondly, it has deleterious effects upon multiple organs. Third, while susceptibility to hypertension is clearly modulated by multiple genetic loci, a region of chromosome 6 that encompasses the SOD2 locus is among those that have been implicated . Fourth, African-Americans provide an excellent opportunity to carry out admixture mapping (reviewed by [27, 28]. It has been estimated that ~75% of the alleles of African-Americans are of African ancestry whilst ~25% are of European ancestry . Gene scans on the parental and admixed populations were successfully employed to identify hypertension loci by Zhu et al. (2005). The most robust finding was a region on the long arm of chromosome 6 that included the SOD2 locus. While the region contains many genes, potentially including more than one susceptibility genes, our observations suggest that one important candidate gene is SOD2.
We demonstrated that the 10T/9T polymorphism in the SOD2 intron 3 clearly results in alternative splicing of exon 4 by an in vitro exon trap assay using minigenes. This difference was masked when endogenous genes were expressed within cells. Over-representation of 9T alleles in African-Americans nevertheless raises the possibility that this intron 3 polymorphism or polymorphisms in linkage disequilibrium may be involved in the altered susceptibilities to a range of disorders that are more common in African-American populations.
Sample collection and DNA isolation
Both peripheral blood and buccal wash samples were obtained from the 1999 NLTCS participants. Informed consent was obtained. Peripheral blood mononuclear cells were isolated by Ficoll-Hypaque gradients (Amersham, CA, USA) gradient. Buccal wash samples were centrifuged at 3,000 × g prior to DNA isolation. Genomic DNA from blood and buccal wash samples were isolated with QIAamp blood isolation kits or QIAamp tissue isolation kits, respectively (Qiagen Inc. Valencia, CA), following protocols specified by the manufacturer. For high throughput DNA isolation, a Qiagen BioRobot 3000 system was employed for DNA manipulation. Isolated genomic DNA samples were kept in 96 well microtiter plates in duplicates, and stored at -80°C.
Single nucleotide polymorphism (SNP) detection
The SNPs in SOD2 were determined through direct sequencing of SOD2 exons and adjacent portions of introns. Briefly, genomic DNA was first amplified by 4 different primer pairs (Table 1), designed using Primer Express software (Applied Biosystems Co., Foster City, CA), to obtain DNA fragments containing exons 1 and 2, exon 3, exon4, and exon 5, respectively. The primer pairs were designed to specifically amplify the designated exon(s) with maximum adjacent upstream and downstream intron sequences. All PCR products were first examined on agarose DNA electrophoresis to verify size and were treated with shrimp alkaline phosphatase and exonuclease I to eliminate excess nucleotides. The subsequent sequencing of the PCR products was performed by a Big Dye-terminator cycle sequencing kit (ABI) on an ABI PRISM 310 automated DNA sequencer with factory-supplied parameters. The sequencing primers for the 4 PCR products are also listed in Table 1. Newly found polymorphisms were confirmed by sequencing in the reverse orientation in all of the samples.
Distributions of allele frequencies within the groups were examined with the chi-square test (one degree of freedom) to test for Hardy-Weinberg equilibrium. Comparisons of the distributions between two groups were tested by the modified Fisher exact test (two degrees of freedom) [30, 31]. The latter was also examined with the standard Pearson chi-square test .
Exon trapping was performed to reveal the potential of SNPs to promote differential splicing of SOD2 exons. Genomic fragments of different alleles of SOD2 were first amplified with a set of primers (forward: gct ggt ccc att atc taa tag c; reverse: cag tag agc atc tct ccc aaa tg) encompassing SOD2 intron 2 to exon 5 by long-range PCR (Roche Diagnostics Co. Indianapolis, IL). The DNA fragments carrying the SNPs in question were cloned by TA cloning techniques (Invitrogen, Hercules, CA), and the SNPs were verified by sequencing. The selected genomic DNA fragments with the intended SNPs were subcloned into a pSPL3 splicing vector for exon trapping analysis (Invitrogen Co. Carlsbad, CA). The nucleotide sequences at the polymorphic site were confirmed by sequencing, and the constructs were verified via the digestion patterns of combinations of restriction enzymes. The exon trapping procedures were performed according to the manufacturer's specifications. The genomic DNA in pSPL3 was first isolated with a plasmid purification kit (Qiagen), and was used to transfect COS-7 cells for transient expression of the reporter system. Exon trapping reporter transcripts were examined by extracting total RNA from COS-7 cells, followed by the reverse transcriptase polymerase chain reaction (RT-PCR) with primers specified in the kit for amplification of the reporter gene. The amplified reporter gene fragments were visualized on 2% agarose gels. Variable lengths of the amplified fragments revealed differential splicing as affected by the inserted alleles in the reporter system.
Establishment of lymphoblastoid cell lines and RT-PCR
Cryopreserved lymphocytes were thawed and immortalized with Epstein-Barr virus using a published protocol . Total RNA was isolated using Trizol (Gibco BRL) and was reverse-transcribed with random hexamers with GeneAmp RNA PCR kit (Perkin Elmer Cetus Inc., Wellesley, MA) following the manufacturer's instructions. PCR reactions were performed as described previously [30, 31]. The products were separated on a 2% agarose-1xTAE gel.
This work was supported by NIH grants U01AG07198, P20AG12852 and P30AG13280. We thank Ms. Elice Y. Kim for editorial assistance.
- Harman D: Aging: A theory based on free radical and radiation chemistry. J Gerontol. 1956, 11: 298-300.View ArticlePubMedGoogle Scholar
- Harman D: Free radical theory of aging. Mutat Res. 1992, 275 (3-6): 257-266.View ArticlePubMedGoogle Scholar
- Finkel T, Holbrook NJ: Oxidants, oxidative stress and the biology of ageing. Nature. 2000, 408 (6809): 239-247. 10.1038/35041687.View ArticlePubMedGoogle Scholar
- Hekimi S, Guarente L: Genetics and the specificity of the aging process. Science. 2003, 299 (5611): 1351-1354. 10.1126/science.1082358.View ArticlePubMedGoogle Scholar
- Schriner SE, Linford NJ, Martin GM, Treuting P, Ogburn CE, Emond M, Coskun PE, Ladiges W, Wolf N, Van Remmen H, Wallace DC, Rabinovitch PS: Extension of murine life span by overexpression of catalase targeted to mitochondria. Science. 2005, 308 (5730): 1909-1911. 10.1126/science.1106653.View ArticlePubMedGoogle Scholar
- Fridovich I: Superoxide radical and superoxide dismutases. Annu Rev Biochem. 1995, 64: 97-112. 10.1146/annurev.bi.64.070195.000525.View ArticlePubMedGoogle Scholar
- Landis GN, Tower J: Superoxide dismutase evolution and life span regulation. Mech Ageing Dev. 2005, 126 (3): 365-379. 10.1016/j.mad.2004.08.012.View ArticlePubMedGoogle Scholar
- Lai CC, Huang WH, Askari A, Wang Y, Sarvazyan N, Klevay LM, Chiu TH: Differential regulation of superoxide dismutase in copper-deficient rat organs. Free Radic Biol Med. 1994, 16 (5): 613-620. 10.1016/0891-5849(94)90061-2.View ArticlePubMedGoogle Scholar
- Tower J: Transgenic methods for increasing Drosophila life span. Mech Ageing Dev. 2000, 118 (1-2): 1-14. 10.1016/S0047-6374(00)00152-4.View ArticlePubMedGoogle Scholar
- Williams MD, Van Remmen H, Conrad CC, Huang TT, Epstein CJ, Richardson A: Increased oxidative damage is correlated to altered mitochondrial function in heterozygous manganese superoxide dismutase knockout mice. J Biol Chem. 1998, 273 (43): 28510-28515. 10.1074/jbc.273.43.28510.View ArticlePubMedGoogle Scholar
- Van Remmen H, Salvador C, Yang H, Huang TT, Epstein CJ, Richardson A: Characterization of the antioxidant status of the heterozygous manganese superoxide dismutase knockout mouse. Arch Biochem Biophys. 1999, 363 (1): 91-97. 10.1006/abbi.1998.1060.View ArticlePubMedGoogle Scholar
- Stessman J, Maaravi Y, Hammerman-Rozenberg R, Cohen A, Nemanov L, Gritsenko I, Gruberman N, Ebstein RP: Candidate genes associated with ageing and life expectancy in the Jerusalem longitudinal study. Mech Ageing Dev. 2005, 126 (2): 333-339. 10.1016/j.mad.2004.08.025.View ArticlePubMedGoogle Scholar
- Wang LI, Miller DP, Sai Y, Liu G, Su L, Wain JC, Lynch TJ, Christiani DC: Manganese superoxide dismutase alanine-to-valine polymorphism at codon 16 and lung cancer risk. J Natl Cancer Inst. 2001, 93 (23): 1818-1821. 10.1093/jnci/93.23.1818.View ArticlePubMedGoogle Scholar
- Liu G, Zhou W, Park S, Wang LI, Miller DP, Wain JC, Lynch TJ, Su L, Christiani DC: The SOD2 Val/Val genotype enhances the risk of nonsmall cell lung carcinoma by p53 and XRCC1 polymorphisms. Cancer. 2004, 101 (12): 2802-2808. 10.1002/cncr.20716.View ArticlePubMedGoogle Scholar
- Zotova EV, Chistiakov DA, Savost'ianov KV, Bursa TR, Galeev IV, Strokov IA, Nosikov VV: [Association of the SOD2 Ala(-9)Val and SOD3 Arg213Gly polymorphisms with diabetic polyneuropathy in patients with diabetes mellitus type 1]. Mol Biol (Mosk). 2003, 37 (3): 404-408. 10.1023/A:1024287327107.View ArticleGoogle Scholar
- Hiroi S, Harada H, Nishi H, Satoh M, Nagai R, Kimura A: Polymorphisms in the SOD2 and HLA-DRB1 genes are associated with nonfamilial idiopathic dilated cardiomyopathy in Japanese. Biochem Biophys Res Commun. 1999, 261 (2): 332-339. 10.1006/bbrc.1999.1036.View ArticlePubMedGoogle Scholar
- Mitrunen K, Sillanpaa P, Kataja V, Eskelinen M, Kosma VM, Benhamou S, Uusitupa M, Hirvonen A: Association between manganese superoxide dismutase (MnSOD) gene polymorphism and breast cancer risk. Carcinogenesis. 2001, 22 (5): 827-829. 10.1093/carcin/22.5.827.View ArticlePubMedGoogle Scholar
- Ambrosone CB, Freudenheim JL, Thompson PA, Bowman E, Vena JE, Marshall JR, Graham S, Laughlin R, Nemoto T, Shields PG: Manganese superoxide dismutase (MnSOD) genetic polymorphisms, dietary antioxidants, and risk of breast cancer. Cancer Res. 1999, 59 (3): 602-606.PubMedGoogle Scholar
- Yen JH, Tsai WC, Lin CH, Ou TT, Hu CJ, Liu HW: Manganese superoxide dismutase gene polymorphisms in psoriatic arthritis. Dis Markers. 2003, 19 (6): 263-265.View ArticlePubMedGoogle Scholar
- Zhu X, Luke A, Cooper RS, Quertermous T, Hanis C, Mosley T, Gu CC, Tang H, Rao DC, Risch N, Weder A: Admixture mapping for hypertension loci with genome-scan markers. Nat Genet. 2005, 37 (2): 177-181. 10.1038/ng1510.View ArticlePubMedGoogle Scholar
- Manton KG, Corder L, Stallard E: Chronic disability trends in elderly United States populations: 1982-1994. Proc Natl Acad Sci U S A. 1997, 94 (6): 2593-2598. 10.1073/pnas.94.6.2593.View ArticlePubMedPubMed CentralGoogle Scholar
- Pagani F, Baralle FE: Genomic variants in exons and introns: identifying the splicing spoilers. Nat Rev Genet. 2004, 5 (5): 389-396. 10.1038/nrg1327.View ArticlePubMedGoogle Scholar
- Chu CS, Trapnell BC, Curristin S, Cutting GR, Crystal RG: Genetic basis of variable exon 9 skipping in cystic fibrosis transmembrane conductance regulator mRNA. Nat Genet. 1993, 3 (2): 151-156. 10.1038/ng0293-151.View ArticlePubMedGoogle Scholar
- Hu Q, Cool BH, Wang B, Hearn MG, Martin GM: A candidate molecular mechanism for the association of an intronic polymorphism of FE65 with resistance to very late onset dementia of the Alzheimer type. Hum Mol Genet. 2002, 11 (4): 465-475. 10.1093/hmg/11.4.465.View ArticlePubMedGoogle Scholar
- Krushkal J, Ferrell R, Mockrin SC, Turner ST, Sing CF, Boerwinkle E: Genome-wide linkage analyses of systolic blood pressure using highly discordant siblings. Circulation. 1999, 99 (11): 1407-1410.View ArticlePubMedGoogle Scholar
- McKeigue PM: Prospects for admixture mapping of complex traits. Am J Hum Genet. 2005, 76 (1): 1-7. 10.1086/426949.View ArticlePubMedGoogle Scholar
- Darvasi A, Shifman S: The beauty of admixture. Nat Genet. 2005, 37 (2): 118-119. 10.1038/ng0205-118.View ArticlePubMedGoogle Scholar
- Destro-Bisol G, Maviglia R, Caglia A, Boschi I, Spedini G, Pascali V, Clark A, Tishkoff S: Estimating European admixture in African Americans by using microsatellites and a microsatellite haplotype (CD4/Alu). Hum Genet. 1999, 104 (2): 149-157. 10.1007/s004390050928.View ArticlePubMedGoogle Scholar
- Castro E, Edland SD, Lee L, Ogburn CE, Deeb SS, Brown G, Panduro A, Riestra R, Tilvis R, Louhija J, Penttinen R, Erkkola R, Wang L, Martin GM, Oshima J: Polymorphisms at the Werner locus: II. 1074Leu/Phe, 1367Cys/Arg, longevity, and atherosclerosis. Am J Med Genet. 2000, 95 (4): 374-380. 10.1002/1096-8628(20001211)95:4<374::AID-AJMG14>3.0.CO;2-4.View ArticlePubMedGoogle Scholar
- Castro E, Ogburn CE, Hunt KE, Tilvis R, Louhija J, Penttinen R, Erkkola R, Panduro A, Riestra R, Piussan C, Deeb SS, Wang L, Edland SD, Martin GM, Oshima J: Polymorphisms at the Werner locus: I. Newly identified polymorphisms, ethnic variability of 1367Cys/Arg, and its stability in a population of Finnish centenarians. Am J Med Genet. 1999, 82 (5): 399-403. 10.1002/(SICI)1096-8628(19990219)82:5<399::AID-AJMG8>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Conover WJ: Practical Nonparametric Statistics. 1980, John Wiley and Sons, pp 162-169. 2nd editionGoogle Scholar
- Beck JC, Beiswanger CM, John EM, Satariano E, West D: Successful transformation of cryopreserved lymphocytes: a resource for epidemiological studies. Cancer Epidemiol Biomarkers Prev. 2001, 10 (5): 551-554.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/8/7/prepub
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.