Genes In Cotton Research Articles

Identification, Characterization, and Expression Patterns of TCP Genes and microRNA319 in Cotton.

The TEOSINTE BRANCHED 1, CYCLOIDEA, and PROLIFERATING CELL FACTORS (TCP) gene family is a group of plant-specific transcription factors that have versatile functions in developmental processes and stress responses. In this study, a total of 73 TCP genes in upland cotton were identified and characterizated. Phylogenetic analysis classified them into three subgroups: 50 belonged to PCF, 16 to CIN, and 7 to CYC/TB1. GhTCP genes are randomly distributed in 22 of the 26 chromosomes in cotton. Expression patterns of GhTCPs were analyzed in 10 tissues, including different developmental stages of ovule and fiber, as well as under heat, salt, and drought stresses. Transcriptome analysis showed that 44 GhTCP genes exhibited varied transcript accumulation patterns in the tested tissues and 41 GhTCP genes were differentially expressed in response to heat, salt, and drought stresses. Furthermore, three GhTCP genes of the CIN clade were found to contain miR319-binding sites. An anti-correlation expression of GhTCP21 and GhTCP54 was analyzed with miR319 under salt and drought stress. Our results lay the foundation for understanding the complex mechanisms of GhTCP-mediated developmental processes and abiotic stress-signaling transduction pathways in cotton.

Natural variation in a CENTRORADIALIS homolog contributed to cluster fruiting and early maturity in cotton

BackgroundPlant architecture and the vegetative-reproductive transition have major impacts on the agronomic success of crop plants, but genetic mechanisms underlying these traits in cotton (Gossypium spp.) have not been identified.ResultsWe identify four natural mutations in GoCEN-Dt associated with cluster fruiting (cl) and early maturity. The situ hybridization shows that GhCEN is preferentially expressed in cotton shoot apical meristems (SAM) of the main stem and axillary buds. Constitutive GhCEN-Dt overexpression suppresses the transition of the cotton vegetative apex to a reproductive shoot. Silencing GoCEN leads to early flowering and determinate growth, and in tetraploids causes the main stem to terminate in a floral bud, a novel phenotype that exemplifies co-adaptation of polyploid subgenomes and suggests new research and/or crop improvement approaches. Natural cl variations are enriched in cottons adapted to high latitudes with short frost-free periods, indicating that mutants of GoCEN have been strongly selected for early maturity.ConclusionWe show that the cotton gene GoCEN-Dt, a homolog of Antirrhinum CENTRORADIALIS, is responsible for determinate growth habit and cluster fruiting. Insight into the genetic control of branch and flower differentiation offers new approaches to develop early maturing cultivars of cotton and other crops with plant architecture appropriate for mechanical harvesting.

LncRNAs in polyploid cotton interspecific hybrids are derived from transposon neofunctionalization

BackgroundInterspecific hybridization and whole genome duplication are driving forces of genomic and organism diversification. But the effect of interspecific hybridization and whole genome duplication on the non-coding portion of the genome in particular remains largely unknown. In this study, we examine the profile of long non-coding RNAs (lncRNAs), comparing them with that of coding genes in allotetraploid cotton (Gossypium hirsutum), its putative diploid ancestors (G. arboreum; G. raimondii), and an F1 hybrid (G. arboreum × G. raimondii, AD).ResultsWe find that most lncRNAs (80%) that were allelic expressed in the allotetraploid genome. Moreover, the genome shock of hybridization reprograms the non-coding transcriptome in the F1 hybrid. Interestingly, the activated lncRNAs are predominantly transcribed from demethylated TE regions, especially from long interspersed nuclear elements (LINEs). The DNA methylation dynamics in the interspecies hybridization are predominantly associated with the drastic expression variation of lncRNAs. Similar trends of lncRNA bursting are also observed in the progress of polyploidization. Additionally, we find that a representative novel lncRNA XLOC_409583 activated after polyploidization from a LINE in the A subgenome of allotetraploid cotton was involved in control of cotton seedling height.ConclusionOur results reveal that the processes of hybridization and polyploidization enable the neofunctionalization of lncRNA transcripts, acting as important sources of increased plasticity for plants.

Recursive Paleohexaploidization Shaped the Durian Genome.

