A sensitive transcriptome analysis method that can detect unknown transcripts

Characterization of Chinese hamster ovary cells with disparate chromosome numbers: Reduction of the amount of mRNA relative to total protein

Chromosomes in Chinese hamster ovary (CHO) cells are labile. We have shown that high-chromosome-number CHO cells have greater potential to become robust producers of recombinant proteins. One explanation being the increase in transgene integration sites. However, high-chromosome-number cell clones produce more IgG3 following culture of single-cell clones, even under conditions that yield the same number of integrations as cells with normal chromosome numbers. Here, we characterized high-chromosome-number cells by transcriptome analysis. RNA standards were used to normalize transcriptomes of cells that had different chromosome numbers. Our results demonstrate that the mRNA ratio of β-actin and many other genes in high-chromosome-number cells to that in normal-chromosome-number cells per cell (normalized to RNA standards) was smaller than the equivalent genomic size and cell volume ratios. Many genes encoding membrane proteins are more highly expressed in high-chromosome-number cells, probably due to differences in cell size caused by the increase in chromosomes. In addition, genes related to histone modification and lipid metabolism are differentially expressed. The reduced transcript level required per protein produced in total and the different intracellular signal transductions might be key factors for antibody production.

Enhanced IgG1 production by overexpression of nuclear factor kappa B inhibitor zeta (NFKBIZ) in Chinese hamster ovary cells

Several engineering strategies have been employed to improve the production of therapeutic recombinant proteins in Chinese hamster ovary (CHO) cell lines. We have focused on unfolded protein response-based engineering and reported that ATF4 overexpression increases protein production. In this study, transcriptome analysis of ATF4-overexpressed CHO cells was performed using high-coverage expression profiling, to search for another key factor contributing to recombinant protein production. We observed the upregulated expression of transcription factor, nuclear factor (NF)-kappa-B inhibitor zeta (NFKBIZ or Iκbζ), in ATF4-overexpressed cells. A total of 1917 bp of CHO NFKBIZ cDNA was cloned, and two stable cell lines overexpressing NFKBIZ were constructed. We investigated the effects of NFKBIZ on IgG1 production in CHO cells. Although the two stable cell lines, NFKBIZ-A and -B, had the opposite phenotypes in cell growth, the specific IgG1 production rate of both cell lines was enhanced by 1.2-1.4-fold. In the NFKBIZ-A cell line, the synergistic effect between enhanced viable cell density and improved specific IgG1 production rate brought about a large increase in the final IgG1 titer. Luciferase-based NF-κB signaling assay results suggest that altered p50/p50 signaling seems to be due to the opposite phenotypes in cell growth. No difference was observed in the translational levels and intracellular assembly states of IgG1 between mock and two NFKBIZ cell lines, indicating that the secretion machinery of correctly folded IgG1 was enhanced in NFKBIZ-overexpressing cell lines.

Discovery of novel cold-induced CISP genes encoding small RNA-binding proteins related to cold adaptation in barley

To adapt to cold conditions, barley plants rely on specific mechanisms, which have not been fully understood. In this study, we characterized a novel barley cold-induced gene identified using a PCR-based high coverage gene expression profiling method. The identified gene encodes a small protein that we named CISP1 (Cold-induced Small Protein 1). Homology searches of sequence databases revealed that CISP1 homologs (CISP2 and CISP3) exist in barley genome. Further database analyses showed that the CISP1 homologs were widely distributed in cold-tolerant plants such as wheat and rye. Quantitative reverse transcription PCR analyses indicated that the expression of barley CISP genes was markedly increased in roots exposed to cold conditions. In situ hybridization analyses showed that the CISP1 transcripts were localized in the root tip and lateral root primordium. We also demonstrated that the CISP1 protein bound to RNA. Taken together, these findings indicate that CISP1 and its homologs encoding small RNA-binding proteins may serve as RNA chaperones playing a vital role in the cold adaptation of barley root. This is the first report describing the likely close relationship between root-specific genes and the cold adaptation process, as well as the potential function of the identified genes.

Differential transcript profiling alters regulatory gene expression during the development of Gossypium arboreum, G.stocksii and somatic hybrids

Polyploidy or genome doubling (i.e., the presence of two or more diploid parental genome sets within an organism) are very important in higher plants. Of particular interest are the mechanisms in the new microenvironment of the common nucleus, where doubled regulatory networks interact to generate a viable genetic system capable of regulating growth, development and responses to the environment. To determine the effects of whole genome merging and doubling on the global gene expression architecture of a new polyploid, derived from protoplast fusion of the A₁A₁ genome of Gossypium arboreum and the E₁E₁ genome of Gossypium stocksii, we monitored gene expression through cDNA-AFLP in the somatic hybrids (G. arboreum + G. stocksii). The genomic expression patterns of the somatic hybrids revealed that changes in expression levels mainly involved regulatory genes (31.8% of the gene expression profiles), and the AA and EE genomes contributed equally to genome-wide expression in the newly formed AAEE genome from additivity and dominance effects. These results provide a novel perspective on polyploid gene regulation and hint at the underlying genetic basis of allopolyploid adaption in the new microenvironmental nucleus.

Application of cDNA-AFLP to biomarker exploration in a non-model species Grandidierella japonica

Biomarkers of exposure can be used to identify specific contaminants that are adversely affecting aquatic organisms. However, it remains prohibitively costly to investigate multiple novel biomarkers of exposure in a non-model species, despite the development of next-generation sequencing technology. In this study, we focused on the use of cDNA-amplified fragment length polymorphism (AFLP) as a cost-effective biomarker discovery tool to test whether it could identify biomarkers of exposure in the non-model amphipod species Grandidierella japonica. Loci were identified that were differentially expressed in amphipods exposed to reference chemicals (Cu, Zn, and nicotine) and to an environmental sample (road dust) at sublethal concentrations. Eight loci were shown to respond consistently to nicotine at different concentrations, but not to Cu or Zn. Some of the loci also responded to an environmental road dust sample containing nicotine. These findings suggest that loci identified using cDNA-AFLP could be used as biomarkers of nicotine exposure in environmental samples with complex matrices. Further studies with other organisms and toxicants are needed, but we have demonstrated that the use of cDNA-AFLP to identify biomarkers for ecotoxicological studies of non-model species is at least feasible.

Methylated site display (MSD)-AFLP, a sensitive and affordable method for analysis of CpG methylation profiles

Background

