Cellular programming of plant gene imprinting.

Gene imprinting, the differential expression of maternal and paternal alleles, independently evolved in mammals and in flowering plants. A unique feature of flowering plants is a double-fertilization event in which the sperm fertilize not only the egg, which forms the embryo, but also the central cell, which develops into the endosperm (an embryo-supporting tissue). The distinctive mechanisms of gene imprinting in the endosperm, which involve DNA demethylation and histone methylation, begin in the central cell and sperm prior to fertilization. Flowering plants might have coevolved double fertilization and imprinting to prevent parthenogenetic development of the endosperm.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae.

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

The catalytic core of DEMETER guides active DNA demethylation in Arabidopsis.

The Arabidopsis DEMETER (DME) DNA glycosylase demethylates the maternal genome in the central cell prior to fertilization and is essential for seed viability. DME preferentially targets small transposons that flank coding genes, influencing their expression and initiating plant gene imprinting. DME also targets intergenic and heterochromatic regions, but how it is recruited to these differing chromatin landscapes is unknown. The C-terminal half of DME consists of 3 conserved regions required for catalysis in vitro. We show that this catalytic core guides active demethylation at endogenous targets, rescuing dme developmental and genomic hypermethylation phenotypes. However, without the N terminus, heterochromatin demethylation is significantly impeded, and abundant CG-methylated genic sequences are ectopically demethylated. Comparative analysis revealed that the conserved DME N-terminal domains are present only in flowering plants, whereas the domain architecture of DME-like proteins in nonvascular plants mainly resembles the catalytic core, suggesting that it might represent the ancestral form of the 5mC DNA glycosylase found in plant lineages. We propose a bipartite model for DME protein action and suggest that the DME N terminus was acquired late during land plant evolution to improve specificity and facilitate demethylation at heterochromatin targets.

Application of 5-Methylcytosine DNA Glycosylase to the Quantitative Analysis of DNA Methylation.

In higher eukaryotes DNA methylation is a prominent epigenetic mark important for chromatin structure and gene expression. Thus, profiling DNA methylation is important for predicting gene expressions associated with specific traits or diseases. DNA methylation is achieved by DNA methyltransferases and can be actively removed by specific enzymes in a replication-independent manner. DEMETER (DME) is a bifunctional 5-methylcytosine (5mC) DNA glycosylase responsible for active DNA demethylation that excises 5mC from DNA and cleaves a sugar-phosphate bond generating a single strand break (SSB). In this study, DME was used to analyze DNA methylation levels at specific epialleles accompanied with gain or loss of DNA methylation. DME treatment on genomic DNA generates SSBs in a nonsequence-specific fashion proportional to 5mC density, and thus DNA methylation levels can be easily measured when combined with the quantitative PCR (qPCR) method. The DME-qPCR analysis was applied to measure DNA methylation levels at the FWA gene in late-flowering Arabidopsis mutants and the CNR gene during fruit ripening in tomato. Differentially methylated epialleles were successfully distinguished corresponding to their expression levels and phenotypes. DME-qPCR is proven a simple yet effective method for quantitative DNA methylation analysis, providing advantages over current techniques based on methylation-sensitive restriction digestion.

Revealing biomass heterosis in the allodiploid xBrassicoraphanus, a hybrid between Brassica rapa and Raphanus sativus, through integrated transcriptome and metabolites analysis.

BACKGROUND: Heterosis is biologically important but the molecular basis of the phenomenon is poorly understood. We characterized intergeneric hybrids between B. rapa cv. Chiifu and R. sativus cv. WK10039 as an extreme example of heterosis. Taking advantage of clear heterosis phenotypes and the genetic distance between parents, we performed transcriptome and metabolite analysis to decipher the molecular basis of heterosis. RESULTS: The heterosis was expressed as fresh weight in the field and as inflorescence stem length in the glass house. Flowering time, distributed as a normal segregating population, ranged from the early flowering of one parent to the late flowering of the other, in contrast to the homogeneous flowering time in a typical F1 population, indicating unstable allelic interactions. The transcriptome and metabolome both indicated that sugar metabolism was altered, suggesting that the change in metabolism was linked to the heterosis. Because alleles were not shared between the hybridized genomes, classic models only partly explain this heterosis, indicating that other mechanisms are involved. CONCLUSION: The differential expression of genes for primary and secondary metabolism, along with the altered metabolite profiles, suggests that heterosis could involve a change in balance between primary and secondary metabolism.