The durian (Durio zibethinus) genome has recently become available, and analysis of this genome reveals two paleopolyploidization events previously inferred as shared with cotton (Gossypium spp.). Here, we reanalyzed the durian genome in comparison with other well-characterized genomes. We found that durian and cotton were actually affected by different polyploidization events: hexaploidization in durian ∼19-21 million years ago (mya) and decaploidization in cotton ∼13-14 mya. Previous interpretations of shared polyploidization events may have resulted from the elevated evolutionary rates in cotton genes due to the decaploidization and insufficient consideration of the complexity of plant genomes. The decaploidization elevated evolutionary rates of cotton genes by ∼64% compared to durian and explained a previous ∼4-fold over dating of the event. In contrast, the hexaploidization in durian did not prominently elevate gene evolutionary rates, likely due to its long generation time. Moreover, divergent evolutionary rates probably explain 98.4% of reconstructed phylogenetic trees of homologous genes being incongruent with expected topology. The findings provide further insight into the roles played by polypoidization in the evolution of genomes and genes, and they suggest revisiting existing reconstructed phylogenetic trees.

Genome-wide analysis of Rf-PPR-like (RFL) genes and a new InDel marker development for Rf1 gene in cytoplasmic male sterile CMS-D2 Upland cotton

BackgroundCytoplasmic male sterility in flowering plants is a convenient way to use heterosis via hybrid breeding and may be restored by nuclear restorer-of-fertility (Rf) genes. In most cases, Rf genes encoded pentatricopeptide repeat (PPR) proteins and several Rf genes are present in clusters of similar Rf-PPR-like (RFL) genes. However, the Rf genes in cotton were not fully characterized until now.ResultsIn total, 35 RFL genes were identified in G. hirsutum, 16 in G. arboreum, and 24 in G. raimondii. Additionally, four RFL-rich regions were identified; the RFL-rich region in Gh_D05 is the probable location of Rf-PPR genes in cotton and will be studied further in the future. Furthermore, an insertion sequence was identified in the promoter sequence of Gh_D05G3392 gene in the restorer line, as compared with the CMS-D2 line and maintainer lines. An InDel-R marker was then developed and could be used to distinguish the restorer line carrying Rf1 from other genotypes without the Rf1 allele.ConclusionIn this study, genome-wide identification and analysis of RFL genes have identified the candidate Rf-PPR genes for CMS in Gossypium. The identification and analysis of RFL genes and sequence variation analysis will be useful for cloning Rf genes in the future and also for three-line hybrid breeding in cotton.

The comparison of four mitochondrial genomes reveals cytoplasmic male sterility candidate genes in cotton

BackgroundThe mitochondrial genomes of higher plants vary remarkably in size, structure and sequence content, as demonstrated by the accumulation and activity of repetitive DNA sequences. Incompatibility between mitochondrial genome and nuclear genome leads to non-functional male reproductive organs and results in cytoplasmic male sterility (CMS). CMS has been used to produce F1 hybrid seeds in a variety of plant species.ResultsHere we compared the mitochondrial genomes (mitogenomes) of Gossypium hirsutum sterile male lines CMS-2074A and CMS-2074S, as well as their restorer and maintainer lines. First, we noticed the mitogenome organization and sequences were conserved in these lines. Second, we discovered the mitogenomes of 2074A and 2074S underwent large-scale substitutions and rearrangements. Actually, there were five and six unique chimeric open reading frames (ORFs) in 2074A and 2074S, respectively, which were derived from the recombination between unique repetitive sequences and nearby functional genes. Third, we found out four chimeric ORFs that were differentially transcribed in sterile line (2074A) and fertile-restored line.ConclusionsThese four novel and recombinant ORFs are potential candidates that confer CMS character in 2074A. In addition, our observations suggest that CMS in cotton is associated with the accelerated rates of rearrangement, and that novel expression products are derived from recombinant ORFs.

Molecular Evolution and Stress and Phytohormone Responsiveness of SUT Genes in Gossypium hirsutum.