It has been pointed out that environmental factors or chemicals can cause diseases that are developmental in origin. To detect abnormal epigenetic alterations in DNA methylation, convenient and cost-effective methods are required for such research, in which multiple samples are processed simultaneously. We here present methylated site display (MSD), a unique technique for the preparation of DNA libraries. By combining it with amplified fragment length polymorphism (AFLP) analysis, we developed a new method, MSD-AFLP.

Results

Methylated site display libraries consist of only DNAs derived from DNA fragments that are CpG methylated at the 5′ end in the original genomic DNA sample. To test the effectiveness of this method, CpG methylation levels in liver, kidney, and hippocampal tissues of mice were compared to examine if MSD-AFLP can detect subtle differences in the levels of tissue-specific differentially methylated CpGs. As a result, many CpG sites suspected to be tissue-specific differentially methylated were detected. Nucleotide sequences adjacent to these methyl-CpG sites were identified and we determined the methylation level by methylation-sensitive restriction endonuclease (MSRE)-PCR analysis to confirm the accuracy of AFLP analysis. The differences of the methylation level among tissues were almost identical among these methods. By MSD-AFLP analysis, we detected many CpGs showing less than 5% statistically significant tissue-specific difference and less than 10% degree of variability. Additionally, MSD-AFLP analysis could be used to identify CpG methylation sites in other organisms including humans.

Conclusion

MSD-AFLP analysis can potentially be used to measure slight changes in CpG methylation level. Regarding the remarkable precision, sensitivity, and throughput of MSD-AFLP analysis studies, this method will be advantageous in a variety of epigenetics-based research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12867-017-0083-2) contains supplementary material, which is available to authorized users.

Suppression Subtractive Hybridization Versus Next-Generation Sequencing in Plant Genetic Engineering: Challenges and Perspectives

Suppression subtractive hybridization (SSH) is an effective method to identify different genes with different expression levels involved in a variety of biological processes. This method has often been used to study molecular mechanisms of plants in complex relationships with different pathogens and a variety of biotic stresses. Compared to other techniques used in gene expression profiling, SSH needs relatively smaller amounts of the initial materials, with lower costs, and fewer false positives present within the results. Extraction of total RNA from plant species rich in phenolic compounds, carbohydrates, and polysaccharides that easily bind to nucleic acids through cellular mechanisms is difficult and needs to be considered. Remarkable advancement has been achieved in the next-generation sequencing (NGS) field. As a result of progress within fields related to molecular chemistry and biology as well as specialized engineering, parallelization in the sequencing reaction has exceptionally enhanced the overall read number of generated sequences per run. Currently available sequencing platforms support an earlier unparalleled view directly into complex mixes associated with RNA in addition to DNA samples. NGS technology has demonstrated the ability to sequence DNA with remarkable swiftness, therefore allowing previously unthinkable scientific accomplishments along with novel biological purposes. However, the massive amounts of data generated by NGS impose a substantial challenge with regard to data safe-keeping and analysis. This review examines some simple but vital points involved in preparing the initial material for SSH and introduces this method as well as its associated applications to detect different novel genes from different plant species. This review evaluates general concepts, basic applications, plus the probable results of NGS technology in genomics, with unique mention of feasible potential tools as well as bioinformatics.

Genetic Approaches to Study Plant Responses to Environmental Stresses: An Overview

The assessment of gene expression levels is an important step toward elucidating gene functions temporally and spatially. Decades ago, typical studies were focusing on a few genes individually, whereas now researchers are able to examine whole genomes at once. The upgrade of throughput levels aided the introduction of systems biology approaches whereby cell functional networks can be scrutinized in their entireties to unravel potential functional interacting components. The birth of systems biology goes hand-in-hand with huge technological advancements and enables a fairly rapid detection of all transcripts in studied biological samples. Even so, earlier technologies that were restricted to probing single genes or a subset of genes still have their place in research laboratories. The objective here is to highlight key approaches used in gene expression analysis in plant responses to environmental stresses, or, more generally, any other condition of interest. Northern blots, RNase protection assays, and qPCR are described for their targeted detection of one or a few transcripts at a once. Differential display and serial analysis of gene expression represent non-targeted methods to evaluate expression changes of a significant number of gene transcripts. Finally, microarrays and RNA-seq (next-generation sequencing) contribute to the ultimate goal of identifying and quantifying all transcripts in a cell under conditions or stages of study. Recent examples of applications as well as principles, advantages, and drawbacks of each method are contrasted. We also suggest replacing the term “Next-Generation Sequencing (NGS)” with another less confusing synonym such as “RNA-seq”, “high throughput sequencing”, or “massively parallel sequencing” to avoid confusion with any future sequencing technologies.

Large-scale tag/PCR-based gene expression profiling

An intriguing enigma in molecular biology is how genes within a single genome are differentially expressed in different cell types of a multicellular organism, or in response to different developmental or environmental queues in a single cell type. Quantification of transcript levels on a genome-wide scale, often termed transcript profiling, provides a powerful approach to identifying protein-coding and non-coding RNAs functionally relevant to a given biological process. Indeed, transcriptome analysis has been a key area of biological inquiry for decades and successfully produced discoveries in a multitude of processes and disease states, and in an increasingly large number of organisms. The evolution of technologies with increasing levels of informational content, ranging from hybridization-based technologies such as Northern blot analysis and microarrays to tag/polymerase chain reaction (PCR)- and sequence-based technologies including differential display and SAGE, along with the next-generation sequencing, has provided hope for revealing the molecular details of biological systems as they respond to change. This review is an overview of selected high throughput tag/PCR-based methods for genome-wide expression profiling amenable to high-throughput automated operation in any standard laboratory.

Transcriptomic and proteomic analyses of resistant host responses in Arachis diogoi challenged with late leaf spot pathogen, Phaeoisariopsis personata

Late leaf spot is a serious disease of peanut caused by the imperfect fungus, Phaeoisariopsis personata. Wild diploid species, Arachis diogoi. is reported to be highly resistant to this disease and asymptomatic. The objective of this study is to investigate the molecular responses of the wild peanut challenged with the late leaf spot pathogen using cDNA-AFLP and 2D proteomic study. A total of 233 reliable, differentially expressed genes were identified in Arachis diogoi. About one third of the TDFs exhibit no significant similarity with the known sequences in the data bases. Expressed sequence tag data showed that the characterized genes are involved in conferring resistance in the wild peanut to the pathogen challenge. Several genes for proteins involved in cell wall strengthening, hypersensitive cell death and resistance related proteins have been identified. Genes identified for other proteins appear to function in metabolism, signal transduction and defence. Nineteen TDFs based on the homology analysis of genes associated with defence, signal transduction and metabolism were further validated by quantitative real time PCR (qRT-PCR) analyses in resistant wild species in comparison with a susceptible peanut genotype in time course experiments. The proteins corresponding to six TDFs were differentially expressed at protein level also. Differentially expressed TDFs and proteins in wild peanut indicate its defence mechanism upon pathogen challenge and provide initial breakthrough of genes possibly involved in recognition events and early signalling responses to combat the pathogen through subsequent development of resistivity. This is the first attempt to elucidate the molecular basis of the response of the resistant genotype to the late leaf spot pathogen, and its defence mechanism.

