The association of single nucleotide polymorphisms with endometrial cancer risk and prognosis

Endometrial cancer is one of the most common female diseases in developed nations and is the most commonly diagnosed gynaecological cancer in Australia. The disease is commonly classified by histology: endometrioid or non-endometrioid endometrial cancer. While non-endometrioid endometrial cancers are accepted to be high-grade, aggressive cancers, endometrioid cancers (comprising 80% of all endometrial cancers diagnosed) generally carry a favourable patient prognosis. However, endometrioid endometrial cancer patients endure significant morbidity due to surgery and radiotherapy used for disease treatment, and patients with recurrent disease have a 5-year survival rate of less than 50%. Genetic analysis of women with endometrial cancer could uncover novel markers associated with disease risk and/or prognosis, which could then be used to identify women at high risk and for the use of specialised treatments. Proteases are widely accepted to play an important role in the development and progression of cancer. This PhD project hypothesised that SNPs from two protease gene families, the matrix metalloproteases (MMPs, including their tissue inhibitors, TIMPs) and the tissue kallikrein-related peptidases (KLKs) would be associated with endometrial cancer susceptibility and/or prognosis. In the first part of this study, optimisation of the genotyping techniques was performed. Results from previously published endometrial cancer genetic association studies were attempted to be validated in a large, multicentre replication set (maximum cases n = 2,888, controls n = 4,483, 3 studies). The rs11224561 progesterone receptor SNP (PGR, A/G) was observed to be associated with increased endometrial cancer risk (per A allele OR 1.31, 95% CI 1.12-1.53; p-trend = 0.001), a result which was initially reported among a Chinese sample set. Previously reported associations for the remaining 8 SNPs investigated for this section of the PhD study were not confirmed, thereby reinforcing the importance of validation of genetic association studies. To examine the effect of SNPs from the MMP and KLK families on endometrial cancer risk, we selected the most significantly associated MMP and KLK SNPs from genome-wide association study analysis (GWAS) to be genotyped in the GWAS replication set (cases n = 4,725, controls n = 9,803, 13 studies). The significance of the MMP24 rs932562 SNP was unchanged after incorporation of the stage 2 samples (Stage 1 per allele OR 1.18, p = 0.002; Combined Stage 1 and 2 OR 1.09, p = 0.002). The rs10426 SNP, located 3' to KLK10 was predicted by bioinformatic analysis to effect miRNA binding. This SNP was observed in the GWAS stage 1 result to exhibit a recessive effect on endometrial cancer risk, a result which was not validated in the stage 2 sample set (Stage 1 OR 1.44, p = 0.007; Combined Stage 1 and 2 OR 1.14, p = 0.08). Investigation of the regions imputed surrounding the MMP, TIMP and KLK genes did not reveal any significant targets for further analysis. Analysis of the case data from the endometrial cancer GWAS to identify genetic variation associated with cancer grade did not reveal SNPs from the MMP, TIMP or KLK genes to be statistically significant. However, the representation of SNPs from the MMP, TIMP and KLK families by the GWAS genotyping platform used in this PhD project was examined and observed to be very low, with the genetic variation of four genes (MMP23A, MMP23B, MMP28 and TIMP1) not captured at all by this technique. This suggests that comprehensive candidate gene association studies will be required to assess the role of SNPs from these genes with endometrial cancer risk and prognosis. Meta-analysis of gene expression microarray datasets curated as part of this PhD study identified a number of MMP, TIMP and KLK genes to display differential expression by endometrial cancer status (MMP2, MMP10, MMP11, MMP13, MMP19, MMP25 and KLK1) and histology (MMP2, MMP11, MMP12, MMP26, MMP28, TIMP2, TIMP3, KLK6, KLK7, KLK11 and KLK12). In light of these findings these genes should be prioritised for future targeted genetic association studies. Two SNPs located 43.5 Mb apart on chromosome 15 were observed from the GWAS analysis to be associated with increased endometrial cancer grade, results that were validated in silico in two independent datasets. One of these SNPs, rs8035725 is located in the 5' untranslated region of a MYC promoter binding protein DENND4A (Stage 1 OR 1.15, p = 9.85 x 10P -5 P, combined Stage 1 and in silico validation OR 1.13, p = 5.24 x 10P -6 P). This SNP has previously been reported to alter the expression of PTPLAD1, a gene involved in the synthesis of very long fatty acid chains and in the Rac1 signaling pathway. Meta-analysis of gene expression microarray data found PTPLAD1 to display increased expression in the aggressive non-endometrioid histology compared with endometrioid endometrial cancer, suggesting that the causal SNP underlying the observed genetic association may influence expression of this gene. Neither rs8035725 nor significant SNPs identified by imputation were predicted bioinformatically to affect transcription factor binding sites, indicating that further studies are required to assess their potential effect on other regulatory elements. The other grade- associated SNP, rs6606792, is located upstream of an inferred pseudogene, ELMO2P1 (Stage 1 OR 1.12, p = 5 x 10P -5 P; combined Stage 1 and in silico validation OR 1.09, p = 3.56 x 10P -5 P). Imputation of the ±1 Mb region surrounding this SNP revealed a cluster of significantly associated variants which are predicted to abolish various transcription factor binding sites, and would be expected to decrease gene expression. ELMO2P1 was not included on the microarray platforms collected for this PhD, and so its expression could not be investigated. However, the high sequence homology of ELMO2P1 with ELMO2, a gene important to cell motility, indicates that ELMO2 could be the parent gene for ELMO2P1 and as such, ELMO2P1 could function to regulate the expression of ELMO2. Increased expression of ELMO2 was seen to be associated with increasing endometrial cancer grade, as well as with aggressive endometrial cancer histological subtypes by microarray meta-analysis. Thus, it is hypothesised that SNPs in linkage disequilibrium with rs6606792 decrease the transcription of ELMO2P1, reducing the regulatory effect of ELMO2P1 on ELMO2 expression. Consequently, ELMO2 expression is increased, cell motility is enhanced leading to an aggressive endometrial cancer phenotype. In summary, these findings have identified several areas of research for further study. The results presented in this thesis provide evidence that a SNP in PGR is associated with risk of developing endometrial cancer. This PhD study also reports two independent loci on chromosome 15 to be associated with increased endometrial cancer grade, and furthermore, genes associated with these SNPs to be differentially expressed according in aggressive subtypes and/or by grade. The studies reported in this thesis support the need for comprehensive SNP association studies on prioritised MMP, TIMP and KLK genes in large sample sets. Until these studies are performed, the role of MMP, TIMP and KLK genetic variation remains unclear. Overall, this PhD study has contributed to the understanding of genetic variation involvement in endometrial cancer susceptibility and prognosis. Importantly, the genetic regions highlighted in this study could lead to the identification of novel gene targets to better understand the biology of endometrial cancer and also aid in the development of therapeutics directed at treating this disease.