ROS1-Dependent DNA Demethylation Is Required for ABA-Inducible NIC3 Expression.

DNA methylation plays an important role in diverse developmental processes in many eukaryotes, including the response to environmental stress. Abscisic acid (ABA) is a plant hormone that is up-regulated under stress. The involvement of DNA methylation in the ABA response has been reported but is poorly understood. DNA demethylation is a reverse process of DNA methylation and often induces structural changes of chromatin leading to transcriptional activation. In Arabidopsis (Arabidopsis thaliana), active DNA demethylation depends on the activity of REPRESSOR OF SILENCING 1 (ROS1), which directly excises 5-methylcytosine from DNA. Here we showed that ros1 mutants were hypersensitive to ABA during early seedling development and root elongation. Expression levels of some ABA-inducible genes were decreased in ros1 mutants, and more than 60% of their proximal regions became hypermethylated, indicating that a subset of ABA-inducible genes are under the regulation of ROS1-dependent DNA demethylation. Notable among them is NICOTINAMIDASE 3 (NIC3) that encodes an enzyme that converts nicotinamide to nicotinic acid in the NAD+ salvage pathway. Many enzymes in this pathway are known to be involved in stress responses. The nic3 mutants display hypersensitivity to ABA, whereas overexpression of NIC3 restores normal ABA responses. Our data suggest that NIC3 is responsive to ABA but requires ROS1-mediated DNA demethylation at the promoter as a prerequisite to transcriptional activation. These findings suggest that ROS1-induced active DNA demethylation maintains the active state of NIC3 transcription in response to ABA.

Deficiency of rice hexokinase HXK5 impairs synthesis and utilization of starch in pollen grains and causes male sterility.

There is little known about the function of rice hexokinases (HXKs) in planta. We characterized hxk5-1, a Tos17 mutant of OsHXK5 that is up-regulated in maturing pollen, a stage when starch accumulates. Progeny analysis of self-pollinated heterozygotes of hxk5-1 and reciprocal crosses between the wild-type and heterozygotes revealed that loss of HXK5 causes male sterility. Homozygous hxk5-1, produced via anther culture, and additional homozygous hxk5-2, hxk5-3 and hxk5-4 lines created by CRISPR/Cas9 confirmed the male-sterile phenotype. In vitro pollen germination ability and in vivo pollen tube growth rate were significantly reduced in the hxk5 mutant pollen. Biochemical analysis of anthers with the mutant pollen revealed significantly reduced hexokinase activity and starch content, although they were sufficient to produce some viable seed. However, the mutant pollen was unable to compete successfully against wild-type pollen. Expression of the catalytically inactive OsHXK5-G113D did not rescue the hxk5 male-sterile phenotype, indicating that its catalytic function was responsible for pollen fertility, rather than its role in sugar sensing and signaling. Our results demonstrate that OsHXK5 contributes to a large portion of the hexokinase activity necessary for the starch utilization pathway during pollen germination and tube growth, as well as for starch biosynthesis during pollen maturation.

AP endonucleases process 5-methylcytosine excision intermediates during active DNA demethylation in Arabidopsis.

DNA methylation is a primary epigenetic modification regulating gene expression and chromatin structure in many eukaryotes. Plants have a unique DNA demethylation system in that 5-methylcytosine (5mC) is directly removed by DNA demethylases, such as DME/ROS1 family proteins, but little is known about the downstream events. During 5mC excision, DME produces 3'-phosphor-α, ß-unsaturated aldehyde and 3'-phosphate by successive ß- and δ-eliminations, respectively. The kinetic studies revealed that these 3'-blocking lesions persist for a significant amount of time and at least two different enzyme activities are required to immediately process them. We demonstrate that Arabidopsis AP endonucleases APE1L, APE2 and ARP have distinct functions to process such harmful lesions to allow nucleotide extension. DME expression is toxic to E. coli due to excessive 5mC excision, but expression of APE1L or ARP significantly reduces DME-induced cytotoxicity. Finally, we propose a model of base excision repair and DNA demethylation pathway unique to plants.