Mutational breeding and genetic engineering in the development of high grain protein content

Cereals are the most important crops in the world for both human consumption and animal feed. Improving their nutritional values, such as high protein content, will have significant implications, from establishing healthy lifestyles to helping remediate malnutrition problems worldwide. Besides providing a source of carbohydrate, grain is also a natural source of dietary fiber, vitamins, minerals, specific oils, and other disease-fighting phytocompounds. Even though cereal grains contain relatively little protein compared to legume seeds, they provide protein for the nutrition of humans and livestock that is about 3 times that of legumes. Most cereal seeds lack a few essential amino acids; therefore, they have imbalanced amino acid profiles. Lysine (Lys), threonine (Thr), methionine (Met), and tryptophan (Trp) are among the most critical and are a limiting factor in many grain crops for human nutrition. Tremendous research has been put into the efforts to improve these essential amino acids. Development of high protein content can be outlined in four different approaches through manipulating seed protein bodies, modulating certain biosynthetic pathways to overproduce essential and limiting amino acids, increasing nitrogen relocation to the grain through the introduction of transgenes, and exploiting new genetic variance. Various technologies have been employed to improve protein content including conventional and mutational breeding, genetic engineering, marker-assisted selection, and genomic analysis. Each approach involves a combination of these technologies. Advancements in nutrigenomics and nutrigenetics continue to improve public knowledge at a rapid pace on the importance of specific aspects of food nutrition for optimum fitness and health. An understanding of the molecular basis for human health and genetic predisposition to certain diseases through human genomes enables individuals to personalize their nutritional requirements. It is critically important, therefore, to improve grain protein quality. Highly nutritious grain can be tailored to functional foods to meet the needs for both specific individuals and human populations as a whole.

Differential gene expression in response to Papaya ringspot virus infection in Cucumis metuliferus using cDNA-amplified fragment length polymorphism analysis

A better understanding of virus resistance mechanisms can offer more effective strategies to control virus diseases. Papaya ringspot virus (PRSV), Potyviridae, causes severe economical losses in papaya and cucurbit production worldwide. However, no resistance gene against PRSV has been identified to date. This study aimed to identify candidate PRSV resistance genes using cDNA-AFLP analysis and offered an open architecture and transcriptomic method to study those transcripts differentially expressed after virus inoculation. The whole genome expression profile of Cucumis metuliferus inoculated with PRSV was generated using cDNA-amplified fragment length polymorphism (cDNA-AFLP) method. Transcript derived fragments (TDFs) identified from the resistant line PI 292190 may represent genes involved in the mechanism of PRSV resistance. C. metuliferus susceptible Acc. 2459 and resistant PI 292190 lines were inoculated with PRSV and subsequently total RNA was isolated for cDNA-AFLP analysis. More than 400 TDFs were expressed specifically in resistant line PI 292190. A total of 116 TDFs were cloned and their expression patterns and putative functions in the PRSV-resistance mechanism were further characterized. Subsequently, 28 out of 116 candidates which showed two-fold higher expression levels in resistant PI 292190 than those in susceptible Acc. 2459 after virus inoculation were selected from the reverse northern blot and bioinformatic analysis. Furthermore, the time point expression profiles of these candidates by northern blot analysis suggested that they might play roles in resistance against PRSV and could potentially provide valuable information for controlling PRSV disease in the future.

Investigation of transcriptional responses of juvenile mouse bone marrow to power frequency magnetic fields

To seek alterations in gene transcription in bone marrow cells following in vivo exposure of juvenile mice to power frequency magnetic fields, young (21-24-day old) C57BL/6 mice were exposed to a 100μT 50Hz magnetic field for 2h. Transcription was analysed by three methods, High Coverage Expression Profiling (HiCEP), Illumina microarrays and quantitative real-time polymerase chain reaction (QRT-PCR). A pilot HiCEP experiment with 6 exposed (E) and 6 non-exposed (NE) mice identified four candidate responsive transcripts (two unknown transcripts (AK152075 and F10-NED), phosphatidylinositol binding clathrin assembly protein (Picalm) and exportin 7 (Xpo7)). A larger experiment compared 19 E and 15 NE mice using two independent QRT-PCR assays and repeated microarray assays. No significant field-dependent changes were seen, although Picalm showed a trend to significance in one QRT-PCR assay (E/NE=0.91; P=0.06). However, the study was underpowered to detect an effect of this magnitude (52% power at P=0.05). These data indicate the current experimental constraints in detecting small changes in transcription that may occur in response to magnetic fields. These constraints result from technical limitations in the accuracy of assays and biological variation, which together were sufficient to account statistically for the number of differentially expressed transcripts identified in the pilot experiment.

Identification of differential expression genes in leaves of rice (Oryza sativa L.) in response to heat stress by cDNA-AFLP analysis

High temperature impedes the growth and productivity of various crop species. To date, rice (Oryza sativa L.) has not been exploited to understand the molecular basis of its abnormally high level of temperature tolerance. To identify transcripts induced by heat stress, twenty-day-old rice seedlings of different rice cultivars suffering from heat stress were treated at different times, and differential gene expression analyses in leaves were performed by cDNA-AFLP and further verified by real-time RT-PCR. In aggregate, more than three thousand different fragments were indentified, and 49 fragments were selected for the sequence and differential expressed genes were classified functionally into different groups. 6 of 49 fragments were measured by real-time RT-PCR. In addition, the variations of three different polyamine contents in response to heat stress through high-performance liquid chromatography (HPLC) analysis were also performed. The results and their direct and indirect relationships to heat stress tolerance mechanism were discussed.

Chicken hemogen homolog is involved in the chicken-specific sex-determining mechanism