Saliva-derived DNA performs well in large-scale, high-density single-nucleotide polymorphism microarray studies

As of June 2009, 361 genome-wide association studies (GWAS) had been referenced by the HuGE database. GWAS require DNA from many thousands of individuals, relying on suitable DNA collections. We recently performed a multiple sclerosis (MS) GWAS where a substantial component of the cases (24%) had DNA derived from saliva. Genotyping was done on the Illumina genotyping platform using the Infinium Hap370CNV DUO microarray. Additionally, we genotyped 10 individuals in duplicate using both saliva- and blood-derived DNA. The performance of blood- versus saliva-derived DNA was compared using genotyping call rate, which reflects both the quantity and quality of genotyping per sample and the “GCScore,” an Illumina genotyping quality score, which is a measure of DNA quality. We also compared genotype calls and GCScores for the 10 sample pairs. Call rates were assessed for each sample individually. For the GWAS samples, we compared data according to source of DNA and center of origin. We observed high concordance in genotyping quality and quantity between the paired samples and minimal loss of quality and quantity of DNA in the saliva samples in the large GWAS sample, with the blood samples showing greater variation between centers of origin. This large data set highlights the usefulness of saliva DNA for genotyping, especially in high-density single-nucleotide polymorphism microarray studies such as GWAS.

