Oil palm suspension cultures were initiated by transferring the gel-like friable embryogenic tissue onto liquid medium supplemented with auxins. In this study, transcripts that were differentially expressed in oil palm suspension cells cultured at different auxin concentrations were examined using suppression subtractive hybridization. Total RNA was first isolated from oil palm suspension cells proliferated in liquid medium with different hormone concentrations for 6 months. Four different hormone combinations: T1 (0.1 mg/l 2,4-D and 1.0 mg/l NAA), T2 (0.4 mg/l 2,4-D and 1.0 mg/l NAA), T3 (1.0 mg/l NAA), and T4 (0.4 mg/l 2,4-D) were used for the treatments. The first and second subtractions were performed using samples T1 and T2 in forward and reverse order. The other two subtractions were forward and reverse subtractions of T3 and T4, respectively. Reverse northern analyses showed that 14.13% of these clones were preferentially expressed in T1, 13.70% in T2, 14.75% in T3, and 15.70% in T4. Among the 294 cDNA clones that were sequenced, 61 contigs (assembled from 165 sequences) and 129 singletons were obtained. Among the 61 contigs, 10 contigs consist of sequences from treatment T1, 8 contigs were from treatment T2, 10 contigs were contains sequences of treatment T3 and 13 contigs contains sequences of treatment T4. Northern analyses of five transcripts that were shown to be differentially expressed in the oil palm suspension cells by reverse northern analysis revealed that transcripts 16A1 (a putative lignostilbene-alpha,beta-dioxygenase, EgLSD) and 16H12 (a putative ethylene responsive 6, EgER6) were differentially expressed in oil palm suspension cells treated with different levels of auxin.
We have characterized an oil palm (Elaeis guineensis Jacq.) constitutive promoter that is derived from a translationally control tumor protein (TCTP) gene. The TCTP promoter was fused transcriptionally with the gusA reporter gene and transferred to monocot and dicot systems in order to study its regulatory role in a transient expression study. It was found that the 5' region of TCTP was capable of driving the gusA expression in all the oil palm tissues tested, including immature embryo, embryogenic callus, embryoid, young leaflet from mature palm, green leaf, mesocarp and stem. It could also be used in dicot systems as it was also capable of driving gusA expression in tobacco leaves. The results indicate that the TCTP promoter could be used for the production of recombinant proteins that require constitutive expression in the plant system.
This study was directed at the understanding of the function of CCoAOMT isolated from Acacia auriculiformis x Acacia mangium. Full length cDNA of the Acacia hybrid CCoAOMT (AhCCoAOMT) was 1024-bp long, containing 750-bp coding regions, with one major open reading frame of 249 amino acids. On the other hand, full length genomic sequence of the CCoAOMT (AhgflCCoAOMT) was 2548 bp long, containing three introns and four exons with a 5' untranslated region (5'UTR) of 391 bp in length. The 5'UTR of the characterized CCoAOMT gene contains various regulatory elements. Southern analysis revealed that the Acacia hybrid has more than three copies of the CCoAOMT gene. Real-time PCR showed that this gene was expressed in root, inner bark, leaf, flower and seed pod of the Acacia hybrid. Downregulation of the homologous CCoAOMT gene in tobacco by antisense (AS) and intron-containing hairpin (IHP) constructs containing partial AhCCoAOMT led to reduction in lignin content. Expression of the CCoAOMT in AS line (pART-HAS78-03) and IHP line (pART-HIHP78-06) was reduced respectively by 37 and 75% compared to the control, resulting in a decrease in the estimated lignin content by 24 and 56%, respectively. AhCCoAOMT was found to have altered not only S and G units but also total lignin content, which is of economic value to the pulp industry. Subsequent polymorphism analysis of this gene across eight different genetic backgrounds each of A. mangium and A. auriculiformis revealed 47 single nucleotide polymorphisms (SNPs) in A. auriculiformis CCoAOMT and 30 SNPs in A. mangium CCoAOMT.
Backcrossing together with simple sequence repeat marker strategy was adopted to improve popular Malaysian chilli Kulai (Capsicum annuum L.) for heat tolerance. The use of molecular markers in backcross breeding and selection contributes significantly to overcoming the main drawbacks such as increase linkage drag and time consumption, in the ancient manual breeding approach (conventional), and speeds up the genome recovery of the recurrent parent. The strategy was adopted to introgress heat shock protein gene(s) from AVPP0702 (C. annuum L.), which are heat-tolerant, into the genetic profile of Kulai, a popular high-yielding chilli but which is heat sensitive. The parents were grown on seed trays, and parental screening was carried out with 252 simple sequence repeat markers. The selected parents were crossed and backcrossed to generate F1 hybrids and backcross generations. Sixty-eight markers appeared to be polymorphic and were used to assess the backcross generation; BC1F1, BC2F1 and BC3F1. The average recipient allele of the selected four BC1F1 plants was 80.75% which were used to produce the BC2F1 generation. BC1-P7 was the best BC1F1 plant because it had the highest recovery at 83.40% and was positive to Hsp-linked markers (Hsp70-u2 and AGi42). After three successive generations of backcrossing, the average genome recovery of the recurrent parent in the selected plants in BC3F1 was 95.37%. Hsp gene expression analysis was carried out on BC1F1, BC2F1 and BC3F1 selected lines. The Hsp genes were found to be up-regulated when exposed to heat treatment. The pattern of Hsp expression in the backcross generations was similar to that of the donor parent. This confirms the successful introgression of a stress-responsive gene (Hsp) into a Kulai chilli pepper variety. Furthermore, the yield performance viz. plant height, number of fruits, fruit length and weight and total yield of the improved plant were similar with the recurrent parent except that the plant height was significantly lower than the Kulai (recurrent) parent.
We applied an integrative approach using multiple methods to verify cytosine methylation in the chloroplast DNA of the multicellular brown alga Saccharina japonica. Cytosine DNA methylation is a heritable process which plays important roles in regulating development throughout the life cycle of an organism. Although methylation of nuclear DNA has been studied extensively, little is known about the state and role of DNA methylation in chloroplast genomes, especially in marine algae. Here, we have applied an integrated approach encompassing whole-genome bisulfite sequencing, methylated DNA immunoprecipitation, gene co-expression networks and photophysiological analyses to provide evidence for the role of chloroplast DNA methylation in a marine alga, the multicellular brown alga Saccharina japonica. Although the overall methylation level was relatively low in the chloroplast genome of S. japonica, gametophytes exhibited higher methylation levels than sporophytes. Gene-specific bisulfite-cloning sequencing provided additional evidence for the methylation of key photosynthetic genes. Many of them were highly expressed in sporophytes whereas genes involved in transcription, translation and biosynthesis were strongly expressed in gametophytes. Nucleus-encoded photosynthesis genes were co-expressed with their chloroplast-encoded counterparts potentially contributing to the higher photosynthetic performance in sporophytes compared to gametophytes where these co-expression networks were less pronounced. A nucleus-encoded DNA methyltransferase of the DNMT2 family is assumed to be responsible for the methylation of the chloroplast genome because it is predicted to possess a plastid transit peptide.
A high leaf vein density is both an essential feature of C4 photosynthesis and a foundation trait to C4 evolution, ensuring the optimal proportion and proximity of mesophyll and bundle sheath cells for permitting the rapid exchange of photosynthates. Two rice mutant populations, a deletion mutant library with a cv. IR64 background (12,470 lines) and a T-DNA insertion mutant library with a cv. Tainung 67 background (10,830 lines), were screened for increases in vein density. A high throughput method with handheld microscopes was developed and its accuracy was supported by more rigorous microscopy analysis. Eight lines with significantly increased leaf vein densities were identified to be used as genetic stock for the global C4 Rice Consortium. The candidate population was shown to include both shared and independent mutations and so more than one gene controlled the high vein density phenotype. The high vein density trait was found to be linked to a narrow leaf width trait but the linkage was incomplete. The more genetically robust narrow leaf width trait was proposed to be used as a reliable phenotypic marker for finding high vein density variants in rice in future screens.
Drought tolerance is an important quantitative trait with multipart phenotypes that are often further complicated by plant phenology. Different types of environmental stresses, such as high irradiance, high temperatures, nutrient deficiencies, and toxicities, may challenge crops simultaneously; therefore, breeding for drought tolerance is very complicated. Interdisciplinary researchers have been attempting to dissect and comprehend the mechanisms of plant tolerance to drought stress using various methods; however, the limited success of molecular breeding and physiological approaches suggests that we rethink our strategies. Recent genetic techniques and genomics tools coupled with advances in breeding methodologies and precise phenotyping will likely reveal candidate genes and metabolic pathways underlying drought tolerance in crops. The WRKY transcription factors are involved in different biological processes in plant development. This zinc (Zn) finger protein family, particularly members that respond to and mediate stress responses, is exclusively found in plants. A total of 89 WRKY genes in japonica and 97 WRKY genes in O. nivara (OnWRKY) have been identified and mapped onto individual chromosomes. To increase the drought tolerance of rice (Oryza sativa L.), research programs should address the problem using a multidisciplinary strategy, including the interaction of plant phenology and multiple stresses, and the combination of drought tolerance traits with different genetic and genomics approaches, such as microarrays, quantitative trait loci (QTLs), WRKY gene family members with roles in drought tolerance, and transgenic crops. This review discusses the newest advances in plant physiology for the exact phenotyping of plant responses to drought to update methods of analysing drought tolerance in rice. Finally, based on the physiological/morphological and molecular mechanisms found in resistant parent lines, a strategy is suggested to select a particular environment and adapt suitable germplasm to that environment.
A cross between IR64 (high-yielding but drought-susceptible) and Aday Sel (drought-tolerant) rice cultivars yielded a stable line with enhanced grain yield under drought screening field trials at International Rice Research Institute. The major effect qDTY4.1 drought tolerance and yield QTL was detected in the IR77298-14-1-2-10 Backcrossed Inbred Line (BIL) and its IR87705-7-15-B Near Isogenic Line (NIL) with 93.9% genetic similarity to IR64. Although rice yield is extremely susceptible to water stress at reproductive stage, currently, there is only one report on the detection of drought-responsive microRNAs in inflorescence tissue of a Japonica rice line. In this study, more drought-responsive microRNAs were identified in the inflorescence tissues of IR64, IR77298-14-1-2-10 and IR87705-7-15-B via next-generation sequencing. Among the 32 families of inflorescence-specific non-conserved microRNAs that were identified, 22 families were up-regulated in IR87705-7-15-B. Overall 9 conserved and 34 non-conserved microRNA families were found as drought-responsive in rice inflorescence with 5 conserved and 30 non-conserved families induced in the IR87705-7-15-B. The observation of more drought-responsive non-conserved microRNAs may imply their prominence over conserved microRNAs in drought response mechanisms of rice inflorescence. Gene Ontology annotation analysis on the target genes of drought-responsive microRNAs identified in IR87705-7-15-B revealed over-representation of biological processes including development, signalling and response to stimulus. Particularly, four inflorescence-specific microRNAs viz. osa-miR5485, osa-miR5487, osa-miR5492 and osa-miR5517, and two non-inflorescence specific microRNAs viz. osa-miR169d and osa-miR169f.2 target genes that are involved in flower or embryonic development. Among them, osa-miR169d, osa-miR5492 and osa-miR5517 are related to flowering time control. It is also worth mentioning that osa-miR2118 and osa-miR2275, which are implicated in the biosynthesis of rice inflorescence-specific small interfering RNAs, were induced in IR87705-7-15-B but repressed in IR77298-14-1-2-10. Further, gene search within qDTY4.1 QTL region had identified multiple copies of NBS-LRR resistance genes (potential target of osa-miR2118), subtilisins and genes implicated in stomatal movement, ABA metabolism and cuticular wax biosynthesis.
Developing drought-tolerant rice varieties with higher yield under water stressed conditions provides a viable solution to serious yield-reduction impact of drought. Understanding the molecular regulation of this polygenic trait is crucial for the eventual success of rice molecular breeding programmes. microRNAs have received tremendous attention recently due to its importance in negative regulation. In plants, apart from regulating developmental and physiological processes, microRNAs have also been associated with different biotic and abiotic stresses. Hence here we chose to analyze the differential expression profiles of microRNAs in three drought treated rice varieties: Vandana (drought-tolerant), Aday Sel (drought-tolerant) and IR64 (drought-susceptible) in greenhouse conditions via high-throughput sequencing.
Calamus palustris Griff. is an economically important dioecious rattan species in Southeast Asia. However, dioecy and onset of flowering at 3-4 years old render uncertainties in desired female:male seedling ratios to establish a productive seed orchard for this rattan species. We constructed a subtractive library for male floral tissue to understand the genetic mechanism for gender determination in C. palustris. The subtractive library produced 1536 clones with 1419 clones of high quality. Reverse Northern screening showed 313 clones with differential expression, and sequence analyses clustered them into 205 unigenes, including 32 contigs and 173 singletons. The subtractive library was further validated with reverse transcription-quantitative polymerase chain reaction analysis. Homology identification classified the unigenes into 12 putative functional proteins with 83% unigenes showing significant match to proteins in databases. Functional annotations of these unigenes revealed genes involved in male flower development, including MADS-box genes, pollen-related genes, phytohormones for flower development, and male flower organ development. Our results showed that the male floral genes may play a vital role in sex determination in C. palustris. The identified genes can be exploited to understand the molecular basis of sex determination in C. palustris.
Casuarina equisetifolia L. is an important multi-purpose, fast growing and widely planted tree species native to tropical and subtropical coastlines of Australia, Southeast Asia, Malaysia, Melanesia, Polynesia and New Caledonia. It is a nitrogen-fixing tree mainly used for charcoal making, construction poles, landscaping, timber, pulp, firewood, windbreaks, shelterbelts, soil erosion and sand dune stabilization. Casuarina wood is presently used for paper and pulp production. Raw material with reduced lignin is highly preferred to increase the pulp yield. Hence, understanding the molecular regulation of wood formation in this tree species is vital for selecting industrially suitable phenotypes for breeding programs. The lignin biosynthetic pathway has been extensively studied in tree species like Eucalypts, poplars, pines, Picea, Betula and Acacia sp. However, studies on wood formation at molecular level is presently lacking in casuarinas. Hence, in the present study, the transcriptome of the developing secondary tissues of 15 years old Casuarina equiseitfolia subsp. equisetifolia was sequenced, de novo assembled, annotated and mapped to functional pathways. Transcriptome sequencing generated a total of 26,985 transcripts mapped to 31 pathways. Mining of the annotated data identified nine genes involved in lignin biosynthesis pathway and relative expression of the transcripts in four tissues including scale-like leaves, needle-like brachlets, wood and root were documented. The expression of CeCCR1 and CeF5H were found to be significantly high in wood tissues, while maximum expression of CeHCT was documented in stem. Additionally, CeTUBA and CeH2A were identified as the most stable reference transcript for normalization of qRT-PCR data in C. equisetifolia. The present study is the first wood genomic resource in C. equisetifolia, which will be valuable for functional genomics research in this genus.
Hevea brasiliensis remains the primary crop commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. Here, we described the transcriptional events related to jasmonic acid (JA)- and linolenic acid (LA)-induced secondary laticifer differentiation (SLD) in H. brasiliensis clone RRIM 600 based on RNA-seq approach. Histochemical approach proved that JA- and LA-treated samples resulted in SLD in H. brasiliensis when compared to ethephon and untreated control. RNA-seq data resulted in 86,614 unigenes, of which 2,664 genes were differentially expressed in JA and LA-induced secondary laticifer harvested from H. brasiliensis bark samples. Among these, 450 genes were unique to JA and LA as they were not differentially expressed in ethephon-treated samples compared with the untreated samples. Most transcription factors from the JA- and LA-specific dataset were classified under MYB, APETALA2/ethylene response factor (AP2/ERF), and basic-helix-loop-helix (bHLH) gene families that were involved in tissue developmental pathways, and we proposed that Bel5-GA2 oxidase 1-KNOTTED-like homeobox complex are likely involved in JA- and LA-induced SLD in H. brasiliensis. We also discovered alternative spliced transcripts, putative novel transcripts, and cis-natural antisense transcript pairs related to SLD event. This study has advanced understanding on the transcriptional regulatory network of SLD in H. brasiliensis.
Chlorella is a popular microalga with robust physiological and biochemical characteristics, which can be cultured under various conditions. The exploration of the small RNA content of Chlorella could improve strategies for the enhancement of metabolite production from this microalga. In this study, stress was introduced to the Chlorella sorokiniana culture to produce high-value metabolites such as carotenoids and phenolic content. The small RNA transcriptome of C. sorokiniana was sequenced, focusing on microRNA (miRNA) content. From the analysis, 98 miRNAs were identified in cultures subjected to normal and stress conditions. The functional analysis result showed that the miRNA targets found were most often involved in the biosynthesis of secondary metabolites, followed by protein metabolism, cell cycle, and porphyrin and chlorophyll metabolism. Furthermore, the biosynthesis of secondary metabolites such as carotenoids, terpenoids, and lipids was found mostly in stress conditions. These results may help to improve our understanding of regulatory mechanisms of miRNA in the biological and metabolic process of Chlorella species. It is important and timely to determine the true potential of this microalga species and to support the potential for genetic engineering of microalgae as they receive increasing focus for their development as an alternative source of biofuel, food, and health supplements.
Oil palm is one of the most productive oil-producing crops and can store up to 90% oil in its fruit mesocarp. Oil palm fruit is a sessile drupe consisting of a fleshy mesocarp from which palm oil is extracted. Biochemical changes in the mesocarp cell walls, polyamines, and hormones at different ripening stages of oil palm fruits were studied, and the relationship between the structural and the biochemical metabolism of oil palm fruits during ripening is discussed. Time-course analysis of the changes in expression of polyamines, hormones, and cell-wall-related genes and metabolites provided insights into the complex processes and interactions involved in fruit development. Overall, a strong reduction in auxin-responsive gene expression was observed from 18 to 22 weeks after pollination. High polyamine concentrations coincided with fruit enlargement during lipid accumulation and latter stages of maturation. The trend of abscisic acid (ABA) concentration was concordant with GA₄ but opposite to the GA₃ profile such that as ABA levels increase the resulting elevated ABA/GA₃ ratio clearly coincides with maturation. Polygalacturonase, expansin, and actin gene expressions were also observed to increase during fruit maturation. The identification of the master regulators of these coordinated processes may allow screening for oil palm variants with altered ripening profiles.
Black pepper (Piper nigrum L.) is one of the most popular and oldest spices in the world with culinary uses and various pharmacological properties. In order to satisfy the growing worldwide demand for black pepper, improved productivity of pepper is highly desirable. A primary constraint in black pepper production is the non-synchronous nature of flower development and non-uniform fruit ripening within a spike. The uneven ripening of pepper berries results in a high labour requirement for selective harvesting contributes to low productivity and affects the quality of the pepper products. In Malaysia, there are a few recommended varieties for black pepper planting, each having some limitations in addition to the useful characteristics. Therefore, a comparative study of different black pepper varieties will provide a better understanding of the mechanisms regulates fruit development and ripening. Plant hormones are known to influence the fruit development process and their roles in black pepper flower and fruit development were inferred based on the probe-based gene expression analysis and the quantification of the multiple plant hormones using high-performance liquid chromatography coupled to tandem mass spectrometry (HPLC-MS/MS). In this study, jasmonic acid and salicylic acid were found to play roles in flowering and fruit setting, whereas auxin, gibberellin and cytokinins are important for fruit growth. Abscisic acid has positive role in fruit maturation and ripening in the development process. Distinct pattern of plant hormones related gene expression profiles with the hormones accumulation profiles suggested a complex network of regulation is involved in the signaling process and crosstalk between plant hormones was another layer of regulation in the black pepper fruit development mechanisms. The current study provides clues to help in elucidating the timing of the action of each specific plant hormone during fruit development and ripening which could be applied to enhance our ability to control the ripening process, leading to improving procedures for the production and post-harvest handling of pepper fruits.
Plant cells contain a diverse repertoire of RNA-binding proteins (RBPs) that coordinate a network of post-transcriptional regulation. RBPs govern diverse developmental processes by modulating the gene expression of specific transcripts. Recent gene annotation and RNA sequencing clearly showed that heterogeneous nuclear ribonucleoprotein (hnRNP)-like proteins which form a family of RBPs, are also expressed in higher plants and serve specific plant functions. In addition to their involvement in post-transcriptional regulation from mRNA capping to translation, they are also involved in telomere regulation, gene silencing and regulation in chloroplast. Here, we review the involvement of plant hnRNP-like proteins in post-transcription regulation of RNA processes and their functional roles in control of plant developmental processes especially plant-specific functions including flowering, chloroplastic-specific mRNA regulation, long-distance phloem transportation and plant responses to environmental stresses.
The spontaneous occurrence of resistance to the herbicide glyphosate in weed species has been an extremely infrequent event, despite over 20 years of extensive use. Recently, a glyphosate-resistant biotype of goosegrass (Eleusine indica) was identified in Malaysia exhibiting an LD(50) value approximately 2- to 4-fold greater than the sensitive biotype collected from the same region. A comparison of the inhibition of 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) activity by glyphosate in extracts prepared from the resistant (R) and sensitive (S) biotypes revealed an approximately 5-fold higher IC(50)(glyphosate) for the (R) biotype. Sequence comparisons of the predicted EPSPS mature protein coding regions from both biotypes revealed four single-nucleotide differences, two of which result in amino acid changes. One of these changes, a proline to serine substitution at position 106 in the (R) biotype, corresponds to a substitution previously identified in a glyphosate-insensitive EPSPS enzyme from Salmonella typhimurium. Kinetic data generated for the recombinant enzymes suggests that the second substitution identified in the (R) EPSPS does not contribute significantly to its reduced glyphosate sensitivity. Escherichia coli aroA- (EPSPS deficient) strains expressing the mature EPSPS enzyme from the (R) biotype exhibited an approximately 3-fold increase in glyphosate tolerance relative to strains expressing the mature EPSPS from the (S) biotype. These results provide the first evidence for an altered EPSPS enzyme as an underlying component of evolved glyphosate resistance in any plant species.
Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca(2+) transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca(2+) transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.
Glyphosate-resistant populations of Eleusine indica are widespread in several states of Malaysia. A whole-plant bioassay confirmed that eight out of the 17 populations tested were resistant to glyphosate at double the recommended rate of 2.44 kg ha-1. Screening with allele-specific PCR (AS-PCR) revealed that resistant plants contained an EPSPS gene with either the homozygous S/S-106 or the heterozygous P/S-106 alleles. All susceptible plants contained only the homozygous P/P-106 allele. In addition, DNA sequences of the full-length EPSPS gene from one susceptible (SB) and four resistant (R2, R6, R8 and R11) populations revealed an amino acid substitution of T102I in all the resistant plants, while another substitution of P381L was only found in resistant populations R6 and R11. The significance of the P381L mutation was examined by Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) and residue interaction network (RIN) analyses, which suggests the P381L mutation may contribute to resistance. Mutations at 102 and 106 occur widely in the EPSPS gene of glyphosate-resistant E. indica populations from Malaysia with the TIPS mutation. In addition, the P381L mutation could also contribute to resistance.
Transgenic technologies have been applied to a wide range of biological research. However, information on the potential epigenetic effects of transgenic technology is still lacking. Here, we show that the transgenic process can simultaneously induce both genetic and epigenetic changes in rice. We analyzed genetic, epigenetic, and phenotypic changes in plants subjected to tissue culture regeneration, using transgenic lines expressing the same coding sequence from two different promoters in transgenic lines of two rice cultivars: Wuyunjing7 (WYJ7) and Nipponbare (NP). We determined the expression of OsNAR2.1 in two overexpression lines generated from the two cultivars, and in the RNA interference (RNAi) OsNAR2.1 line in NP. DNA methylation analyses were performed on wild-type cultivars (WYJ7 and NP), regenerated lines (CK, T0 plants), segregation-derived wild-type from pOsNAR2.1-OsNAR2.1 (SDWT), pOsNAR2.1-OsNAR2.1, pUbi-OsNAR2.1, and RNAi lines. Interestingly, we observed global methylation decreased in the T0 regenerated line of WYJ7 (CK-WJY7) and pOsNAR2.1-OsNAR2.1 lines but increased in pUbi-OsNAR2.1 and RNAi lines of NP. Furthermore, the methylation pattern in SDWT returned to the WYJ7 level after four generations. Phenotypic changes were detected in all the generated lines except for SDWT. Global methylation was found to decrease by 13% in pOsNAR2.1-OsNAR2.1 with an increase in plant height of 4.69% compared with WYJ7, and increased by 18% in pUbi-OsNAR2.1 with an increase of 17.36% in plant height compared with NP. This suggests an absence of a necessary link between global methylation and the phenotype of transgenic plants with OsNAR2.1 gene over-expression. However, epigenetic changes can influence phenotype during tissue culture, as seen in the massive methylation in CK-WYJ7, T0 regenerated lines, resulting in decreased plant height compared with the wild-type, in the absence of a transformed gene. We conclude that in the transgenic lines the phenotype is mainly determined by the nature and function of the transgene after four generations of transformation, while the global epigenetic modification is dependent on the genetic background. Our research suggests an innovative insight in explaining the reason behind the occurrence of transgenic plants with random and undesirable phenotypes.