5-methylcytosine recognition by Arabidopsis thaliana DNA glycosylases DEMETER and DML3.

Methylation of cytosine to 5-methylcytosine (5mC) is important for gene expression, gene imprinting, X-chromosome inactivation, and transposon silencing. Active demethylation in animals is believed to proceed by DNA glycosylase removal of deaminated or oxidized 5mC. In plants, 5mC is removed from the genome directly by the DEMETER (DME) family of DNA glycosylases. Arabidopsis thaliana DME excises 5mC to activate expression of maternally imprinted genes. Although the related Repressor of Silencing 1 (ROS1) enzyme has been characterized, the molecular basis for 5mC recognition by DME has not been investigated. Here, we present a structure-function analysis of DME and the related DME-like 3 (DML3) glycosylases for 5mC and its oxidized derivatives. Relative to 5mC, DME and DML3 exhibited robust activity toward 5-hydroxymethylcytosine, limited activity for 5-carboxylcytosine, and no activity for 5-formylcytosine. We used homology modeling and mutational analysis of base excision and DNA binding to identify residues important for recognition of 5mC within the context of DNA and inside the enzyme active site. Our results indicate that the 5mC binding pocket is composed of residues from discrete domains and is responsible for discrimination against 5mC derivatives, and suggest that DME, ROS1, and DML3 utilize subtly different mechanisms to probe the DNA duplex for cytosine modifications.

Excision of 5-hydroxymethylcytosine by DEMETER family DNA glycosylases.

In plants and animals, 5-methylcytosine (5mC) serves as an epigenetic mark to repress gene expression, playing critical roles for cellular differentiation and transposon silencing. Mammals also have 5-hydroxymethylcytosine (5hmC), resulting from hydroxylation of 5mC by TET family-enzymes. 5hmC is abundant in mouse Purkinje neurons and embryonic stem cells, and regarded as an important intermediate for active DNA demethylation in mammals. However, the presence of 5hmC in plants has not been clearly demonstrated. In Arabidopsis, the DEMETER (DME) family DNA glycosylases efficiently remove 5mC, which results in DNA demethylation and transcriptional activation of target genes. Here we show that DME and ROS1 have a significant 5hmC excision activity in vitro, although we detected no 5hmC in Arabidopsis, suggesting that it is very unlikely for plants to utilize 5hmC as a DNA demethylation intermediate. Our results indicate that both plants and animals have 5mC in common but DNA demethylation systems have independently evolved with distinct mechanisms.

Domain structure of the DEMETER 5-methylcytosine DNA glycosylase.

DNA glycosylases initiate the base excision repair (BER) pathway by excising damaged, mismatched, or otherwise modified bases. Animals and plants independently evolved active BER-dependent DNA demethylation mechanisms important for epigenetic reprogramming. One such DNA demethylation mechanism is uniquely initiated in plants by DEMETER (DME)-class DNA glycosylases. Arabidopsis DME family glycosylases contain a conserved helix-hairpin-helix domain present in both prokaryotic and eukaryotic DNA glycosylases as well as two domains A and B of unknown function that are unique to this family. Here, we employed a mutagenesis approach to screen for DME residues critical for DNA glycosylase activity. This analysis revealed that amino acids clustered in all three domains, but not in the intervening variable regions, are required for in vitro 5-methylcytosine excision activity. Amino acids in domain A were found to be required for nonspecific DNA binding, a prerequisite for 5-methylcytosine excision. In addition, mutational analysis confirmed the importance of the iron-sulfur cluster motif to base excision activity. Thus, the DME DNA glycosylase has a unique structure composed of three essential domains that all function in 5-methylcytosine excision.

Epigenome editing: targeted manipulation of epigenetic modifications in plants.