A rapid and sensitive UPLC-ESI MS method for analysis of isofraxidin, a natural antistress compound, and its metabolites in rat plasma

Isofraxidin is one of the main bioactive constituents in the root of Acanthopanax senticosus, which has antifatigue, antistress, and immuno-accommondating effects. In this study, an ultraperformance LC (UPLC)-ESI MS method was developed for analyzing isofraxidin and its metabolites in rat plasma. The analysis was performed on a UPLC coupled with ESI MS (quadropole MS tandem TOF MS). The lower LOD (LLOD) for isofraxidin was 0.25 ng/mL, the intraday precision was less than 10%, the interday precision was less than 10%, and the extraction recovery was more than 80%. Isofraxidin and two metabolites (M1 and M2) were detected in rat plasma after oral administration of isofraxidin, and the molecular polarities of M1 and M2 were both increased compared to isofraxidin. The metabolites were identified as 5,6-dihydroxyl-7-methoxycoumarin and 5-hydroxyl-6,7-dimethoxycoumarin when subjected to parent ion spectra, product ion spectra, and extract mass and element composition analyses.

The complete nucleotide sequence of subterranean clover mottle virus

The complete nucleotide sequence of Subterranean clover mottle virus (SCMoV) genomic RNA has been determined. The SCMoV genome is 4,258 nucleotides in length. It shares most nucleotide and amino acid sequence identity with the genome of Lucerne transient streak virus (LTSV). SCMoV RNA encodes four overlapping open reading frames and has a genome organisation similar to that of Cocksfoot mottle virus (CfMV). ORF1 and ORF4 are predicted to encode single proteins. ORF2 is predicted to encode two proteins that are derived from a -1 translational frameshift between two overlapping reading frames (ORF2a and ORF2b). A search of amino acid databases did not find a significant match for ORF1 and the function of this protein remains unclear. ORF2a contains a motif typical of chymotrypsin-like serine proteases and ORF2b has motifs characteristically present in positive-stranded RNA-dependent RNA polymerases. ORF4 is likely to be expressed from a subgenomic RNA and encodes the viral coat protein. The ORF2a/ORF2b overlapping gene expression strategy used by SCMoV and CfMV is similar to that of the poleroviruses and differ from that of other published sobemoviruses. These results suggest that the sobemoviruses could now be divided into two distinct subgroups based on those that express the RNA-dependent RNA polymerase from a single, in-frame polyprotein, and those that express it via a -1 translational frameshifting mechanism.

Nucleotide sequence of coat protein gene for GPV isolate of barley yellow dwarf virus and construction of expression plasmid for plant

GPV is a Chinese serotype isolate of barley yellow dwarf virus (BYDV) that has no reaction with antiserum of MAV, PAV, SGV, RPV and RMV The sequence of the coat protein (CP) of GPV isolate of BYDV was identified and its amino acid sequence was deduced. The coding region for the putative GPV CP is 603 bases nucleotides and encodes a Mr 22 218 (22 ku) protein. The same as MAV, PAV and RPV, GPV contained a second ORF within the coat protein coding region. This protein of 17 024 Mr (17 ku) is thought to correspond to the Virion protein genome linked (Vpg). Sequence comparisons of the CP coding region between the GPV isolate of BYDV and other isolates of BYDV have been done. The nucleotide and amino acid sequence homology of GPV has a greater identity to the sequence of RPV than those of PAV and MAV. The GPV CP sequence stored 83.7% of nucleotide similarity and 77.5% of deduced amino acid similarity, whereas that of the PAV and MAV shared 56.9%, 53.2% and 44.1%, 43.8% respectively. According to BYDV-GPV CP sequence, two primers were designed. The cDNA of CP was produced by RT-PCR. Full-length cDNA of CP was inserted into plasmid to construct expression plasmids named pPPI1, pPPI2 and pPPI5 based on different promoters. The recombinant plasmids were identified by using α-32P-dATP labelled CP probe, α-32P-ATP labelled GPV RNA probe and sequencing to confirm real GPV CP gene cDNA in plasmids.