Using a comprehensive transcriptome analysis, a Z chromosome-linked chicken homolog of hemogen (cHEMGN) was identified and shown to be specifically involved in testis differentiation in early chicken embryos. Hemogen [Hemgn in mice, EDAG (erythroid differentiation-associated gene protein) in humans] was recently characterized as a hematopoietic tissue-specific gene encoding a transcription factor that regulates the proliferation and differentiation of hematopoietic cells in mammals. In chicken, cHEMGN was expressed not only in hematopoietic tissues but also in the early embryonic gonad of male chickens. The male-specific expression was identified in the nucleus of (pre)Sertoli cells after the sex determination period and before the expression of SOX9 (SRY-box 9). The expression of cHEMGN was induced in ZW embryonic gonads that were masculinized by aromatase inhibitor treatment. ZW embryos overexpressing cHEMGN, generated by infection with retrovirus carrying cHEMGN, showed masculinized gonads. These findings suggest that cHEMGN is a transcription factor specifically involved in chicken sex determination.

Transcriptional profiling by cDNA-AFLP analysis showed differential transcript abundance in response to water stress in Populus hopeiensis

BACKGROUND

Drought is one of the main environmental factors limiting tree growth and productivity of plantation forests worldwide. Populus hopeiensis Hu et Chow is one of the most important commercial plantation tree species in China. However, the genes controlling drought tolerance in this species have not been identified or characterized. Here, we conducted differential expression analyses and identified a number of genes that were up- or downregulated in P. hopeiensis during water stress. To the best of our knowledge, this is the first comprehensive study of differentially expressed genes in water-stressed P. hopeiensis.

RESULTS

Using the cDNA-AFLP detection technique, we used 256 primer combinations to identify differentially expressed genes in P. hopeiensis during water stress. In total, 415 transcript derived-fragments (TDFs) were obtained from 10× deep sequencing of 473 selected TDFs. Of the 415 TDFs, 412 were annotated by BLAST searches against various databases. The majority of these genes encoded products involved in ion transport and compartmentalization, cell division, metabolism, and protein synthesis. The TDFs were clustered into 12 groups on the basis of their expression patterns. Of the 415 reliable TDFs, the sequences of 35 were homologous to genes that play roles in short or long-term resistance to drought stress. Some genes were further selected for validation of cDNA-AFLP expression patterns using real-time PCR analyses. The results confirmed the expression patterns that were detected using the cDNA-AFLP technique.

CONCLUSION

The cDNA-AFLP technique is an effective and powerful tool for identifying candidate genes that are differentially expressed under water stress. We demonstrated that 415 TDFs were differentially expressed in water-stressed poplar. The products of these genes are involved in various biological processes in the drought response of poplar. The results of this study will aid in the identification of candidate genes of future experiments aimed at understanding this response of poplar.

Mixture of differentially tagged Tol2 transposons accelerates conditional disruption of a broad spectrum of genes in mouse embryonic stem cells

Among the insertional mutagenesis techniques used in the current international knockout mouse project (KOMP) on the inactivation of all mouse genes in embryonic stem (ES) cells, random gene trapping has been playing a major role. Gene-targeting experiments have also been performed to individually and conditionally knockout the remaining ‘difficult-to-trap’ genes. Here, we show that transcriptionally silent genes in ES cells are severely underrepresented among the randomly trapped genes in KOMP. Our conditional poly(A)-trapping vector with a common retroviral backbone also has a strong bias to be integrated into constitutively transcribed genome loci. Most importantly, conditional gene disruption could not be successfully accomplished by using the retrovirus vector because of the frequent development of intra-vector deletions/rearrangements. We found that one of the cut and paste-type DNA transposons, Tol2, can serve as an ideal platform for gene-trap vectors that ensures identification and conditional disruption of a broad spectrum of genes in ES cells. We also solved a long-standing problem associated with multiple vector integration into the genome of a single cell by incorporating a mixture of differentially tagged Tol2 transposons. We believe our strategy indicates a straightforward approach to mass-production of conditionally disrupted alleles for genes in the target cells.

Crucial role of c-Myc in the generation of induced pluripotent stem cells

c-Myc transduction has been considered previously to be nonessential for induced pluripotent stem cell (iPSC) generation. In this study, we investigated the effects of c-Myc transduction on the generation of iPSCs from an inbred mouse strain using a genome integration-free vector to exclude the effects of the genetic background and the genomic integration of exogenous genes. Our findings reveal a clear difference between iPSCs generated using the four defined factors including c-Myc (4F-iPSCs) and those produced without c-Myc (3F-iPSCs). Molecular and cellular analyses did not reveal any differences between 3F-iPSCs and 4F-iPSCs, as reported previously. However, a chimeric mice formation test indicated clear differences, whereby few highly chimeric mice and no germline transmission was observed using 3F-iPSCs. Similar differences were also observed in the mouse line that has been widely used in iPSC studies. Furthermore, the defect in 3F-iPSCs was considerably improved by trichostatin A, a histone deacetyl transferase inhibitor, indicating that c-Myc plays a crucial role in iPSC generation through the control of histone acetylation. Indeed, low levels of histone acetylation were observed in 3F-iPSCs. Our results shed new light on iPSC generation mechanisms and strongly recommend c-Myc transduction for preparing high-quality iPSCs.

Robust ordered mRNA differential display: an improved method for global gene expression profiling

Global gene expression profiling (GGEP) plays a pivotal role in biological research. We developed an improved GGEP method called "robust ordered mRNA differential display (RoDD)" by combining mRNA differential display (DD) and complementary DNA amplified fragment length polymorphisms (cDNA-AFLP) using elaborately designed primers and a poly (dT:A) replacement technique. Redundancy was minimized by bead-based isolation and coverage was improved by using restriction enzymes that recognized 4-bp sites. This method offers the common virtues of gel-based methods along with the reliability of cDNA-AFLP. The most significant advantage of RoDD over current gel-based methods is greatly improved coverage and minimized redundancy.

Gene regulation by cytokinin in Arabidopsis