Sucrose transporters (SUTs) play key roles in allocating the translocation of assimilates from source to sink tissues. Although the characteristics and biological roles of SUTs have been intensively investigated in higher plants, this gene family has not been functionally characterized in cotton. In this study, we performed a comprehensive analysis of SUT genes in the tetraploid cotton Gossypium hirsutum. A total of 18 G. hirsutum SUT genes were identified and classified into three groups based on their evolutionary relationships. Up to eight SUT genes in G. hirsutum were placed in the dicot-specific SUT1 group, while four and six SUT genes were, respectively, clustered into SUT4 and SUT2 groups together with members from both dicot and monocot species. The G. hirsutum SUT genes within the same group displayed similar exon/intron characteristics, and homologous genes in G. hirsutum At and Dt subgenomes, G. arboreum, and G. raimondii exhibited one-to-one relationships. Additionally, the duplicated genes in the diploid and polyploid cotton species have evolved through purifying selection, suggesting the strong conservation of SUT loci in these species. Expression analysis in different tissues indicated that SUT genes might play significant roles in cotton fiber elongation. Moreover, analyses of cis-acting regulatory elements in promoter regions and expression profiling under different abiotic stress and exogenous phytohormone treatments implied that SUT genes, especially GhSUT6A/D, might participate in plant responses to diverse abiotic stresses and phytohormones. Our findings provide valuable information for future studies on the evolution and function of SUT genes in cotton.

Zinc Finger-Homeodomain Transcriptional Factors (ZHDs) in Upland Cotton (Gossypium hirsutum): Genome-Wide Identification and Expression Analysis in Fiber Development.

Zinc finger-homeodomain (ZHD) genes encode a family of plant-specific transcription factors that not only participate in the regulation of plant growth and development but also play an important role in the response to abiotic stress. The ZHD gene family has been studied in several model plants, including Solanum lycopersicum, Zea mays, Oryza sativa, and Arabidopsis thaliana. However, a comprehensive study of the genes of the ZHD family and their roles in fiber development and pigmentation in upland cotton has not been completed. To address this gap, we selected a brown fiber cultivar for our study; brown color in cotton is one of the most desired colors in the textile industry. The natural colored fibers require less processing and little dying, thereby eliminating dye costs and chemical residues. Using bioinformatics approaches, we identified 37 GhZHD genes from Gossypium hirsutum and then divided these genes into seven groups based on their phylogeny. The GhZHD genes were mostly conserved in each subfamily with minor variations in motif distribution and gene structure. These genes were largely distributed on 19 of the 26 upland cotton chromosomes. Among the Gossypium genomes, the paralogs and orthologs of the GhZHD genes were identified and further characterized. Furthermore, among the paralogs, we observed that the ZHD family duplications in Gossypium genomes (G. hirsutum, G. arboreum, and G. raimondii) were probably derived from segmental duplication or genome-wide duplication (GWD) events. Through a combination of qRT-PCR and proanthocyanidins (PA) accumulation analyses in brown cotton fibers, we concluded that the candidate genes involved in early fiber development and fiber pigment synthesis include the following: GhZHD29, GhZHD35, GhZHD30, GhZHD31, GhZHD11, GhZHD27, GhZHD18, GhZHD15, GhZHD16, GhZHD22, GhZHD6, GhZHD33, GhZHD13, GhZHD5, and GhZHD23. This study delivers insights into the evolution of the GhZHD genes in brown cotton, serves as a valuable resource for further studies, and identifies the conditions necessary for improving the quality of brown cotton fiber.

Genome-wide identification and expression analyses of TCP transcription factor genes in Gossypium barbadense

Sea-island cotton (Gossypium barbadense) has drawn great attention in the textile industry for its comprehensive resistance and superior fiber properties. However, the mechanisms involved in fiber growth and development are unclear. As TCP transcription factors play important roles in plant growth and development, this study investigated the TCP family genes in G. barbadense (GbTCP). We identified 75 GbTCP genes, of which 68 had no introns. Phylogenetic analyses categorized the GbTCP transcription factors into 11 groups. Genomic analyses showed that 66 genes are located on 21 chromosomes. Phylogenetic analyses of G. arboreum, G. raimondii, G. hirsutum, G. barbadense, Theobroma cacao, Arabidopsis thaliana, Oryza sativa, Sorghum bicolor, and Zea mays, Picea abies, Sphagnum fallax and Physcomitrella patens, categorized 373 TCP genes into two classes (Classes I and II). By studying the structures of TCP genes in sea-island cotton, we identified genes from the same evolutionary branches that showed similar motif patterns. qRT-PCR results suggested that the GbTCPs had different expression patterns in fibers at various developmental stages of cotton, with several showing specific expression patterns during development. This report helps lay the foundation for future investigations of TCP functions and molecular mechanisms in sea-island cotton fiber development.