Nucleotide sequences of an Australian and a Canadian isolate of potato leafroll luteovirus and their relationships with two European isolates

The genomes of an Australian and a Canadian isolate of potato leafroll virus have been cloned and sequenced. The sequences of both isolates are similar (about 93%), but the Canadian isolate (PLRV-C) is more closely related (about 98% identity) to a Scottish (PLRV-S) and a Dutch isolate (PLRV-N) than to the Australian isolate (PLRV-A). The 5'-terminal 18 nucleotide residues of PLRV-C, PLRV-A, PLRV-N and beet western yellows virus have 17 residues in common. In contrast, PLRV-S shows no obvious similarity in this region. PLRV-A and PLRV-C genomic sequences have localized regions of marked diversity, in particular a 600 nucleotide residue sequence in the polymerase gene. These data provide a world-wide perspective on the molecular biology of PLRV strains and their comparison with other luteoviruses and related RNA plant viruses suggests that there are two major subgroups in the plant luteoviruses.

Single nucleotide polymorphisms and mRNA expression of VEGF-A in papillary thyroid carcinoma: Potential markers for aggressive phenotypes

BACKGROUND AND OBJECTIVES Polymorphisms of the VEGF gene are known to affect the biological behaviour of cancers but have seldom been studied in thyroid cancer. The aim of the current study is to evaluate the prevalence and relevance of VEGF-A polymorphisms and mRNA expression in papillary thyroid carcinoma (PTC). MATERIALS AND METHODS Genomic DNA and total RNA were isolated from paraffin-embedded tissue from 91 PTC (51 conventional PTC and 40 follicular variant) and 78 control thyroid tissues. Three DNA polymorphisms (+936C > T, +405C > G and -141A > C) in the 3' and 5' untranslated region (3'-UTR, 5'-UTR) of VEGF-A were studied using PCR and RFLP. Also, the mRNA expression of VEGF-A in these tissues was studied by real-time PCR. RESULTS Distribution of polymorphisms in the 5'-UTR (VEGF-VEGF -141A > C and +405C > G) and 3'-UTR (VEGF +936C > T) were all significantly different in PTC and benign thyroid tissue (p = 0.0001, 0.001 and 0.028 respectively). The VEGF -141 C allele was more common in PTC with lymph node metastases (p = 0.026). VEGF + 405 Galleles andVEGF +936 CC genotype were more common in PTC of advanced pathological staging (p = 0.018 and 0.017 respectively). Also, increased expression of VEGF-A mRNA was noted in PTC compared to control (p = 0.009). Within the group of patients with conventional PTC, those with lymph nodal metastases had a higher level of VEGF-A mRNA expression than other patients (p = 0.0003). CONCLUSION These findings suggest that VEGF polymorphisms and mRNA expression may predict the aggressiveness behaviour of thyroid cancer.

Determination of glutathione transferase (GSTT1-1) activities in different tissues based on formation of radioactive metabolites using 35S-glutathione

A new system has been developed to determine enzyme activities of glutathione transferase θ (GSTT1-1) based on radiometric product detection resulting from the enzymic reaction of methyl chloride with 35S-labelled glutathione. In principle, the method is universally applicable for determination of glutathione transferase activities towards a multiplicity of substrates. The method distinguishes between erythrocyte GSTT1-1 activities of human 'non-conjugators', 'low conjugators' and 'high conjugators'. Application to cytosol preparations of livers and kidneys of male and female Fischer 344 and B6C3F1 mice reveals differential GSTT1-1 activities in hepatic and renal tissues. These ought to be considered in species-specific modellings of organ toxicities of chlorinated hydrocarbons.