BACKGROUND: Epigenetic modifications play important roles in diverse cellular processes such as X chromosome inactivation, cell differentiation, development and senescence. DNA methylation and histone modifications are major epigenetic modifications that regulate chromatin structure and gene expression without DNA sequence changes. Epigenetic alterations may induce phenotypic changes stable enough for mitotic or meiotic inheritance. Moreover, the reversibility of epigenetic marks makes the manipulation of chromatin and epigenetic signature an attractive strategy for therapeutic and breeding purposes. Targeted epigenetic manipulation, or epigenome editing, at the gene of interest commonly utilizes specific epigenetic modifiers fused with a targeting module of the conventional genome editing system. OBJECTIVE: This review aims to summarize essential epigenetic components and introduce currently available epigenetic mutants and the corresponding epialleles in plants. Furthermore, advances in epigenome editing technology are discussed while proposing its potential application to plant breeding. CONCLUSIONS: Epimutations associated with useful traits may provide a valuable resource for crop development. It is important to explore epimutations in a variety of crop species while understanding the fundamental aspects of epigenetic regulation of agronomically important traits such as yield, quality, disease resistance and stress tolerance. In the end, plant breeding programs through epigenome editing may help not only to expand the use of limited genetic resources but also to alleviate consumers' concerns about genetically manipulated crops.

Dynamic changes in DNA methylation occur in TE regions and affect cell proliferation during leaf-to-callus transition in Arabidopsis.

Plant somatic cells can be reprogrammed into pluripotent cell mass, called callus, through a two-step in vitro tissue culture method. Incubation on callus-inducing medium triggers active cell proliferation to form a pluripotent callus. Notably, DNA methylation is implicated during callus formation, but a detailed molecular process regulated by DNA methylation remains to be fully elucidated. Here, we compared genome-wide DNA methylation profiles between leaf and callus tissues in Arabidopsis using whole-genome bisulphite-sequencing. Global distribution of DNA methylation showed that CHG methylation was increased, whereas CHH methylation was reduced especially around transposable element (TE) regions during the leaf-to-callus transition. We further analysed differentially expressed genes around differentially methylated TEs (DMTEs) during the leaf-to-callus transition and found that genes involved in cell cycle regulation were enriched and also constituted a coexpression gene network along with pluripotency regulators. In addition, a conserved DNA sequence analysis for upstream cis-elements led us to find a putative transcription factor associated with cell fate transition. CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) was newly identified as a regulator of plant regeneration, and consistently, the cca1lhy mutant displayed altered phenotypes in callus proliferation. Overall, these results suggest that DNA methylation coordinates cell cycle regulation during callus formation, and CCA1 may act as a key upstream coordinator at least in part in the processes.

Endosperm gene imprinting and seed development.

Imprinting occurs in the endosperm of flowering plants. Endosperm, produced by fertilization of the central cell in the female gametophyte, is essential for embryo and seed development. Several imprinted genes play an important role in endosperm development. The mechanism of gene imprinting involves DNA methylation and histone modification. DNA methylation is actively removed at the imprinted alleles to be activated. Histone methylation mediated by the Polycomb group complex provides another layer of epigenetic regulation at the silenced alleles. Endosperm gene imprinting can be uncoupled from seed development when fertilization of the central cell is prevented. Imprinting may be a mechanism to ensure fertilization of the central cell thereby preventing parthenogenic development of the endosperm.

Subgenome Discrimination in Brassica and Raphanus Allopolyploids Using Microsatellites.

Intergeneric crosses between Brassica species and Raphanus sativus have produced crops with prominent shoot and root systems of Brassica and R. sativus, respectively. It is necessary to discriminate donor genomes when studying cytogenetic stability in distant crosses to identify homologous chromosome pairing, and microsatellite repeats have been used to discriminate subgenomes in allopolyploids. To identify genome-specific microsatellites, we explored the microsatellite content in three Brassica species (B. rapa, AA, B. oleracea, CC, and B. nigra, BB) and R. sativus (RR) genomes, and validated their genome specificity by fluorescence in situ hybridization. We identified three microsatellites showing A, C, and B/R genome specificity. ACBR_msat14 and ACBR_msat20 were detected in the A and C chromosomes, respectively, and ACBR_msat01 was detected in B and R genomes. However, we did not find a microsatellite that discriminated the B and R genomes. The localization of ACBR_msat20 in the 45S rDNA array in ×Brassicoraphanus 977 corroborated the association of the 45S rDNA array with genome rearrangement. Along with the rDNA and telomeric repeat probes, these microsatellites enabled the easy identification of homologous chromosomes. These data demonstrate the utility of microsatellites as probes in identifying subgenomes within closely related Brassica and Raphanus species for the analysis of genetic stability of new synthetic polyploids of these genomes.