Identification and expression analysis of phosphatidy ethanolamine-binding protein (PEBP) gene family in cotton

The phosphatidy ethanolamine-binding proteins (PEBP) play an important role in controlling flower development and phase change. Here, a total of 61 PEBP genes were identified, in which 20, 21, 10, and 10 were from tetraploid Gossypium hirsutum (AD1) and G. barbadense (AD2), and diploid G. raimondii (D5) and G. arboreum (A2), respectively. In G. hirsutum, 20 identified PEBP genes were unevenly distributed on 12 chromosomes. The identified PEBP genes were classified into four groups (TFL1, MFT, FT and FT-like). Among those, FT-like group are unique to cotton. The majority of PEBP genes had similar intron/exon distribution, whereas the divergence of PEBP genes suggests the possibility of functional diversification. The expression of PEBP genes varied among different tissues. This study brings new insights into the integrated genome-wide identification of PEBP genes in cotton and provides a foundation for breeding cotton cultivars with early maturation.

Genome-Wide Analysis Elucidates the Role of CONSTANS-like Genes in Stress Responses of Cotton.

The CONSTANS (CO)-like gene family has been well studied for its role in the regulation of plant flowering time. However, their role remains poorly understood in cotton. To better understand the possible roles of CO-like in cotton, we performed a comprehensive genome-wide analysis of CO-like genes in cotton. Phylogenetic tree analysis showed that CO-like genes naturally clustered into three groups. Segmental duplication and whole genome duplication (WGD), which occurred before polyploidy, were important contributors to its expansion within the At (“t” indicates tetraploid) and Dt subgenomes, particularly in Group III. Long-terminal repeat retroelements were identified as the main transposable elements accompanying 18 genes. The genotype of GhCOL12_Dt displayed low diversity; it was a candidate involved in domestication. Selection pressure analyses indicated that relaxed purifying selection might have provided the main impetus during the evolution of CO-like genes in upland cotton. In addition, the high expression in the torus and calycle indicated that CO-like genes might affect flowering. The genes from Group II, and those from Group III involved in segmental duplication or WGD, might play important roles in response to drought and salt stress. Overall, this comprehensive genome-wide study of the CO-like gene family would facilitate further detailed studies in cotton.

Loop-Mediated Isothermal Amplification for Detection of Endogenous Sad1 Gene in Cotton: An Internal Control for Rapid Onsite GMO Testing.

Confirming the integrity of seed samples in powdered form is important prior to conducting a genetically modified organism (GMO) test. Rapid onsite methods may provide a technological solution to check for genetically modified (GM) events at ports of entry. In India, Bt cotton is the commercialized GM crop with four approved GM events; however, 59 GM events have been approved globally. GMO screening is required to test for authorized GM events. The identity and amplifiability of test samples could be ensured first by employing endogenous genes as an internal control. A rapid onsite detection method was developed for an endogenous reference gene, stearoyl acyl carrier protein desaturase (Sad1) of cotton, employing visual and real-time loop-mediated isothermal amplification (LAMP). The assays were performed at a constant temperature of 63°C for 30min for visual LAMP and 62ºC for 40min for real-time LAMP. Positive amplification was visualized as a change in color from orange to green on addition of SYBR® Green or detected as real-time amplification curves. Specificity of LAMP assays was confirmed using a set of 10 samples. LOD for visual LAMP was up to 0.1%, detecting 40 target copies, and for real-time LAMP up to 0.05%, detecting 20 target copies. The developed methods could be utilized to confirm the integrity of seed powder prior to conducting a GMO test for specific GM events of cotton. LAMP assays for the endogenous Sad1 gene of cotton have been developed to be used as an internal control for onsite GMO testing in cotton.

Retraction Note to: Overexpression of a cotton annexin gene, GhAnn1, enhances drought and salt stress tolerance in transgenic cotton

The authors are retracting this article because irregularities were involved in the reuse of stomatal images in Figure 9.

A genome-wide analysis of SWEET gene family in cotton and their expressions under different stresses