The plant hormone cytokinin realizes at least part of its signaling output through the regulation of gene expression. A great part of the early transcriptional regulation is mediated by type-B response regulators, which are transcription factors of the MYB family. Other transcription factors, such as the cytokinin response factors of the AP2/ERF family, have also been shown to be involved in this process. Additional transcription factors mediate distinct parts of the cytokinin response through tissue- and cell-specific downstream transcriptional cascades. In Arabidopsis, only a single cytokinin response element, to which type-B response regulators bind, has been clearly proven so far, which has 5'-GAT(T/C)-3' as a core sequence. This motif has served to construct a synthetic cytokinin-sensitive two-component system response element, which is useful for monitoring the cellular cytokinin status. Insight into the extent of transcriptional regulation has been gained by genome-wide gene expression analyses following cytokinin treatment and from plants having an altered cytokinin content or signaling. This review presents a meta analysis of such microarray data resulting in a core list of cytokinin response genes. Genes encoding type-A response regulators displayed the most stable response to cytokinin, but a number of cytokinin metabolism genes (CKX4, CKX5, CYP735A2, UGT76C2) also belong to them, indicating homeostatic mechanisms operating at the transcriptional level. The cytokinin core response genes are also the target of other hormones as well as biotic and abiotic stresses, documenting crosstalk of the cytokinin system with other hormonal and environmental signaling pathways. The multiple links of cytokinin to diverse functions, ranging from control of meristem activity, hormonal crosstalk, nutrient acquisition, and various stress responses, are also corroborated by a compilation of genes that have been repeatedly found by independent gene expression profiling studies. Such functions are, at least in part, supported by genetic studies.

ERK signaling controls blastema cell differentiation during planarian regeneration

The robust regenerative ability of planarians depends on a population of somatic stem cells called neoblasts, which are the only mitotic cells in adults and are responsible for blastema formation after amputation. The molecular mechanism underlying neoblast differentiation associated with blastema formation remains unknown. Here, using the planarian Dugesia japonica we found that DjmkpA, a planarian mitogen-activated protein kinase (MAPK) phosphatase-related gene, was specifically expressed in blastema cells in response to increased extracellular signal-related kinase (ERK) activity. Pharmacological and genetic [RNA interference (RNAi)] approaches provided evidence that ERK activity was required for blastema cells to exit the proliferative state and undergo differentiation. By contrast, DjmkpA RNAi induced an increased level of ERK activity and rescued the differentiation defect of blastema cells caused by pharmacological reduction of ERK activity. These observations suggest that ERK signaling plays an instructive role in the cell fate decisions of blastema cells regarding whether to differentiate or not, by inducing DjmkpA as a negative regulator of ERK signaling during planarian regeneration.

Identification of imprinted genes subject to parent-of-origin specific expression in Arabidopsis thaliana seeds

BACKGROUND

Epigenetic regulation of gene dosage by genomic imprinting of some autosomal genes facilitates normal reproductive development in both mammals and flowering plants. While many imprinted genes have been identified and intensively studied in mammals, smaller numbers have been characterized in flowering plants, mostly in Arabidopsis thaliana. Identification of additional imprinted loci in flowering plants by genome-wide screening for parent-of-origin specific uniparental expression in seed tissues will facilitate our understanding of the origins and functions of imprinted genes in flowering plants.

RESULTS

cDNA-AFLP can detect allele-specific expression that is parent-of-origin dependent for expressed genes in which restriction site polymorphisms exist in the transcripts derived from each allele. Using a genome-wide cDNA-AFLP screen surveying allele-specific expression of 4500 transcript-derived fragments, we report the identification of 52 maternally expressed genes (MEGs) displaying parent-of-origin dependent expression patterns in Arabidopsis siliques containing F1 hybrid seeds (3, 4 and 5 days after pollination). We identified these MEGs by developing a bioinformatics tool (GenFrag) which can directly determine the identities of transcript-derived fragments from (i) their size and (ii) which selective nucleotides were added to the primers used to generate them. Hence, GenFrag facilitates increased throughput for genome-wide cDNA-AFLP fragment analyses. The 52 MEGs we identified were further filtered for high expression levels in the endosperm relative to the seed coat to identify the candidate genes most likely representing novel imprinted genes expressed in the endosperm of Arabidopsis thaliana. Expression in seed tissues of the three top-ranked candidate genes, ATCDC48, PDE120 and MS5-like, was confirmed by Laser-Capture Microdissection and qRT-PCR analysis. Maternal-specific expression of these genes in Arabidopsis thaliana F1 seeds was confirmed via allele-specific transcript analysis across a range of different accessions. Differentially methylated regions were identified adjacent to ATCDC48 and PDE120, which may represent candidate imprinting control regions. Finally, we demonstrate that expression levels of these three genes in vegetative tissues are MET1-dependent, while their uniparental maternal expression in the seed is not dependent on MET1.

CONCLUSIONS

Using a cDNA-AFLP transcriptome profiling approach, we have identified three genes, ATCDC48, PDE120 and MS5-like which represent novel maternally expressed imprinted genes in the Arabidopsis thaliana seed. The extent of overlap between our cDNA-AFLP screen for maternally expressed imprinted genes, and other screens for imprinted and endosperm-expressed genes is discussed.

Environmental protection: researches in National Institute of Radiological Sciences

Some studies for radiological protection of the environment have been made at the National Institute of Radiological Sciences (NIRS). Transfer of radionuclides and related elements has been investigated for dose estimation of non-human biota. A parameter database and radionuclide transfer models have been also developed for the Japanese environments. Dose (rate)-effect relationships for survival, growth and reproduction have been investigated in conifers, Arabidopsis, fungi, earthworms, springtails, algae, duckweeds, daphnia and medaka. Also genome-wide gene expression analysis has been carried out by high coverage expression profiling (HiCEP). Effects on aquatic microbial communities have been studied in experimental ecosystem models, i.e., microcosms. Some effects were detected at a dose rate of 1 Gy day(-1) and were likely to arise from interspecies interactions. The results obtained at NIRS have been used in development of frameworks for environmental protection by some international bodies, and will contribute to environmental protection in Japan and other Asian countries.

Profiling differential gene expression of symbiotic and aposymbiotic corals using a high coverage gene expression profiling (HiCEP) analysis

Coral generally harbors zooxanthellae (genus Symbiodinium) in the body for mutualistic symbiosis, which favors the host through effects on growth, stress response, and nutrient utilization. However, little is known about the molecular mechanisms by which the partners establish and regulate the endosymbiosis. In this study, we conducted a comprehensive transcriptome analysis in the coral Acropora tenuis using a high coverage gene expression profiling (HiCEP) method, to assess the genes that are involved in the coral-zooxanthellae symbiosis. For this purpose, we compared between aposymbiotic juveniles and those inoculated with a cultured monoclonal Symbiodinium species in two different clades (PL-TS-1 or CCMP2467). Among the 765 genes that exhibited different expression profiles between the two groups, 462 were upregulated and 303 downregulated by the symbiosis with somewhat variable responses to the two different symbionts. Among the responsive genes, we could annotate 33 genes by bioinformatic analyses and confirmed that their expression is actually altered in the same direction in the symbiotic individuals using real-time polymerase chain reaction. Functional analyses of the annotated genes indicate that they are involved in carbohydrate and lipid metabolism, intracellular signal transduction, and membrane transport of ions in the host corals as expected from the endosymbiosis of zooxanthellae.