Reduced fertility caused by meiotic defects and micronuclei formation during microsporogenesis in xBrassicoraphanus.

BACKGROUND: Hybridization and polyploidization events are important driving forces in plant evolution. Allopolyploids formed between different species can be naturally or artificially created but often suffer from genetic instability and infertility in successive generations. xBrassicoraphanus is an intergeneric allopolyploid obtained from a cross between Brassica rapa and Raphanus sativus, providing a useful resource for genetic and genomic study in hybrid species. OBJECTIVE: The current study aims to understand the cause of hybrid sterility and pollen abnormality in different lines of synthetic xBrassicoraphanus from the cytogenetic perspective. METHODS: Alexander staining was used to assess the pollen viability. Cytogenetic analysis was employed to monitor meiotic chromosome behaviors in pollen mother cells (PMCs). Origins of parental chromosomes in xBrassicoraphanus meiocytes were determined by genome in situ hybridization analysis. RESULTS: The xBrassicoraphanus lines BB#4 and BB#6 showed high rates of seed abortion and pollen deformation. Abnormal chromosome behaviors were observed in their PMCs, frequently forming univalents and inter-chromosomal bridges during meiosis. A positive correlation also exists between meiotic defects and the formation of micronuclei, which is conceivably responsible for unbalanced gamete production and pollen sterility. CONCLUSION: These results suggest that unequal segregation of meiotic chromosomes, due in part to non-homologous interactions, is responsible for micronuclei and unbalanced gamete formation, eventually leading to pollen degeneration and inferior fertility in unstable xBrassicoraphanus lines.

Meiotic Chromosome Stability and Suppression of Crossover Between Non-homologous Chromosomes in xBrassicoraphanus, an Intergeneric Allotetraploid Derived From a Cross Between Brassica rapa and Raphanus sativus.

Hybridization and polyploidization are major driving forces in plant evolution. Allopolyploids can be occasionally formed from a cross between distantly related species but often suffer from chromosome instability and infertility. xBrassicoraphanus is an intergeneric allotetraploid (AARR; 2n = 38) derived from a cross between Brassica rapa (AA; 2n = 20) and Raphanus sativus (RR; 2n = 18). xBrassicoraphanus is fertile and genetically stable, while retaining complete sets of both B. rapa and R. sativus chromosomes. Precise control of meiotic recombination is essential for the production of balanced gametes, and crossovers (COs) must occur exclusively between homologous chromosomes. Many interspecific hybrids have problems with meiotic division at early generations, in which interactions between non-homologous chromosomes often bring about aneuploidy and unbalanced gamete formation. We analyzed meiotic chromosome behaviors in pollen mother cells (PMCs) of allotetraploid and allodiploid F1 individuals of newly synthesized xBrassicoraphanus. Allotetraploid xBrassicoraphanus PMCs showed a normal diploid-like meiotic behavior. By contrast, allodiploid xBrassicoraphanus PMCs displayed abnormal segregation of chromosomes mainly due to the absence of homologous pairs. Notably, during early stages of meiosis I many of allodiploid xBrassicoraphanus chromosomes behave independently with few interactions between B. rapa and R. sativus chromosomes, forming many univalent chromosomes before segregation. Chromosomes were randomly assorted at later stages of meiosis, and tetrads with unequal numbers of chromosomes were formed at completion of meiosis. Immunolocalization of HEI10 protein mediating meiotic recombination revealed that COs were more frequent in synthetic allotetraploid xBrassicoraphanus than in allodiploid, but less than in the stabilized line. These findings suggest that structural dissimilarity between B. rapa and R. sativus chromosomes prevents non-homologous interactions between the parental chromosomes in allotetraploid xBrassicoraphanus, allowing normal diploid-like meiosis when homologous pairing partners are present. This study also suggests that CO suppression between non-homologous chromosomes is required for correct meiotic progression in newly synthesized allopolyploids, which is important for the formation of viable gametes and reproductive success in the hybrid progeny.