BackgroundThe SWEET (Sugars will eventually be exported transporters) gene family plays multiple roles in plant physiological activities and development process. It participates in reproductive development and in the process of sugar transport and absorption, plant senescence and stress responses and plant-pathogen interaction. However, thecomprehensive analysis of SWEET genes has not been reported in cotton.ResultsIn this study, we identified 22, 31, 55 and 60 SWEET genes from the sequenced genomes of Gossypium arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. Phylogenetic tree analysis showed that the SWEET genes could be divided into four groups, which were further classified into 14 sub-clades. Further analysis of chromosomal location, synteny analysis and gene duplication suggested that the orthologs showed a good collinearity and segmental duplication events played a crucial role in the expansion of the family in cotton. Specific MtN3_slv domains were highly conserved between Arabidopsis and cotton by exon-intron organization and motif analysis. In addition, the expression pattern in different tissues indicated that the duplicated genes in cotton might have acquired new functions as a result of sub-functionalization or neo-functionalization. The expression pattern of SWEET genes showed that the different genes were induced by diverse stresses. The identification and functional analysis of SWEET genes in cotton may provide more candidate genes for genetic modification.ConclusionSWEET genes were classified into four clades in cotton. The expression patterns suggested that the duplicated genes might have experienced a functional divergence. This work provides insights into the evolution of SWEET genes and more candidates for specific genetic modification, which will be useful in future research.

Characterization and functional analysis of GhWRKY42, a group IId WRKY gene, in upland cotton (Gossypium hirsutum L.)

BackgroundWRKY transcription factors (TFs) participate in various physiological processes of plants. Although WRKY genes have been well studied in model plants, knowledge of the functional roles of these genes is still extremely limited in cotton.ResultsIn this study, a group IId WRKY gene from cotton, GhWRKY42, was isolated and characterized. Our data showed that GhWRKY42 localized to the nucleus. A transactivation assay in yeast demonstrated that GhWRKY42 was not a transcriptional activator. A β-glucuronidase (GUS) activity assay revealed that the promoter of GhWRKY42 showed fragment deletion activity in Nicotiana tabacum and was mainly expressed in the roots, stems and leaves of ProGhWRKY42::GUS transgenic Arabidopsis plants. Quantitative real-time PCR (qRT-PCR) analysis indicated that GhWRKY42 was up-regulated during leaf senescence and was induced after exposure to abiotic stresses. Constitutive expression of GhWRKY42 in Arabidopsis led to a premature aging phenotype, which was correlated with an increased number of senescent leaves, reduced chlorophyll content and elevated expression of senescence-associated genes (SAGs). In addition, virus-induced gene silencing (VIGS) was used to silence the endogenous GhWRKY42 gene in cotton, and this silencing reduced plant height.ConclusionsOur findings indicate that GhWRKY42 is involved in abiotic stress responses, premature leaf senescence and stem development. This work establishes a solid foundation for further functional analysis of the GhWRKY42 gene in cotton.

Genome-wide analysis of nucleotide binding site-leucine-rich repeats (NBS-LRR) disease resistance genes in Gossypium hirsutum

Nucleotide binding site-leucine-rich repeats (NBS-LRR) disease resistance genes (R genes) are the main class of R genes in plants. In order to fully understand the structure characteristic and composition of NBS-LRR R genes in upland cotton, and further reveal their roles in against to various pathogens, all NBS-LRR R genes in upland cotton were gathered, then their physical located distribution on chromosomes were analyzed. Next, all these upland cotton NBS-LRR R genes were performed by homologous comparison analysis compared with Arabidopsis thaliana (A. thaliana) NBS-LRR R genes. The results showed that cotton NBS-LRR R genes are distributed on 26 chromosomes with clustered or single permutation characters. These R genes can be divided into two subfamilies of group CC-NBS-LRR (315 genes) and TIR-NBS-LRR (122 genes). The CC-NBS-LRR subtype can be further classified into five groups: A, B, C, D and NL-A. The TIR-NBS-LRR subtype contains seven groups: B, C, D, E, G, TN-A and putative genes. Some cotton NBS-LRR R genes have homologous genes with those of A. thaliana, but most of them presented obvious differences compared with Arabidopsis NBS-LRR R genes. The basic information on the classification of cotton NBS-LRR R genes would be helpful for further ongoing their functional exploration and finally revealing their roles in against to various pathogens.

Genome-wide characterization, identification, and expression analysis of the WD40 protein family in cotton.