Association study of MTHFD1 coding polymorphisms R134K and R653Q with migraine susceptibility

Objective There is evidence that folate metabolism has a role in migraine pathophysiology, particularly in the migraine with aura subtype. In this study we investigate whether two non-synonymous single nucleotide polymorphisms (SNPs), rs1950902 (C401T; R134K) and rs2236225 (G1958A; R653Q), in MTHDF1 are associated with migraine in an Australian case-control population. Background Increased plasma levels of homocysteine (HCy), one of the metabolites produced in the folate pathway, has been found to be a risk factor for migraine. There is also a genetic link, as a common polymorphism (C667T) that reduces the catalytic activity of MTHFR, the enzyme that catalyses the formation of HCy, is associated with an increase in risk of the migraine with aura (MA) subtype. MTHFD1 is a crucial multifunctional enzyme that catalyses three separate reactions of the folate pathway and therefore variants in MTHFD1 may also influence migraine susceptibility. Methods The R134K and R653Q variants in MTHFD1 were genotyped in an Australian cohort of 520 unrelated migraineurs (162 were diagnosed with migraine without aura [MO] and 358 with MA) and 520 matched controls. Data were analysed for association with migraine and for interaction with the MTHFR C667T polymorphism. Results We find no significant differences in genotype or allele frequencies for either SNP between migraineurs and controls, or when either MO or MA cases were compared to controls. In addition these MTHFD1 polymorphisms did not appear to influence the risk of MA conferred by the MTHFR 667T allele. Conclusions We find no evidence for association of the MTHFD1 R134K and R653Q polymorphisms with migraine in our Australian case-control population. However, as folate metabolism appears to be important in migraine, particularly with respect to the aura component, future studies using high throughput methods to expand the number of SNPs in folate-related genes genotyped and investigation of interactions between SNPs may be justified.

Association between single-nucleotide polymorphisms in growth factor genes and quality of life in men with prostate cancer and the general population

Purpose Improved survival for men with prostate cancer has led to increased attention to factors influencing quality of life (QOL). As protein levels of vascular endothelial growth factor (VEGF) and insulin-like growth factor 1 (IGF-1) have been reported to be associated with QOL in people with cancer, we sought to identify whether single-nucleotide polymorphisms (SNPs) of these genes were associated with QOL in men with prostate cancer. Methods Multiple linear regression of two data sets (including approximately 750 men newly diagnosed with prostate cancer and 550 men from the general population) was used to investigate SNPs of VEGF and IGF-1 (10 SNPs in total) for associations with QOL (measured by the SF-36v2 health survey). Results Men with prostate cancer who carried the minor ‘T’ allele for IGF-1 SNP rs35767 had higher mean Role-Physical scale scores (≥0.3 SD) compared to non-carriers (p < 0.05). While this association was not identified in men from the general population, one IGF-1 SNP rs7965399 was associated with higher mean Bodily Pain scale scores in men from the general population that was not found in men with prostate cancer. Men from the general population who carried the rare ‘C’ allele had higher mean Bodily Pain scale scores (≥0.3 SD) than non-carriers (p < 0.05). Conclusions Through identifying SNPs that are associated with QOL in men with prostate cancer and men from the general population, this study adds to the mapping of complex interrelationships that influence QOL and suggests a role for IGF-I in physical QOL outcomes. Future research may identify biomarkers associated with increased risk of poor QOL that could assist in the provision of pre-emptive support for those identified at risk.

Transferability of cupped oyster EST (Expressed Sequence Tag)-derived SNP (Single Nucleotide Polymorphism) markers to related crassostrea and ostrea species