Purification of human transcription factors Nanog and Sox2, each in complex with Skp, an Escherichia coli periplasmic chaperone.

Nanog and Sox2 are key transcriptional factors involved in self-renewal and pluripotency of stem cells in human and other mammals. Nanog and Sox2 contain homeodomain (HD) and high-mobility group (HMG) DNA-binding domain, respectively, for targeting them to their regulatory regions and the other regions with transactivation function by providing sites for recruiting other transcriptional regulators. To gain insights in the biochemical and biophysical characteristics of the other regions of Nanog and Sox2, we have tried to overproduce and purify full length wild-type human Nanog and Sox2 expressed in Escherichia coli. Interestingly, we found that Nanog and Sox2 were individually stabilized by tight interaction with Skp, an E. coli periplasmic chaperone, thereby enabling stable over-expression and purification of Nanog and Sox2, each in complex with Skp. Purified Skp complexes of Nanog and Sox maintained DNA-binding activity toward its cognate DNA sequence. A similar approach may be applicable for some other mammalian proteins that are unstable or difficult to over-express in E. coli.

A genetic characterization of Korean waxy maize (Zea mays L.) landraces having flowering time variation by RNA sequencing.

Maize is the second-most produced crop in the Korean peninsula and has been continuously cultivated since the middle of the 16th century, when it was originally introduced from China. Even with this extensive cultivation history, the diversity and properties of Korean landraces have not been investigated at the nucleotide sequence level. We collected 12 landraces with various flowering times and performed RNA-seq in the early vegetative stage. The transcriptomes of 12 Korean landraces have been analyzed for their genetic variations in coding sequence and genetic relationships to other maize germplasm. The Korean landraces showed specific genetic characteristics and were closely related to a Chinese inbred line. Flowering-time related gene profiles pointed to multiple causes for the variation of flowering time within Korean landraces; the profiles revealed significant positive and negative correlations among genes, allowing us to infer possible mechanisms for flowering time variation in maize. Our results demonstrate the value of transcriptome-based genetic and gene expression profiles for information on possible breeding resources, which is particularly needed in Korean waxy landraces.

Transcriptional regulation of abscisic acid biosynthesis and signal transduction, and anthocyanin biosynthesis in 'Bluecrop' highbush blueberry fruit during ripening.

Highbush blueberry (Vaccinium corymbosum) fruit accumulate high levels of anthocyanins during ripening, which might be controlled by abscisic acid (ABA), a signal molecule in non-climacteric fruits. For an integrated view of the ripening process from ABA to anthocyanin biosynthesis, we analyzed the transcriptomes of 'Bluecrop' highbush blueberry fruit using RNA-Seq at three ripening stages, categorized based on fruit skin coloration: pale green at ca. 30 days after full bloom (DAFB), reddish purple at ca. 40 DAFB, and dark purple at ca. 50 DAFB. Mapping the trimmed reads against the reference sequences yielded 25,766 transcripts. Of these, 143 transcripts were annotated to encode five ABA biosynthesis enzymes, four ABA signal transduction regulators, four ABA-responsive transcription factors, and 12 anthocyanin biosynthesis enzymes. The analysis of differentially expressed genes between the ripening stages revealed that 11 transcripts, including those encoding nine-cis-epoxycarotenoid dioxygenase, SQUAMOSA-class MADS box transcription factor, and flavonoid 3',5'-hydroxylase, were significantly up-regulated throughout the entire ripening stages. In fruit treated with 1 g L-1 ABA, at least nine transcripts of these 11 transcripts as well as one transcript encoding flavonoid 3'-hydroxylase were up-regulated, presumably promoting anthocyanin accumulation and fruit skin coloration. These results will provide fundamental information demonstrating that ABA biosynthesis and signal transduction, and anthocyanin biosynthesis are closely associated with anthocyanin accumulation and skin coloration in highbush blueberry fruit during ripening.