WD40 repeat proteins are largely distributed across the plant kingdom and play an important role in diverse biological activities. In this work, we performed genome-wide identification, characterization, and expression level analysis of WD40 genes in cotton. A total of 579, 318, and 313 WD40 genes were found in Gossypium hirsutum, G. arboreum, and G. raimondii, respectively. Based on phylogenetic tree analyses, WD40 genes were divided into 11 groups with high similarities in exon/intron features and protein domains within the group. Expression analysis of WD40 genes showed differential expression at different stages of cotton fiber development (0 and 8 DPA) and cotton stem. A number of miRNAs were identified to target WD40 genes that are significantly involved in cotton fiber development during the initiation and elongation stages. These include miR156, miR160, miR162, miR164, miR166, miR167, miR169, miR171, miR172, miR393, miR396, miR398, miR2950, and miR7505. The findings provide a stronger indication of WD40 gene function and their involvement in the regulation of cotton fiber development during the initiation and elongation stages.

Genome-Wide Analysis of Multidrug and Toxic Compound Extrusion (MATE) Family in Gossypium raimondii and Gossypium arboreum and Its Expression Analysis Under Salt, Cadmium, and Drought Stress.

The extrusion of toxins and substances at a cellular level is a vital life process in plants under abiotic stress. The multidrug and toxic compound extrusion (MATE) gene family plays a large role in the exportation of toxins and other substrates. We carried out a genome-wide analysis of MATE gene families in Gossypium raimondii and Gossypium arboreum and assessed their expression levels under salt, cadmium and drought stresses. We identified 70 and 68 MATE genes in G. raimondii and G. arboreum, respectively. The majority of the genes were predicted to be localized within the plasma membrane, with some distributed in other cell parts. Based on phylogenetic analysis, the genes were subdivided into three subfamilies, designated as M1, M2 and M3. Closely related members shared similar gene structures, and thus were highly conserved in nature and have mainly evolved through purifying selection. The genes were distributed in all chromosomes. Twenty-nine gene duplication events were detected, with segmental being the dominant type. GO annotation revealed a link to salt, drought and cadmium stresses. The genes exhibited differential expression, with GrMATE18, GrMATE34, GaMATE41 and GaMATE51 significantly upregulated under drought, salt and cadmium stress, and these could possibly be the candidate genes. Our results provide the first data on the genome-wide and functional characterization of MATE genes in diploid cotton, and are important for breeders of more stress-tolerant cotton genotypes.

Genome- Wide Analysis and Characterization of the TRX Gene Family in Upland Cotton

Thioredoxins (TRX) are small molecules of proteins that are present in all organisms. TRXs play an important role in diverse functions of plant growth and development. In this study, we performed genome-wide, characterization and expression levels of TRX gene family in cotton. A total of 150 GhTRX proteins were identified in upland cotton and classified into five subfamilies based on their domain compositions. Phylogenetic tree analysis divided TRX genes into seven subgroups. GhTRX genes covered all upland cotton chromosomes, with duplicated gene events. Ka/Ks ratio of three gene pairs was less than 1, suggesting purifying selection. The functions of GhTRX genes were studied using gene ontology, protein localization, and promoter analysis. Furthermore, six GhTRX genes were randomly selected to examine their expression level in cotton development and under various exogenous treatments. The genes showed high expressions in various tissues and at different stages of leaf senescence, also showed high expression under abscisic acid, ethylene, drought, and salinity. This study reveals the first report of TRX family genes in upland cotton. However further studies are needed to elucidate their specific functions in cotton plant.

Changes in DNA methylation assessed by genomic bisulfite sequencing suggest a role for DNA methylation in cotton fruiting branch development.

Cotton plant architecture, including fruit branch formation and flowering pattern, influences plant light exploitation, cotton yield and planting cost. DNA methylation has been widely observed at different developmental stages in both plants and animals and is associated with regulation of gene expression, chromatin remodelling, genome protection and other functions. Here, we investigated the global epigenetic reprogramming during the development of fruiting branches and floral buds at three developmental stages: the seedling stage, the pre-squaring stage and the squaring stage. We first identified 22 cotton genes which potentially encode DNA methyltransferases and demethylases. Among them, the homologous genes of CMT, DRM2 and MET1 were upregulated at pre-squaring and squaring stages, suggesting that DNA methylation is involved in the development of floral buds and fruit branches. Although the global methylation at all of three developmental stages was not changed, the CHG-type methylation of non-expressed genes was higher than those of expressed genes. In addition, we found that the expression of the homologous genes of the key circadian rhythm regulators, including CRY, LHY and CO, was associated with changes of DNA methylation at three developmental stages.