Single nucleotide polymorphisms (SNPs) are widely acknowledged as the marker of choice for many genetic and genomic applications because they show co-dominant inheritance, are highly abundant across genomes and are suitable for high-throughput genotyping. Here we evaluated the applicability of SNP markers developed from Crassostrea gigas and C. virginica expressed sequence tags (ESTs) in closely related Crassostrea and Ostrea species. A total of 213 putative interspecific level SNPs were identified from re-sequencing data in six amplicons, yielding on average of one interspecific level SNP per seven bp. High polymorphism levels were observed and the high success rate of transferability show that genic EST-derived SNP markers provide an efficient method for rapid marker development and SNP discovery in closely related oyster species. The six EST-SNP markers identified here will provide useful molecular tools for addressing questions in molecular ecology and evolution studies including for stock analysis (pedigree monitoring) in related oyster taxa.

High-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population

Objective Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. Methods A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. Results We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly149Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly149Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). Conclusion This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

Single nucleotide polymorphisms in clinics: Fantasy or reality for cancer?

Single nucleotide polymorphisms (SNPs) have been classically used for dissecting various human complex disorders using candidate gene studies. During the last decade, large scale SNP analysis i.e. genome-wide association studies (GWAS) have provided an agnostic approach to identify possible genetic loci associated with heterogeneous disease such as cancer susceptibility, prognosis of survival or drug response. Further, the advent of new technologies, including microarray based genotyping as well as high throughput next generation sequencing has opened new avenues for SNPs to be used in clinical practice. It is speculated that the utility of SNPs to understand the mechanisms, biology of variable drug response and ultimately treatment individualization based on the individual’s genome composition will be indispensable in the near future. In the current review, we discuss the advantages and disadvantages of the clinical utility of genetic variants in disease risk-prediction, prognosis, clinical outcome and pharmacogenomics. The lessons and challenges for the utility of SNP based biomarkers are also discussed, including the need for additional functional validation studies.

Variations in the levels of aminoacylation and modified nucleotide content between total tRNAs from chloroplasts and etioplasts in cucumber cotyledons

Total tRNAs isolated from chloroplasts and etioplasts of cucumber cotyledons were compared with respect toamino acid acceptance, isoacceptor distribution and extent of modification. Aminoacylation of the tRNAs with nine different amino acids studied indicated that the relative acceptor activities of chloroplast total tRNAs for four amino acids are significantly higher than etioplast total tRNAs. Two dimensional polyacrylamide gel electrophoresis(2D-PAGE) of chloroplast total tRNAs separated at least 32 spots, while approximately 41 spots were resolved from etioplast total tRNAs. Comparison of the reversed-phase chromatography (RPC-5) profiles of chloroplast and etioplast leucyl-, lysyl-, phenylalanyl-, and valyl-tRNA species showed no qualitative differences in the elution profiles. However, leucyl-, lysyl- and valyl-tRNA species showed quantitative differences in the relative amounts of the isoaccepting species present in chloroplasts and etioplasts. The analysis of modified nucleotides of total tRNAs from the two plastid types indicated that total tRNA from etioplasts was undermodified with respect to ribothymidine, isopentenyladenosine/hydroxy-isopentenyladenosine, 1 -methylguanosine and 2-o-methylguanosine. This indicates that illumination may cause de novo synthesis of chloroplast tRNAmodifying enzymes encoded for by nuclear genes leading to the formation of highly modified tRNAs in chloroplasts. Based on these results, we speculate that the observed decrease in levels of aminoacylation, variations in the relative amounts of certain isoacceptors, and differences in the electrophoretic mobilities of some extra tRNA spots in the etioplast total tRNAs as compared to chloroplast total tRNAs could be due to some partially undermodified etioplast tRNAs. Taken together, the data suggested that the light-induced transformation of etioplasts into chloroplasts is accompanied by increases in the relative levels of some functional chloroplast tRNAs by post transcriptional nucleotide modifications.