mRNA expression of a cadmium-responsive gene is a sensitive biomarker of cadmium exposure in the soil collembolan Folsomia candida

The gene expression of environmental organisms is useful as a biomarker of environmental pollution. One of its advantages is high sensitivity. We identified the cDNA of a novel cadmium-responsive gene in the soil collembolan Folsomia candida. The deduced protein, designated "metallothionein-like motif containing protein" (MTC), was cysteine-rich and contained a metallothionein-like motif with similarity to metallothionein, but had a much longer sequence than metallothionein and contained repeated sequences of amino acids. Expression of MTC mRNA was sensitively induced by cadmium exposure at 0.3 mg/kg of dry food, a concentration at which toxic effects are not observed, but expression was not affected by gamma-ray exposure (an inducer of oxidative stress). These findings suggest that MTC is involved in cadmium-binding processes rather than in oxidative-stress responses. In conclusion, we suggest that gene expression of MTC may be a candidate biomarker for detecting low levels of cadmium contamination in soil.

Identification of genes involved in stem rust resistance from wheat mutant D51 with the cDNA-AFLP technique

Wheat (Triticum aestivum L.) stem rust caused by Puccinia graminis f. sp. tritici is one of the main diseases of wheat worldwide. Wheat mutant line D51, which was derived from the highly susceptible cultivar L6239, shows resistance to the prevailing races 21C3CPH, 21C3CKH, and 21C3CTR of P. graminis f. sp. tritici in China. In this study, we used the cDNA-AFLP technology to identify the genes that are likely involved in the stem rust resistance. EcoRI/MseI selective primers were used to generate approximately 1920 DNA fragments. Seventy five differentially transcribed fragments (3.91%) were identified by comparing the samples of 21C3CPH infected D51 with infected L6239 or uninfected D51. Eleven amplified cDNA fragments were sequenced. Eight showed significant similarity to known genes, including TaLr1 (leaf rust resistance gene), wlm24 (wheat powdery mildew resistance gene), stress response genes and ESTs of environment stress of tall fescue. These identified genes are involved in plant defense response and stem rust resistance and need further research to determine their usefulness in breeding new resistance cultivars.

An improved and cost-effective cDNA-AFLP method to investigate transcription-derived products when high throughput sequencing is not available

Although the cost of high throughput sequencing is decreasing, the cost is still often too high for individual projects targeted at, e.g., genome-wide transcription profiling in non-model organisms. Then, a low-cost alternative is cDNA-AFLP, which we have now considerably modified in order to develop a faster and simpler method to identify and analyze genes involved in specific, possibly adaptive characteristics. Particularly, we wanted to exclude repetitive PCR amplifications, extensive cloning and the presence of overlapping transcripts, which all lower the efficiency of the method.

Recent advances in the biology of heavy-ion cancer therapy

Superb biological effectiveness and dose conformity represent a rationale for heavy-ion therapy, which has thus far achieved good cancer controllability while sparing critical normal organs. Immediately after irradiation, heavy ions produce dense ionization along their trajectories, cause irreparable clustered DNA damage, and alter cellular ultrastructure. These ions, as a consequence, inactivate cells more effectively with less cell-cycle and oxygen dependence than conventional photons. The modes of heavy ion-induced cell death/inactivation include apoptosis, necrosis, autophagy, premature senescence, accelerated differentiation, delayed reproductive death of progeny cells, and bystander cell death. This paper briefly reviews the current knowledge of the biological aspects of heavy-ion therapy, with emphasis on the authors' recent findings. The topics include (i) repair mechanisms of heavy ion-induced DNA damage, (ii) superior effects of heavy ions on radioresistant tumor cells (intratumor quiescent cell population, TP53-mutated and BCL2-overexpressing tumors), (iii) novel capacity of heavy ions in suppressing cancer metastasis and neoangiogenesis, and (iv) potential of heavy ions to induce secondary (especially breast) cancer.

Tourette syndrome and klippel-feil anomaly in a child with chromosome 22q11 duplication

This is the first case description of the association of Klippel-Feil Syndrome (KFS), Tourette Syndrome (TS), Motor Stereotypies, and Obsessive Compulsive Behavior, with chromosome 22q11.2 Duplication Syndrome (22q11DupS). Neuropsychiatric symptoms in persons with 22q11.2 deletion, including obsessive compulsiveness, anxiety, hyperactivity, and one prior case report of TS, have been attributed to low copy number effects on Catechol-O-Methyltransferase (COMT). However, the present unique case of 22q11DupS and TS suggests a more complex relationship, either for low- or high-COMT activity, or for other genes at this locus.

A comprehensive analysis of the soybean genes and proteins expressed under flooding stress using transcriptome and proteome techniques

The inducible genes and proteins were analyzed using transcriptome and proteome techniques to explore the mechanisms underlying soybean response to flooding stress. Soybean seedlings were germinated for 2 days and subjected to flooding for 12 h, and the total RNAs and proteins were extracted from the root and hypocotyl. High-coverage gene expression profiling analysis as transcriptome technique was performed. Ninety-seven out of the 29,388 peaks observed demonstrated a greater than 25-fold change following 12 h of flood-induced stress. Furthermore, 34 proteins out of 799 proteins were changed by 12 h stress. Genes associated with alcohol fermentation, ethylene biosynthesis, pathogen defense, and cell wall loosening were significantly up-regulated. Hemoglobin, acid phosphatase, and Kunitz trypsin protease inhibitor were altered at both transcriptional and translational levels. Reactive oxygen species scavengers and chaperons were changed only at the translational level. It is suggested that the early response of soybean under flooding might be important stress adaptation to ensure survival against not only hypoxia but also the direct damage of cell by water.

Genome-wide gene expression profiling in gamma-irradiated green alga, Pseudokirchneriella subcapitata, by HiCEP

Transcriptome was analyzed in gamma-irradiated green alga Pseudokirchneriella subcapitata by high coverage gene expression profiling (HiCEP). Approximately 7,800 expressed genes were detected. Expression levels of 623-707 genes were affected at 100-300 Gy. Nucleotide sequences of 41 up-regulated genes were determined. The quantitative reverse transcription polymerase chain reaction validated the up-regulation. Two genes had homology to genes related to ionizing radiation. These results indicate usefulness of HiCEP for screening of stress-responsive genes in species that are ecotoxicologically important but for which genomic sequence information is lacking.

Soybean proteome database: a data resource for plant differential omics

The Soybean Proteome Database aims to be a data repository for functional analyses of soybean responses to flooding injury, recognized as a major constraint for establishment and production of this plant. The current release contains 21 reference maps of soybean (Glycine max cv. Enrei) proteins electrophoresed on two-dimensional polyacrylamide gels of which the samples were collected from several organs, tissues and organelles. These reference maps include 7311 detected proteins and 532 identified proteins, or proteins for which a sequence or peptide peak has been determined. The database is searchable by protein properties such as accession number, description and isoelectric point and molecular weight range. The Soybean Proteome Database also integrates multiple "omes". An omics table reveals relationships among 106 mRNAs, 51 proteins and 89 metabolites that vary over time under flooding stress. The tabulated metabolites are anchored to a metabolome network. A unified temporal-profile tag attached to the mRNAs, proteins and metabolites facilitates retrieval of the data based on the temporal expression profiles. A graphical user interface based on dynamic HTML facilitates viewing the metabolome network as well as the profiles of multiple omes in a uniform manner. The entire database is available at http://proteome.dc.affrc.go.jp/Soybean/.

cDNA-AFLP analysis reveals differential gene expression in response to salt stress in foxtail millet (Setaria italica L.)

Plant growth and productivity are affected by various abiotic stresses such as heat, drought, cold, salinity, etc. The mechanism of salt tolerance is one of the most important subjects in plant science as salt stress decreases worldwide agricultural production. In our present study we used cDNA-AFLP technique to compare gene expression profiles of a salt tolerant and a salt-sensitive cultivar of foxtail millet (Seteria italica) in response to salt stress to identify early responsive differentially expressed transcripts accumulated upon salt stress and validate the obtained result through quantitative real-time PCR (qRT-PCR). The expression profile was compared between a salt tolerant (Prasad) and susceptible variety (Lepakshi) of foxtail millet in both control condition (L0 and P0) and after 1 h (L1 and P1) of salt stress. We identified 90 transcript-derived fragments (TDFs) that are differentially expressed, out of which 86 TDFs were classified on the basis of their either complete presence or absence (qualitative variants) and 4 on differential expression pattern levels (quantitative variants) in the two varieties. Finally, we identified 27 non-redundant differentially expressed cDNAs that are unique to salt tolerant variety which represent different groups of genes involved in metabolism, cellular transport, cell signaling, transcriptional regulation, mRNA splicing, seed development and storage, etc. The expression patterns of seven out of nine such genes showed a significant increase of differential expression in tolerant variety after 1 h of salt stress in comparison to salt-sensitive variety as analyzed by qRT-PCR. The direct and indirect relationship of identified TDFs with salinity tolerance mechanism is discussed.

Improved coverage of cDNA-AFLP by sequential digestion of immobilized cDNA

BACKGROUND

cDNA-AFLP is a transcriptomics technique which does not require prior sequence information and can therefore be used as a gene discovery tool. The method is based on selective amplification of cDNA fragments generated by restriction endonucleases, electrophoretic separation of the products and comparison of the band patterns between treated samples and controls. Unequal distribution of restriction sites used to generate cDNA fragments negatively affects the performance of cDNA-AFLP. Some transcripts are represented by more than one fragment while other escape detection, causing redundancy and reducing the coverage of the analysis, respectively.

RESULTS

With the goal of improving the coverage of cDNA-AFLP without increasing its redundancy, we designed a modified cDNA-AFLP protocol. Immobilized cDNA is sequentially digested with several restriction endonucleases and the released DNA fragments are collected in mutually exclusive pools. To investigate the performance of the protocol, software tool MECS (Multiple Enzyme cDNA-AFLP Simulation) was written in Perl. cDNA-AFLP protocols described in the literature and the new sequential digestion protocol were simulated on sets of cDNA sequences from mouse, human and Arabidopsis thaliana. The redundancy and coverage, the total number of PCR reactions, and the average fragment length were calculated for each protocol and cDNA set.

CONCLUSION

Simulation revealed that sequential digestion of immobilized cDNA followed by the partitioning of released fragments into mutually exclusive pools outperformed other cDNA-AFLP protocols in terms of coverage, redundancy, fragment length, and the total number of PCRs. Primers generating 30 to 70 amplicons per PCR provided the highest fraction of electrophoretically distinguishable fragments suitable for normalization. For A. thaliana, human and mice transcriptome, the use of two marking enzymes and three sequentially applied releasing enzymes for each of the marking enzymes is recommended.

Application of HiCEP to screening of radiation stress-responsive genes in the soil microarthropod Folsomia candida (Collembola)

The field of ecotoxicogenomics has received increasing attention for its potential to provide insight into pressing ecological issues. However, its applications are limited due to a lack of genetic sequence information for organisms used in ecotoxicological studies. We used high-coverage expression profiling (HiCEP), a method that requires no prior sequence knowledge, to examine stress-responsive genes and their dose dependence in the springtail Folsomia candida using gamma radiation as the stressor. Radiation-responsive genes and their dose dependency were detected at effective doses for reproduction, and 16 up-regulated transcript-derived fragments (TDFs) were sequenced. Quantitative PCR analysis also found that most of the TDFs were up-regulated. The sequences of the TDFs showed resemblance to known genes, such as glutathione S-transferase and poly(ADP-ribose) polymerase, but most showed no similarity to any genes in the gene databases. These results suggest that HiCEP is effective for discovering differently expressed genes and their dose dependence, even in organisms for which few sequence data are available. The limited length of the TDFs, however, may impede functional annotation of the genes. In conclusion, HiCEP is useful for ecotoxicogenomic studies in which various organisms with few available genomic resources are involved.

Gibberellin-regulated gene in the basal region of rice leaf sheath encodes basic helix-loop-helix transcription factor

Genes regulated by gibberellin (GA) during leaf sheath elongation in rice seedlings were identified using the transcriptome approach. mRNA from the basal regions of leaf sheaths treated with GA3 was analyzed by high-coverage gene expression profiling. 33,004 peaks were detected, and 30 transcripts showed significant changes in the presence of GA3. Among these, basic helix-loop-helix transcription factor (AK073385) was significantly upregulated. Quantitative PCR analysis confirmed that expression of AK073385 was controlled by GA3 in a time- and dose-dependent manner. Basic helix-loop-helix transcription factor (AK073385) is therefore involved in the regulation of gene expression by GA3.

Radiation induction of delayed recombination in Schizosaccharomyces pombe

Ionizing radiation is known to induce delayed chromosome and gene mutations in the descendants of the irradiated tissue culture cells. Molecular mechanisms of such delayed mutations are yet to be elucidated, since high genomic complexity of mammalian cells makes it difficult to analyze. We now tested radiation induction of delayed recombination in the fission yeast Schizosaccharomyces pombe by monitoring the frequency of homologous recombination after X-irradiation. A reporter with 200 bp tandem repeats went through spontaneous recombination at a frequency of 1.0 x 10(-4), and the frequency increased dose-dependently to around 10 x 10(-4) at 500 Gy of X-irradiation. Although the repair of initial DNA damage was thought to be completed before the restart of cell division cycle, the elevation of the recombination frequency persisted for 8-10 cell generations after irradiation (delayed recombination). The delayed recombination suggests that descendants of the irradiated cells keep a memory of the initial DNA damage which upregulates recombination machinery for 8-10 generations even in the absence of DNA double-strand breaks (DSBs). Since radical scavengers were ineffective in inhibiting the delayed recombination, a memory by continuous production of DNA damaging agents such as reactive oxygen species (ROS) was excluded. Recombination was induced in trans in a reporter on chromosome III by a DNA DSB at a site on chromosome I, suggesting the untargeted nature of delayed recombination. Interestingly, Rad22 foci persisted in the X-irradiated population in parallel with the elevation of the recombination frequency. These results suggest that the epigenetic damage memory induced by DNA DSB upregulates untargeted and delayed recombination in S. pombe.

[DNA methylation in cotton hybrids and their parents]

The possible role of methylation in the performance of heterosis has been analyzed in many crops. To further study this possibility, we investigated both the differences in cytosine methylation patterns between cotton heterotic hybrid/nonheterotic hybrids and their parental lines and the change in methylation level from seedling stage to flowering stage by using the methylation-sensitive amplified polymorphism (MSAP) method. The results showed that the number of demethylation loci in highly heterotic hybrids was greater that in lowly heterotic hybrids, and the level of DNA cytosine methylation in cotton at the seedling stage is higher than that at the flowering stage. The altered methylation patterns at low-copy genomic regions can be confirmed by DNA gel blot analysis. A total of 39 fragments that showed different methylation patterns were cloned and sequenced. The methylation status of these genes was modified differentially in hybrid and parents, suggesting that these genes might play a role in the performance of heterosis.

Identification of Genes Responding to Low-Dose Arsenite Using HiCEP

ABSTRACT Chronic ingestion of arsenic in polluted food and water can cause various human disorders including skin and lung cancers. Sensitive biomarkers from human tissue/cells could help to prevent chronic intoxication with low-dose arsenite. Using High-Coverage Expression Profiling (HiCEP), an Amplified Fragment Length Polymorphism (AFLP)-based gene expression profiling technique, we analyzed the expression of approximately 11,000 genes in human lung fibroblasts (HFLIII) and compared the profiles between cells, treated and untreated with 1 muM sodium arsenite (NaAsO(2)). Hundreds of genes appeared upregulated and downregulated more than two-fold, 2 h after the treatment. Marked induction was found (>4.4-fold) in a few genes including HMOX1, INHBA, and ANKRD11. Induction of the HMOX1 was detected with a dose of arsenite at as low as 0.3 muM (0.04 ppm) and reached its maximum at 4 h after the treatment. The arsenite-induced HMOX1 expression was attenuated by the promoted glutathione (GSH) synthesis by N-acetyl-L-cysteine (NAC). However, it was not affected by pretreating the cells with general radical scavengers, consistent with the fact that ionizing radiation at either high-or low-doses has never induced HMOX1 in the same assay system. Thus, induction of HMOX1 gene is highly sensitive and also selective against arsenite in the cells. The present process could provide a useful strategy for exploring biomarkers that might help in assessing the known and unknown risks of any natural and artificial toxic reagents.

GOGOT: a method for the identification of differentially expressed fragments from cDNA-AFLP data

Background

One-dimensional (1-D) electrophoretic data obtained using the cDNA-AFLP method have attracted great interest for the identification of differentially expressed transcript-derived fragments (TDFs). However, high-throughput analysis of the cDNA-AFLP data is currently limited by the need for labor-intensive visual evaluation of multiple electropherograms. We would like to have high-throughput ways of identifying such TDFs.

Results

We describe a method, GOGOT, which automatically detects the differentially expressed TDFs in a set of time-course electropherograms. Analysis by GOGOT is conducted as follows: correction of fragment lengths of TDFs, alignment of identical TDFs across different electropherograms, normalization of peak heights, and identification of differentially expressed TDFs using a special statistic. The output of the analysis is a highly reduced list of differentially expressed TDFs. Visual evaluation confirmed that the peak alignment was performed perfectly for the TDFs by virtue of the correction of peak fragment lengths before alignment in step 1. The validity of the automated ranking of TDFs by the special statistic was confirmed by the visual evaluation of a third party.

Conclusion

GOGOT is useful for the automated detection of differentially expressed TDFs from cDNA-AFLP temporal electrophoretic data. The current algorithm may be applied to other electrophoretic data and temporal microarray data.

ARR1 directly activates cytokinin response genes that encode proteins with diverse regulatory functions

Plant cells respond to cytokinins by changing their gene expression patterns. The histidyl-aspartyl (His-Asp) phosphorelay mediates the signal from cytokinin receptors to type-B response regulators including ARR1, which transactivate cytokinin primary response genes. However, the overall architecture of the signal cascade leading to cytokinin-responsive phenomena is still unclear, mainly because it is not known how the His-Asp phosphorelay is connected to downstream phenomena. To reveal events immediately downstream from the phosphorelay-mediated transcriptional activation, we searched for direct-target genes of ARR1 by exploiting ARR1DeltaDDK-GR, a chimeric transcription factor that transactivates ARR1 direct-target genes in transgenic plants by glucocorticoid induction. We identified 23 direct-target genes, most of which were found to be cytokinin primary response genes. The arr1-1 mutation clearly affected the primary response in at least 17 genes, meaning that they respond primarily to cytokinins through the function of ARR1. The 17 genes encode proteins with diverse functions, including type-A response regulators, cytokinin metabolic enzymes and putative disease resistance response proteins. These results provide novel evidence indicating that the His-Asp phosphorelay is connected to diverse regulatory levels of cytokinin-responsive phenomena through ARR1 direct-target genes.