Sensitivity of 70-mer oligonucleotides and cDNAs for microarray analysis of gene expression in Arabidopsis and its related species

Brassica villosa, a system for studying non-glandular trichomes and genes in the Brassicas

Brassica villosa is a wild Brassica C genome species with very dense trichome coverage and strong resistance to many insect pests of Brassica oilseeds and vegetables. Transcriptome analysis of hairy B. villosa leaves indicated higher expression of several important trichome initiation genes compared with glabrous B. napus leaves and consistent with the Arabidopsis model of trichome development. However, transcripts of the TRY inhibitory gene in hairy B. villosa were surprisingly high relative to B. napus and relative transcript levels of SAD2, EGL3, and several XIX genes were low, suggesting potential ancillary or less important trichome-related roles for these genes in Brassica species compared with Arabidopsis. Several antioxidant, calcium, non-calcium metal and secondary metabolite genes also showed differential expression between these two species. These coincided with accumulation of two alkaloid-like compounds, high levels of calcium, and other metals in B. villosa trichomes that are correlated with the known tolerance of B. villosa to high salt and the calcium-rich natural habitat of this wild species. This first time report on the isolation of large amounts of pure B. villosa trichomes, on trichome content, and on relative gene expression differences in an exceptionally hairy Brassica species compared with a glabrous species opens doors for the scientific community to understand trichome gene function in the Brassicas and highlights the potential of B. villosa as a trichome research platform.

Integrative biological analysis for neuropsychopharmacology

Although advances in psychotherapy have been made in recent years, drug discovery for brain diseases such as schizophrenia and mood disorders has stagnated. The need for new biomarkers and validated therapeutic targets in the field of neuropsychopharmacology is widely unmet. The brain is the most complex part of human anatomy from the standpoint of number and types of cells, their interconnections, and circuitry. To better meet patient needs, improved methods to approach brain studies by understanding functional networks that interact with the genome are being developed. The integrated biological approaches--proteomics, transcriptomics, metabolomics, and glycomics--have a strong record in several areas of biomedicine, including neurochemistry and neuro-oncology. Published applications of an integrated approach to projects of neurological, psychiatric, and pharmacological natures are still few but show promise to provide deep biological knowledge derived from cells, animal models, and clinical materials. Future studies that yield insights based on integrated analyses promise to deliver new therapeutic targets and biomarkers for personalized medicine.

Comprehensive analysis of genic male sterility-related genes in Brassica rapa using a newly developed Br300K oligomeric chip

To identify genes associated with genic male sterility (GMS) that could be useful for hybrid breeding in Chinese cabbage (Brassicarapa ssp. pekinensis), floral bud transcriptome analysis was carried out using a B. rapa microarray with 300,000 probes (Br300K). Among 47,548 clones deposited on a Br300K microarray with seven probes of 60 nt length within the 3' 150 bp region, a total of 10,622 genes were differentially expressed between fertile and sterile floral buds; 4,774 and 5,848 genes were up-regulated over 2-fold in fertile and sterile buds, respectively. However, the expression of 1,413 and 199 genes showed fertile and sterile bud-specific features, respectively. Genes expressed specifically in fertile buds, possibly GMS-related genes, included homologs of several Arabidopsis male sterility-related genes, genes associated with the cell wall and synthesis of its surface proteins, pollen wall and coat components, signaling components, and nutrient supplies. However, most early genes for pollen development, genes for primexine and callose formation, and genes for pollen maturation and anther dehiscence showed no difference in expression between fertile and sterile buds. Some of the known genes associated with Arabidopsis pollen development showed similar expression patterns to those seen in this study, while others did not. BrbHLH89 and BrMYP99 are putative GMS genes. Additionally, 17 novel genes identified only in B. rapa were specifically and highly expressed only in fertile buds, implying the possible involvement in male fertility. All data suggest that Chinese cabbage GMS might be controlled by genes acting in post-meiotic tapetal development that are different from those known to be associated with Arabidopsis male sterility.

Plants know where it hurts: root and shoot jasmonic acid induction elicit differential responses in Brassica oleracea

Plants respond to herbivore attack by rapidly inducing defenses that are mainly regulated by jasmonic acid (JA). Due to the systemic nature of induced defenses, attack by root herbivores can also result in a shoot response and vice versa, causing interactions between above- and belowground herbivores. However, little is known about the molecular mechanisms underlying these interactions. We investigated whether plants respond differently when roots or shoots are induced. We mimicked herbivore attack by applying JA to the roots or shoots of Brassica oleracea and analyzed molecular and chemical responses in both organs. In shoots, an immediate and massive change in primary and secondary metabolism was observed. In roots, the JA-induced response was less extensive and qualitatively different from that in the shoots. Strikingly, in both roots and shoots we also observed differential responses in primary metabolism, development as well as defense specific traits depending on whether the JA induction had been below- or aboveground. We conclude that the JA response is not only tissue-specific but also dependent on the organ that was induced. Already very early in the JA signaling pathway the differential response was observed. This indicates that both organs have a different JA signaling cascade, and that the signal eliciting systemic responses contains information about the site of induction, thus providing plants with a mechanism to tailor their responses specifically to the organ that is damaged.

Unstable transcripts in Arabidopsis allotetraploids are associated with nonadditive gene expression in response to abiotic and biotic stresses

Genome-wide analysis has documented differential gene expression between closely related species in plants and animals and nonadditive gene expression in hybrids and allopolyploids compared to the parents. In Arabidopsis, 15–43% of genes are expressed differently between the related species, Arabidopsis thaliana and Arabidopsis arenosa, the majority of which are nonadditively expressed (differently from mid-parent value) in allotetraploids. Nonadditive gene expression can be caused by transcriptional regulation through chromatin modifications, but the role of posttranscriptional regulation in nonadditive gene expression is largely unknown. Here we reported genome-wide analysis of mRNA decay in resynthesized Arabidopsis allotetraploids. Among ∼26,000 annotated genes, over 1% of gene transcripts showed rapid decay with an estimated half-life of less than 60 minutes, and they are called allotetraploid genes with unstable transcripts (AlloGUTs). Remarkably, 30% of alloGUTs matched the nonadditively expressed genes, and their expression levels were negatively correlated with the decay rate. Compared to all genes, these nonadditively expressed alloGUTs were overrepresented 2-6-fold in the Gene Ontology (GOSlim) classifications in response to abiotic and biotic stresses, signal transduction, and transcription. Interestingly, the AlloGUTs include transcription factor genes that are highly inducible under stress conditions and circadian clock regulators that regulate growth in A. thaliana. These data suggest a role of mRNA stability in homoeologous gene expression in Arabidopsis allopolyploids. The enrichment of nonadditively expressed genes in stress-related pathways were commonly observed in Arabidopsis and other allopolyploids such as wheat and cotton, which may suggest a role for stress-mediated growth vigor in hybrids and allopolyploids.

Towards unambiguous transcript mapping in the allotetraploid Brassica napusThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”

The architecture of the Brassica napus genome is marked by its evolutionary origins. The genome of B. napus was formed from the hybridization of two closely related diploid Brassica species, both of which evolved from an hexaploid ancestor. The extensive whole genome duplication events in its near and distant past result in the allotetraploid genome of B. napus maintaining multiple copies of most genes, which predicts a highly complex and redundant transcriptome that can confound any expression analyses. A stringent assembly of 142 399 B. napus expressed sequence tags allowed the development of a well-differentiated set of reference transcripts, which were used as a foundation to assess the efficacy of available tools for identifying and distinguishing transcripts in B. napus ; including microarray hybridization and 3′ anchored sequence tag capture. Microarray platforms cannot distinguish transcripts derived from the two progenitors or close homologues, although observed differential expression appeared to be biased towards unique transcripts. The use of 3′ capture enhanced the ability to unambiguously identify homologues within the B. napus transcriptome but was limited by tag length. The ability to comprehensively catalogue gene expression in polyploid species could be transformed by the application of cost-efficient next generation sequencing technologies that will capture millions of long sequence tags.

Intraspecific variation in herbivore community composition and transcriptional profiles in field-grown Brassica oleracea cultivars

Intraspecific differences in plant defence traits are often correlated with variation in transcriptional profiles and can affect the composition of herbivore communities on field-grown plants. However, most studies on transcriptional profiling of plant-herbivore interactions have been carried out under controlled conditions in the laboratory or greenhouse and only a few examine intraspecific transcriptional variation. Here, intraspecific variation in herbivore community composition and transcriptional profiles between two Brassica oleracea cultivars grown in the field is addressed. Early in the season, no differences in community composition were found for naturally occurring herbivores, whereas cultivars differed greatly in abundance, species richness, and herbivore community later in the season. Genome-wide transcriptomic analysis using an Arabidopsis thaliana oligonucleotide microarray showed clear differences for the expression levels of 26 genes between the two cultivars later in the season. Several defence-related genes showed higher levels of expression in the cultivar that harboured the lowest numbers of herbivores. Our study shows that herbivore community composition develops differentially throughout the season on the two B. oleracea cultivars grown in the field. The correlation between the differences in herbivore communities and differential expression of particular defence-related genes is discussed.

A separation defect of tapetum cells and microspore mother cells results in male sterility in Brassica napus: the role of abscisic acid in early anther development

Male sterility is an important contributor to heterosis in Brassica napus L. The B. napus line 7-7365ABC is a recessive epistatic genic male sterile (REGMS) three-line system. The 7-7365A line with the genotype Bnms3ms3ms4ms4RfRf is male-sterile, while the 7-7365B line with the genotype BnMs3ms3ms4ms4RfRf is male-fertile, and 7-7365C with homozygous recessive genotypes at the three loci shows male fertility because the loss function of Bnrf gene causes the inhibition of the genetic trait of the double mutant Bnms3 Bnms4. Histological studies addressing male sterility, transcriptional regulation pathways and the role of abscisic acid (ABA) in the anther development of REGMS plants are reported here. In the male-sterile line 7-7365A, tapetum cell and microspore mother cell separation were affected, and this led to failure of microspore release. The activity of polygalacturonase and the expression of the pectin methylesterase gene (AT3g06830) were significantly downregulated. Nine genes were downregulated in 7-7365A compared to 7-7365B and 7-7365C, including genes specifically expressed in tapetum (A3, A9, MS1) and the ABA-responsive gene KIN1. ABA concentration in 7-7365B was significantly higher than in 7-7365A and 7-7365C in young flower buds. Furthermore, temperature treatment made some sterile 7-7365A flowers become fertile. The stamens in these flowers produced viable pollen, and filament elongation was restored to its level in 7-7365C. We propose that ABA might control the expression of genes involved in cell separation during early anther development. The REGMS phenotype could be controlled by a primary pathway of male sterile metabolism positively regulated by the BnMs3 gene and a supplementary pathway negatively regulated by the BnRf gene.

Down-regulation of defense genes and resource allocation into infected roots as factors for compatibility between Fagus sylvatica and Phytophthora citricola

Phytophthora citricola is a wide spread and highly aggressive pathogen of Fagus sylvatica. The hemibiotrophic oomycete infects the roots and establishes a compatible interaction with F. sylvatica. To investigate the transcriptional changes associated with P. citricola infection, 68 custom oligo-microarray measurements were conducted. Hierarchical as well as non-hierarchical clustering was carried out to analyze the expression profiles. Experimental setup includes a time scale covering the biotrophic and necrotrophic stages of interaction as well as comparative analyses of the local and systemic responses. The local reaction of F. sylvatica is characterized by a striking lack of defense gene induction leading to the conclusion that P. citricola escapes the main recognition systems and/or suppresses the host's response. The analysis of the systemic reaction revealed a massive shift in gene expression patterns during the biotrophic phase that is interpreted as evidence of resource allocation into the roots to support the increased sink caused by pathogen growth. Defense genes known to be responsive to salicylic acid (effective against biotrophs), jasmonic acid, and ethylene (effective against necrotrophs and herbivores) are represented on the arrays. All significant changes in gene expression measured for salicylic acid responsive genes were down-regulations in roots and leaves while some jasmonic acid responsive genes showed a very late up-regulation only in leaves, probably caused by the desiccation shortly before plant death. Together, these expression changes could explain the success of the pathogen.

Patterns of differential gene expression in Brassica napus cultivars infected with Sclerotinia sclerotiorum

SUMMARY The fungal pathogen Sclerotinia sclerotiorum infects a broad range of dicotyledonous plant species and causes stem rot in Brassica napus. To elucidate the mechanisms underlying the defence response, the patterns of gene expression in the partially resistant B. napus cultivar ZhongYou 821 (ZY821) and the susceptible cultivar Westar were studied using a B. napus oligonucleotide microarray. Although maximum differential gene expression was observed at 48 h post-inoculation (hpi) in both cultivars, increased transcript levels were detected in cv. ZY821 at the earlier stages of infection (6-12 hpi) for many genes, including those encoding defence-associated proteins, such as chitinases, glucanases, osmotins and lectins, as well as genes encoding transcription factors belonging to the zinc finger, WRKY, APETALA2 (AP2) and MYB classes. In both cultivars, genes encoding enzymes involved in jasmonic acid, ethylene and auxin synthesis were induced, as were those for gibberellin degradation. In addition, changes in the expression of genes encoding enzymes involved in carbohydrate and energy metabolism appeared to be directed towards shuttling carbon reserves to the tricarboxylic acid cycle and generating reactive oxygen species. Transcripts from genes encoding enzymes involved in glucosinolate and phenylpropanoid biosynthesis were highly elevated in both cultivars, suggesting that secondary metabolites are also components of the response to S. sclerotiorum in B. napus.

Analysis of gene expression profile in pollen development of recessive genic male sterile Brassica napus L. line S45A

Male sterility in a near-isogenic line S45AB after 25 generations of subcrossing is controlled by two pairs of duplicate genes. The genotype of S45A is Bnms1Bnms1Bnms2Bnms2, and that of S45B is BnMs1Bnms1Bnms2Bnms2, respectively. Histological observations revealed that abnormal anther development appeared in the tapetum and pollen exine during the tetrad stage. This male sterility was characterized by hypertrophy of the tapetal cells at the tetrad stage and a complete lack of microspore exine after the release of microspores from the tetrads. To elucidate the mechanism of this recessive genic male sterility, the flower bud expression profiles of the S45A and S45B lines were analyzed using an Arabidopsis thaliana ATH1 oligonucleotide array. When compared with the S45B line, 69 genes were significantly downregulated, and 46 genes were significantly upregulated in the S45A line. Real-time polymerase chain reaction (PCR) was then used to verify the results of the microarray analysis, and the majority of the downregulated genes in the S45A line were abundantly and specifically expressed in the anther. The results of the real-time PCR suggest that Bnms1 might be involved in the metabolism of lipid/fatty acids, and the homologous mutation of Bnms1 may either block the biosynthesis of sporopollenin or block sporopollenin from being deposited on the microspore surface, thus, preventing pollen exine formation. The role of Bnms1 in the regulatory network of exine formation is also discussed as well.

Analysis of gene expression in resynthesized Brassica napus Allopolyploids using arabidopsis 70mer oligo microarrays

Background

Studies in resynthesized Brassica napus allopolyploids indicate that homoeologous chromosome exchanges in advanced generations (S_5∶6) alter gene expression through the loss and doubling of homoeologous genes within the rearrangements. Rearrangements may also indirectly affect global gene expression if homoeologous copies of gene regulators within rearrangements have differential affects on the transcription of genes in networks.

Methodology/Principal Findings

We utilized Arabidopsis 70mer oligonucleotide microarrays for exploring gene expression in three resynthesized B. napus lineages at the S_0∶1 and S_5∶6 generations as well as their diploid progenitors B. rapa and B. oleracea. Differential gene expression between the progenitors and additive (midparent) expression in the allopolyploids were tested. The S_5∶6 lines differed in the number of genetic rearrangements, allowing us to test if the number of genes displaying nonadditive expression was related to the number of rearrangements. Estimates using per-gene and common variance ANOVA models indicated that 6–15% of 26,107 genes were differentially expressed between the progenitors. Individual allopolyploids showed nonadditive expression for 1.6–32% of all genes. Less than 0.3% of genes displayed nonadditive expression in all S_0∶1 lines and 0.1–0.2% were nonadditive among all S_5∶6 lines. Differentially expressed genes in the polyploids were over-represented by genes differential between the progenitors. The total number of differentially expressed genes was correlated with the number of genetic changes in S_5∶6 lines under the common variance model; however, there was no relationship using a per-gene variance model, and many genes showed nonadditive expression in S_0∶1 lines.

Conclusions/Significance

Few genes reproducibly demonstrated nonadditive expression among lineages, suggesting few changes resulted from a general response to polyploidization. Furthermore, our microarray analysis did not provide strong evidence that homoeologous rearrangements were a determinant of genome-wide nonadditive gene expression. In light of the inherent limitations of the Arabidopsis microarray to measure gene expression in polyploid Brassicas, further studies are warranted.

Altered circadian rhythms regulate growth vigour in hybrids and allopolyploids

Segregating hybrids and stable allopolyploids display morphological vigor1,2,3, and Arabidopsis allotetraploids are larger than the parents Arabidopsis thaliana and A. arenosa1,4. The mechanisms are unknown. Circadian clocks mediate metabolic pathways and increase fitness in animals and plants5,6,7,8. Here we report that epigenetic modifications of the circadian clock genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY)9,10 and their reciprocal regulators TIMING OF CAB EXPRESSION 1 (TOC1) and GIGANTEA (GI)10,11,12 mediate expression changes in downstream genes and pathways. During the day, epigenetic repression of CCA1 and LHY induced expression of TOC1, GI and downstream genes that contain CCA1 binding site (CBS)13 in chlorophyll and starch metabolic pathways in allotetraploids and F₁ hybrids, which produced more chlorophyll and starch than the parents in the same environment. Mutations in cca1 and cca1 lhy and daily repression of cca1 in TOC1:cca1-RNAi transgenic plants increased expression of downstream genes and chlorophyll and starch content, whereas constitutively expressing CCA1 or ectopically expressing TOC1:CCA1 had the opposite effects. The causal effects of CCA1 on output traits suggest that hybrids and allopolyploids gain advantages from the control of circadian-mediated physiological and metabolic pathways, leading to growth vigor and increased biomass.

Responses of Brassica oleracea cultivars to infestation by the aphid Brevicoryne brassicae: an ecological and molecular approach

Intraspecific variation in resistance or susceptibility to herbivorous insects has been widely studied through bioassays. However, few studies have combined this with a full transcriptomic analysis. Here, we take such an approach to study the interaction between the aphid Brevicoryne brassicae and four white cabbage (Brassica oleracea var. capitata) cultivars. Both under glasshouse and field conditions, two of the cultivars clearly supported a faster aphid population development than the other two, indicating that aphid population development was largely independent of the environmental conditions. Genome-wide transcriptomic analysis using 70-mer oligonucleotide microarrays based on the Arabidopsis thaliana genome showed that only a small number of genes were differentially regulated, and that this regulation was highly cultivar specific. The temporal pattern in the expression behaviour of two B. brassicae-responsive genes in all four cultivars together with targeted studies employing A. thaliana knockout mutants revealed a possible role for a trypsin-and-protease inhibitor in defence against B. brassicae. Conversely, a xyloglucan endotransglucosylase seemed to have no effect on aphid performance. Overall, this study shows clear intraspecific variation in B. brassicae susceptibility among B. oleracea cultivars under glasshouse and field conditions that can be partly explained by certain differences in induced transcriptional changes.

Male-derived butterfly anti-aphrodisiac mediates induced indirect plant defense

Plants can recruit parasitic wasps in response to egg deposition by herbivorous insects-a sophisticated indirect plant defense mechanism. Oviposition by the Large Cabbage White butterfly Pieris brassicae on Brussels sprout plants induces phytochemical changes that arrest the egg parasitoid Trichogramma brassicae. Here, we report the identification of an elicitor of such an oviposition-induced plant response. Eliciting activity was present in accessory gland secretions released by mated female butterflies during egg deposition. In contrast, gland secretions from virgin female butterflies were inactive. In the male ejaculate, P. brassicae females receive the anti-aphrodisiac benzyl cyanide (BC) that reduces the females' attractiveness for subsequent mating. We detected this pheromone in the accessory gland secretion released by mated female butterflies. When applied onto leaves, BC alone induced phytochemical changes that arrested females of the egg parasitoid. Microarray analyses revealed a similarity in induced plant responses that may explain the arrest of T. brassicae to egg-laden and BC-treated plants. Thus, a male-derived compound endangers the offspring of the butterfly by inducing plant defense. Recently, BC was shown to play a role in foraging behavior of T. brassicae, by acting as a cue to facilitate phoretic transport by mated female butterflies to oviposition sites. Our results suggest that the anti-aphrodisiac pheromone incurs fitness costs for the butterfly by both mediating phoretic behavior and inducing plant defense.

RNAi of met1 reduces DNA methylation and induces genome-specific changes in gene expression and centromeric small RNA accumulation in Arabidopsis allopolyploids

Changes in genome structure and gene expression have been documented in both resynthesized and natural allopolyploids that contain two or more divergent genomes. The underlying mechanisms for rapid and stochastic changes in gene expression are unknown. Arabidopsis suecica is a natural allotetraploid derived from the extant A. thaliana and A. arenosa genomes that are homeologous in the allotetraploid. Here we report that RNAi of met1 reduced DNA methylation and altered the expression of approximately 200 genes, many of which encode transposons, predicted proteins, and centromeric and heterochromatic RNAs. Reduced DNA methylation occurred frequently in promoter regions of the upregulated genes, and an En/Spm-like transposon was reactivated in met1-RNAi A. suecica lines. Derepression of transposons, heterochromatic repeats, and centromeric small RNAs was primarily derived from the A. thaliana genome, and A. arenosa homeologous loci were less affected by methylation defects. A high level of A. thaliana centromeric small RNA accumulation was correlated with hypermethylation of A. thaliana centromeres. The greater effects of reduced DNA methylation on transposons and centromeric repeats in A. thaliana than in A. arenosa are consistent with the repression of many genes that are expressed at higher levels in A. thaliana than in A. arenosa in the resynthesized allotetraploids. Moreover, non-CG (CC) methylation in the promoter region of A. thaliana At2g23810 remained in the resynthesized allotetraploids, and the methylation spread within the promoter region in natural A. suecica, leading to silencing of At2g23810. At2g23810 was demethylated and reactivated in met1-RNAi A. suecica lines. We suggest that many A. thaliana genes are transcriptionally repressed in resynthesized allotetraploids, and a subset of A. thaliana loci including transposons and centromeric repeats are heavily methylated and subjected to homeologous genome-specific RNA-mediated DNA methylation in natural allopolyploids.

Global analysis of gene expression in flower buds of Ms-cd1 Brassica oleracea conferring male sterility by using an Arabidopsis microarray

The dominant male sterility gene Ms-cd1 is identified in Brassica oleracea. Electron microscopical observations revealed that abortion of pollen development starts after tetrad formation. This important male sterility phenotype is characterized by lack of degradation of the primary pollen mother cell (PMC) wall and delayed degradation of callose surrounding the tetrads and thus arrest of microspore release. Gene expression of the male sterile and fertile buds was analyzed by heterologous hybridization of Brassica oleracea cRNA onto an Arabidopsis whole genome oligonucleotide microarray. A total of 277 suppressed genes including 40 kinase-, 32 cell wall modification and 29 transport related genes were found to be significantly down regulated >3-fold in the male sterile mutant. The vast majority of the differentially expressed transcripts are found to present late pollen stage specific genes. Kinase genes, cell wall modification genes and ion transport genes were greatly over-represented when compared to their percentage of all flower bud expressed genes and represent 36.5% of the genes suppressed by Ms-cd1. Our results also suggest that Ms-cd1 may blocks an anther developmental pathway with a small number of genes suppressed in tapetum cells which prevent the degradation of callose and PMC wall, which further leads to the suppression of a large number of genes involved in signaling pathways, cell wall modification and ion transport in pollen grains.

Analysis of gene expression profiles in response to Sclerotinia sclerotiorum in Brassica napus

Sclerotinia sclerotiorum is a necrotrophic plant pathogen which causes serious disease in agronomically important crop species. The molecular basis of plant defense to this pathogen is poorly understood. We investigated gene expression changes associated with S. sclerotiorum infection in a partially resistant and a susceptible genotype of oilseed Brassica napus using a whole genome microarray from Arabidopsis. A total of 686 and 1,547 genes were found to be differentially expressed after infection in the resistant and susceptible genotypes, respectively. The number of differentially expressed genes increased over infection time with the majority being up-regulated in both genotypes. The putative functions of the differentially expressed genes included pathogenesis-related (PR) proteins, proteins involved in the oxidative burst, protein kinase, molecule transporters, cell maintenance and development, abiotic stress, as well as proteins with unknown functions. The gene regulation patterns indicated that a large part of the defense response exhibited as a temporal and quantitative difference between the two genotypes. Genes associated with jasmonic acid (JA) and ethylene signal transduction pathways were induced, but no salicylic acid (SA) responsive genes were identified. Candidate defense genes were identified by integration of the early response genes in the partially resistant line with previously mapped quantitative trait loci (QTL). Expression levels of these genes were verified by Northern blot analyses. These results indicate that genes encoding various proteins involved in diverse roles, particularly WRKY transcription factors and plant cell wall related proteins may play an important role in the defense response to S. sclerotiorum disease.

Plant interactions with microbes and insects: from molecular mechanisms to ecology

Plants are members of complex communities and interact both with antagonists and beneficial organisms. An important question in plant defense-signaling research is how plants integrate signals induced by pathogens, beneficial microbes and insects into the most appropriate adaptive response. Molecular and genomic tools are now being used to uncover the complexity of the induced defense signaling networks that have evolved during the arms races between plants and their attackers. Molecular biologists and ecologists are joining forces to place molecular mechanisms of plant defense into an ecological perspective. Here, we review our current understanding of the molecular mechanisms of induced plant defense and their potential ecological relevance in nature.

Sensitivity and speed of induced defense of cabbage (Brassica oleracea L.): dynamics of BoLOX expression patterns during insect and pathogen attack

The lipoxygenase pathway is involved in the early steps of plant responses to herbivorous insects and phytopathogens. Induced defenses in the crucifer Brassica oleracea have been well documented. Here, we have cloned a LIPOXYGENASE (LOX) from B. oleracea (BoLOX). The sequence reveals that the BoLOX protein has a transit peptide for chloroplast targeting, which is characteristic for class 2 LOXs involved in jasmonic acid (JA) biosynthesis which takes place in the chloroplast. Phylogenetic analysis shows that BoLOX is closely related to B. napus BnLOX2fl and Arabidopsis thaliana AtLOX2, which mediates JA biosynthesis. BoLOX also shares functional characteristics with AtLOX2; BoLOX is inducible by wounding, JA treatment, and herbivores such as caterpillars (Pieris rapae, P. brassicae, and Mamestra brassicae), spider mites (Tetranychus urticae), locusts (Schistocerca gregaria), and a bacterial pathogen (Pseudomonas syringae pv. tomato). Of these, Pieris spp. caterpillars also induce AtLOX2 and JA biosynthesis in Arabidopsis. However, the aphid Myzus persicae did not induce BoLOX, which agrees with previous reports that this aphid induces neither AtLOX2 nor JA biosynthesis in Arabidopsis. Quantitative expression analysis of temporal, spatial, and density-dependent BoLOX transcript levels through real-time quantitative polymerase chain reaction demonstrated that BoLOX is maximally expressed after feeding by only two first-instar caterpillars for 24 h. Systemic expression was approximately 10-fold lower than local expression for herbivore-induced responses. The good correlation of BoLOX transcript levels with reports in the literature on induced defenses of B. oleracea is discussed.

Genotypic variation in genome-wide transcription profiles induced by insect feeding: Brassica oleracea--Pieris rapae interactions

Background

Transcriptional profiling after herbivore attack reveals, at the molecular level, how plants respond to this type of biotic stress. Comparing herbivore-induced transcriptional responses of plants with different phenotypes provides insight into plant defense mechanisms. Here, we compare the global gene expression patterns induced by Pieris rapae caterpillar attack in two white cabbage (Brassica oleracea var. capitata) cultivars. The two cultivars are shown to differ in their level of direct defense against caterpillar feeding. Because Brassica full genome microarrays are not yet available, 70-mer oligonucleotide microarrays based on the Arabidopsis thaliana genome were used for this non-model plant.

Results

The transcriptional responses of the two cultivars differed in timing as characterized by changes in their expression pattern after 24, 48 and 72 hours of caterpillar feeding. In addition, they also differed qualitatively. Surprisingly, of all genes induced at any time point, only one third was induced in both cultivars. Analyses of transcriptional responses after jasmonate treatment revealed that the difference in timing did not hold for the response to this phytohormone. Additionally, comparisons between Pieris rapae- and jasmonate-induced transcriptional responses showed that Pieris rapae induced more jasmonate-independent than jasmonate-dependent genes.

Conclusion

The present study clearly shows that global transcriptional responses in two cultivars of the same plant species in response to insect feeding can differ dramatically. Several of these differences involve genes that are known to have an impact on Pieris rapae performance and probably underlie different mechanisms of direct defense, present in the cultivars.

Microarray analysis reveals altered expression of a large number of nuclear genes in developing cytoplasmic male sterile Brassica napus flowers

To gain new insights into the mechanism underlying cytoplasmic male sterility (CMS), we compared the nuclear gene expression profiles of flowers of a Brassica napus CMS line with that of the fertile B. napus maintainer line using Arabidopsis thaliana flower-specific cDNA microarrays. The CMS line used has a B. napus nuclear genome, but has a rearranged mitochondrial (mt) genome consisting of both B. napus and A. thaliana DNA. Gene expression profiling revealed that a large number of genes differed in expression between the two lines. For example, nuclear genes coding for proteins that are involved in protein import into organelles, genes expressed in stamens and pollen, as well as genes implicated in either cell-wall remodeling or architecture, were repressed in the CMS line compared with B. napus. These results show that the mt genome of the CMS line strongly influences nuclear gene expression, and thus reveal the importance of retrograde signalling between the mitochondria and the nucleus. Furthermore, flowers of the CMS line are characterized by a replacement of stamens with carpelloid organs, and thus partially resemble the APETALA3 (AP3) and PISTILLATA (PI) mutants. In accordance with this phenotype, AP3 expression was downregulated in the stamens, shortly before these organs developed carpelloid characteristics, even though it was initiated correctly. Repression of PI succeeded that of AP3 and might be a consequence of a loss of AP3 activity. These results suggest that AP3 expression in stamens depends on proper mt function and a correct nuclear-mt interaction, and that mt alterations cause the male sterility phenotype of the CMS line.

Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids

Polyploidy, or whole-genome duplication (WGD), is an important genomic feature for all eukaryotes, especially many plants and some animals. The common occurrence of polyploidy suggests an evolutionary advantage of having multiple sets of genetic material for adaptive evolution. However, increased gene and genome dosages in autopolyploids (duplications of a single genome) and allopolyploids (combinations of two or more divergent genomes) often cause genome instabilities, chromosome imbalances, regulatory incompatibilities, and reproductive failures. Therefore, new allopolyploids must establish a compatible relationship between alien cytoplasm and nuclei and between two divergent genomes, leading to rapid changes in genome structure, gene expression, and developmental traits such as fertility, inbreeding, apomixis, flowering time, and hybrid vigor. Although the underlying mechanisms for these changes are poorly understood, some themes are emerging. There is compelling evidence that changes in DNA sequence, cis- and trans-acting effects, chromatin modifications, RNA-mediated pathways, and regulatory networks modulate differential expression of homoeologous genes and phenotypic variation that may facilitate adaptive evolution in polyploid plants and domestication in crops.

Exploiting the wild crucifer Thlaspi arvense to identify conserved and novel genes expressed during a plant's response to cold stress

Thlaspi arvense, a wild species from the Brassicaceae family, was shown to have a higher level of freezing tolerance than either of its close relatives, the model plant Arabidopsis thaliana or the crop Brassica napus (canola). Over 600 clones were sequenced from a subtractive cDNA library generated from cold treated T. arvense tissue, establishing that T. arvense shared significant sequence identity with both A. thaliana and B. napus (90-92%). In light of the strong sequence similarity between T. arvense and A. thaliana and to exploit the available genomics resources for Arabidopsis, the efficacy of using long 70 mer oligonucleotide whole genome Arabidopsis microarrays was tested for T. arvense. Gene expression in T. arvense leaf tissue during the very early stages of cold acclimation (or cold stress) was assayed at three time points and compared to an untreated control. This analysis highlights some of the difficulties and benefits of using cross-species microarray analysis. The data suggested that T. arvense responds in a similar fashion to cold stress as the model plant A. thaliana. However, for a number of genes quantitative differences in the level and timing of expression were identified. One of the most notable differences suggested that sulphur assimilation leading to the increased production of the methyl donor S-adenosyl-methionine was playing a role in the response of T. arvense to cold stress.

Natural genetic variation in whole-genome expression in Arabidopsis thaliana: the impact of physiological QTL introgression

A long-standing and fundamental question in biology is how genes influence complex phenotypes. Combining near-isogenic line mapping with genome expression profiling offers a unique opportunity for exploring the functional relationship between genotype and phenotype and for generating candidate genes for future study. We used a whole-genome microarray produced with ink-jet technology to measure the relative expression level of over 21,500 genes from an Arabidopsis thaliana near-isogenic line (NIL) and its recurrent parent. The NIL material contained two introgressions (bottom of chromosome II and top of chromosome III) of the Cvi-1 ecotype in a Ler-2 ecotype genome background. Each introgression 'captures' a Cvi allele of a physiological quantitative trait loci (QTL) that our previous studies have shown increases transpiration and reduces water-use efficiency at the whole-plant level. We used a mixed model anova framework for assessing sources of expression variability and for evaluating statistical significance in our array experiment. We discovered 25 differentially expressed genes in the introgression at a false-discovery rate (FDR) cut-off of 0.20 and identified new candidate genes for both QTL regions. Several differentially expressed genes were confirmed with QRT-PCR (quantitative reverse transcription-polymerase chain reaction) assays. In contrast, we found no statistically significant differentially expressed genes outside of the QTL introgressions after controlling for multiple tests. We discuss these results in the context of candidate genes, cloning QTL, and phenotypic evolution.

Microarray analysis reveals differential gene expression in hybrid sunflower species

This paper describes the creation of a cDNA microarray for annual sunflowers and its use to elucidate patterns of gene expression in Helianthus annuus, Helianthus petiolaris, and the homoploid hybrid species Helianthus deserticola. The array comprises 3743 ESTs (expressed sequence tags) representing approximately 2897 unique genes. It has an average clone/EST identity rate of 91%, is applicable across species boundaries within the annual sunflowers, and shows patterns of gene expression that are highly reproducible according to real-time RT-PCR (reverse transcription-polymerase chain reaction) results. Overall, 12.8% of genes on the array showed statistically significant differential expression across the three species. Helianthus deserticola displayed transgressive, or extreme, expression for 58 genes, with roughly equal numbers exhibiting up- or down-regulation relative to both parental species. Transport-related proteins were strongly over-represented among the transgressively expressed genes, which makes functional sense given the extreme desert floor habitat of H. deserticola. The potential adaptive value of differential gene expression was evaluated for five genes in two populations of early generation (BC2) hybrids between the parental species grown in the H. deserticola habitat. One gene (a G protein-coupled receptor) had a significant association with fitness and maps close to a QTL controlling traits that may be adaptive in the desert habitat.

Numerous and rapid nonstochastic modifications of gene products in newly synthesized Brassica napus allotetraploids

Polyploidization is a widespread process that results in the merger of two or more genomes in a common nucleus. To investigate modifications of gene expression occurring during allopolyploid formation, the Brassica napus allotetraploid model was chosen. Large-scale analyses of the proteome were conducted on two organs, the stem and root, so that >1600 polypeptides were screened. Comparative proteomics of synthetic B. napus and its homozygous diploid progenitors B. rapa and B. oleracea showed that very few proteins disappeared or appeared in the amphiploids (<1%), but a strikingly high number (25-38%) of polypeptides displayed quantitative nonadditive pattern. Nonstochastic gene expression repatterning was found since 99% of the detected variations were reproducible in four independently created amphiploids. More than 60% of proteins displayed a nonadditive pattern closer to the paternal parent B. rapa. Interspecific hybridization triggered the majority of the deviations (89%), whereas very few variations (approximately 3%) were associated with genome doubling and more significant alterations arose from selfing (approximately 9%). Some nonadditive proteins behaved similarly in both organs, while others exhibited contrasted behavior, showing rapid organ-specific regulation. B. napus formation was therefore correlated with immediate and directed nonadditive changes in gene expression, suggesting that the early steps of allopolyploidization repatterning are controlled by nonstochastic mechanisms.

Developmental and gene expression analyses of a cotton naked seed mutant

Cotton fiber development is a fundamental biological phenomenon, yet the molecular basis of fiber cell initiation is poorly understood. We examined molecular and cellular events of fiber cell development in the naked seed mutant (N1N1) and its isogenic line of cotton (Gossypium hirsutum L. cv. Texas Marker-1, TM-1). The dominant mutation not only delayed the process of fiber cell formation and elongation but also reduced the total number of fiber cells, resulting in sparsely distributed short fibers. Gene expression changes in TM-1 and N1N1 mutant lines among four tissues were analyzed using spotted cotton oligo-gene microarrays. Using the Arabidopsis genes, we selected and designed approximately 1,334 70-mer oligos from a subset of cotton fiber ESTs. Statistical analysis of the microarray data indicates that the number of significantly differentially expressed genes was 856 in the leaves compared to the ovules (3 days post-anthesis, DPA), 632 in the petals relative to the ovules (3 DPA), and 91 in the ovules at 0 DPA compared to 3 DPA, all in TM-1. Moreover, 117 and 30 genes were expressed significantly different in the ovules at three and 0 DPA, respectively, between TM-1 and N1N1. Quantitative RT-PCR analysis of 23 fiber-associated genes in seven tissues including ovules, fiber-bearing ovules, fibers, and non-fiber tissues in TM-1 and N1N1 indicates a mode of temporal regulation of the genes involved in transcriptional and translational regulation, signal transduction, and cell differentiation during early stages of fiber development. Suppression of the fiber-associated genes in the mutant may suggest that the N1N1 mutation disrupts temporal regulation of gene expression, leading to a defective process of fiber cell elongation and development.

Genomewide nonadditive gene regulation in Arabidopsis allotetraploids

Polyploidy has occurred throughout the evolutionary history of all eukaryotes and is extremely common in plants. Reunification of the evolutionarily divergent genomes in allopolyploids creates regulatory incompatibilities that must be reconciled. Here we report genomewide gene expression analysis of Arabidopsis synthetic allotetraploids, using spotted 70-mer oligo-gene microarrays. We detected >15% transcriptome divergence between the progenitors, and 2105 and 1818 genes were highly expressed in Arabidopsis thaliana and A. arenosa, respectively. Approximately 5.2% (1362) and 5.6% (1469) genes displayed expression divergence from the midparent value (MPV) in two independently derived synthetic allotetraploids, suggesting nonadditive gene regulation following interspecific hybridization. Remarkably, the majority of nonadditively expressed genes in the allotetraploids also display expression changes between the parents, indicating that transcriptome divergence is reconciled during allopolyploid formation. Moreover, >65% of the nonadditively expressed genes in the allotetraploids are repressed, and >94% of the repressed genes in the allotetraploids match the genes that are expressed at higher levels in A. thaliana than in A. arenosa, consistent with the silencing of A. thaliana rRNA genes subjected to nucleolar dominance and with overall suppression of the A. thaliana phenotype in the synthetic allotetraploids and natural A. suecica. The nonadditive gene regulation is involved in various biological pathways, and the changes in gene expression are developmentally regulated. In contrast to the small effects of genome doubling on gene regulation in autotetraploids, the combination of two divergent genomes in allotetraploids by interspecific hybridization induces genomewide nonadditive gene regulation, providing a molecular basis for de novo variation and allopolyploid evolution.

Conservation and divergence of light-regulated genome expression patterns during seedling development in rice and Arabidopsis

Genome-wide 70-mer oligonucleotide microarrays of rice (Oryza sativa) and Arabidopsis thaliana were used to profile genome expression changes during light-regulated seedling development. We estimate that the expression of approximately 20% of the genome in both rice and Arabidopsis seedlings is regulated by white light. Qualitatively similar expression profiles from seedlings grown under different light qualities were observed in both species; however, a quantitatively weaker effect on genome expression was observed in rice. Most metabolic pathways exhibited qualitatively similar light regulation in both species with a few species-specific differences. Global comparison of expression profiles between rice and Arabidopsis reciprocal best-matched gene pairs revealed a higher correlation of genome expression patterns in constant light than in darkness, suggesting that the genome expression profile of photomorphogenesis is more conserved. Transcription factor gene expression under constant light exposure was poorly conserved between the two species, implying a faster-evolving rate of transcription factor gene expression in light-grown plants. Organ-specific expression profiles during seedling photomorphogenesis provide genome-level evidence for divergent light effects in different higher plant organs. Finally, overrepresentation of specific promoter motifs in root- and leaf-specific light-regulated genes in both species suggests that these cis-elements are important for gene expression responses to light.

Development of anonymous cDNA microarrays to study changes to the Senecio floral transcriptome during hybrid speciation

Interspecific hybridization is an important process through which abrupt speciation can occur. In recent years, genetic changes associated with hybrid speciation have been identified through a variety of techniques, including AFLP/SSR mapping, GISH/FISH and cDNA-AFLP differential display. However, progress in using microarray technology to analyse whole genome/transcriptome changes associated with hybrid speciation has been limited due to the lack of extensive sequence data for many hybrid species and the difficulties in extrapolating results from commercially available microarrays for model species onto nonmodel hybrid taxa. Increasingly therefore researchers studying nonmodel systems are turning to the development of 'anonymous' cDNA microarrays, where the time and cost of producing microarrays is reduced by printing unsequenced cDNA clones, and sequencing only those clones that display interesting expression patterns. Here we describe the creation, testing and preliminary use of anonymous cDNA microarrays to study changes in floral transcriptome associated with allopolyploid speciation in the genus Senecio. We report a comparison of gene expression between the allohexaploid hybrid, Senecio cambrensis, its parental taxa Senecio squalidus (diploid) and Senecio vulgaris (tetraploid), and the intermediate triploid (sterile) hybrid Senecioxbaxteri. Anonymous microarray analysis revealed dramatic differences in floral gene expression between these four taxa and demonstrates the power of this technique for studies of the genetic impact of hybridization in nonmodel flowering plants.

Designing experiments using spotted microarrays to detect gene regulation differences within and among species

Comparative studies of genome-wide gene expression must account for variation not only among species, but also within species. Such studies are necessarily large in scale, because they incorporate experiments on multiple individuals of multiple species in multiple developmental stages in multiple environmental conditions. If the experiments are carefully designed and performed, the data they provide are worth the effort. We describe the utility of spotted microarrays for these studies and highlight experimental design criteria that will maximize inferential and statistical power. We conclude with a discussion of experimental protocols that are designed for investigations of differential gene expression and their pitfalls.

OsPTF1, a novel transcription factor involved in tolerance to phosphate starvation in rice

We report here on a novel transcription factor with a basic helix-loop-helix domain for tolerance to inorganic phosphate (Pi) starvation in rice (Oryza sativa). The gene is designated OsPTF1. The expression of OsPTF1 is Pi starvation induced in roots while constitutively expressed in shoots, as shown by northern-blot analysis. Overexpression of OsPTF1 enhanced tolerance to Pi starvation in transgenic rice. Tillering ability, root and shoot biomass, and phosphorus content of transgenic rice plants were about 30% higher than those of the wild-type plants in Pi-deficient conditions in hydroponic experiments. In soil pot and field experiments, more than 20% increase in tiller number, panicle weight, and phosphorus content was observed in transgenic plants compared to wild-type plants at low-Pi levels. In Pi-deficient conditions, transgenic rice plants showed significantly higher total root length and root surface area, which results in a higher instantaneous Pi uptake rate over their wild-type counterparts. Microarray analysis for transgenic plants overexpressing OsPTF1 has been performed to investigate the downstream regulation of OsPTF1.

A tiling microarray expression analysis of rice chromosome 4 suggests a chromosome-level regulation of transcription

The complete genome sequence of cultivated rice (Oryza sativa) provides an unprecedented opportunity to understand the biology of this model cereal. An essential and necessary step in this effort is the determination of the coding information and expression patterns of each sequenced chromosome. Here, we report an analysis of the transcriptional activity of rice chromosome 4 using a tiling path microarray based on PCR-generated genomic DNA fragments. Six representative rice organ types were examined using this microarray to catalog the transcribed regions of rice chromosome 4 and to reveal organ- and developmental stage-specific transcription patterns. This analysis provided expression support for 82% of the gene models in the chromosome. Transcriptional activities in 1643 nonannotated regions were also detected. Comparison with cytologically defined chromatin features indicated that in juvenile-stage rice the euchromatic region is more actively transcribed than is the transposon-rich heterochromatic portion of the chromosome. Interestingly, increased transcription of transposon-related gene models in certain heterochromatic regions was observed in mature-stage rice organs and in suspension-cultured cells. These results suggest a close correlation between transcriptional activity and chromosome organization and the developmental regulation of transcription activity at the chromosome level.

Three microarray platforms: an analysis of their concordance in profiling gene expression

Background

Microarrays for the analysis of gene expression are of three different types: short oligonucleotide (25–30 base), long oligonucleotide (50–80 base), and cDNA (highly variable in length). The short oligonucleotide and cDNA arrays have been the mainstay of expression analysis to date, but long oligonucleotide platforms are gaining in popularity and will probably replace cDNA arrays. As part of a validation study for the long oligonucleotide arrays, we compared and contrasted expression profiles from the three formats, testing RNA from six different cell lines against a universal reference standard.

Results

The three platforms had 6430 genes in common. In general, correlation of gene expression levels across the platforms was good when defined by concordance in the direction of expression difference (upregulation or downregulation), scatter plot analysis, principal component analysis, cell line correlation or quantitative RT-PCR. The overall correlations (r values) between platforms were in the range 0.7 to 0.8, as determined by analysis of scatter plots. When concordance was measured for expression ratios significant at p-values of <0.05 and at expression threshold levels of 1.5 and 2-fold, the agreement among the platforms was very high, ranging from 93% to 100%.

Conclusion

Our results indicate that the long oligonucleotide platform is highly suitable for expression analysis and compares favorably with the cDNA and short oligonucleotide varieties. All three platforms can give similar and reproducible results if the criterion is the direction of change in gene expression and minimal emphasis is placed on the magnitude of change.

Methods for genome-wide analysis of gene expression changes in polyploids

Polyploidy is an evolutionary innovation, providing extra sets of genetic material for phenotypic variation and adaptation. It is predicted that changes of gene expression by genetic and epigenetic mechanisms are responsible for novel variation in nascent and established polyploids (Liu and Wendel, 2002; Osborn et al., 2003; Pikaard, 2001). Studying gene expression changes in allopolyploids is more complicated than in autopolyploids, because allopolyploids contain more than two sets of genomes originating from divergent, but related, species. Here we describe two methods that are applicable to the genome-wide analysis of gene expression differences resulting from genome duplication in autopolyploids or interactions between homoeologous genomes in allopolyploids. First, we describe an amplified fragment length polymorphism (AFLP)--complementary DNA (cDNA) display method that allows the discrimination of homoeologous loci based on restriction polymorphisms between the progenitors. Second, we describe microarray analyses that can be used to compare gene expression differences between the allopolyploids and respective progenitors using appropriate experimental design and statistical analysis. We demonstrate the utility of these two complementary methods and discuss the pros and cons of using the methods to analyze gene expression changes in autopolyploids and allopolyploids. Furthermore, we describe these methods in general terms to be of wider applicability for comparative gene expression in a variety of evolutionary, genetic, biological, and physiological contexts.

Reversible histone acetylation and deacetylation mediate genome-wide, promoter-dependent and locus-specific changes in gene expression during plant development

Histone acetylation and deacetylation activate or repress transcription, yet the physiological relevance of reversible changes in chromatin structure and gene expression is poorly understood. We have shown that disrupting the expression of AtHD1 that encodes a putative Arabidopsis thaliana histone deacetylase induces a variety of developmental abnormalities. However, causal effects of the AtHD1 disruption on chromatin structure and gene expression are unknown. Using Arabidopsis spotted oligo-gene microarray analysis, here we report that >7% of the transcriptome was up- or downregulated in A. thaliana plants containing a T-DNA insertion in AtHD1 (athd1-t1), indicating that AtHD1 provides positive and negative control of transcriptional regulation. Remarkably, genes involved in ionic homeostasis and protein synthesis were ectopically expressed, whereas genes in ionic homeostasis, protein transport, and plant hormonal regulation were repressed in athd1-t1 leaves or flowers, suggesting a role of AtHD1 in developmental and environmental regulation of gene expression. Moreover, defective AtHD1 induced site-specific and reversible acetylation changes in H3-Lys9, H4-Lys12, and H4 tetra-lysines (residues 5, 8, 12, and 16) in homozygous recessive and heterozygous plants. Transcriptional activation was locus specific and often associated with specific acetylation sites in the vicinity of promoters, whereas gene repression did not correlate with changes in histone acetylation or correlated directly with H3-Lys9 methylation but not with DNA methylation. The data suggest that histone acetylation and deacetylation are promoter dependent, locus specific, and genetically reversible, which provides a general mechanism for reversible gene regulation responsive to developmental and environmental changes.

A rapid method for computationally inferring transcriptome coverage and microarray sensitivity

MOTIVATION

There are many different gene expression technologies, including cDNA and oligo-based microarrays, SAGE and MPSS. For each organism of interest, coverage of the transcriptome and the genome will be different. We address the question of what level of coverage is required to exploit the sensitivity of the different technologies, and what is the sensitivity of the different approaches in the experimental study.

RESULTS

We estimate the transcriptome coverage by randomly sampling transcripts from a pre-defined tag-to-gene mapping function. For a given microarray experiment, we locate the thresholds in intensities that define the distribution of transcript abundance. These values are compared against the distribution obtained by applying the same thresholds to the intensities from differentially expressed genes. The ratio of these two distributions meets at the equilibrium defining sensitivity. We conclude that a collection of approximately 340,000 sequences is adequate for microarrays, but not large enough for maximum utilization of tag-based technologies. In the absence of large-scale sequencing, the majority of the tags detected by the latter approaches will remain unidentified until the genome sequence is available.

The development of an Arabidopsis model system for genome-wide analysis of polyploidy effects

Arabidopsis is a model system not only for studying numerous aspects of plant biology, but also for understanding mechanisms of the rapid evolutionary process associated with genome duplication and polyploidization. Although in animals interspecific hybrids are often sterile and aneuploids are related to disease syndromes, both Arabidopsis autopolyploids and allopolyploids occur in nature and can be readily formed in the laboratory, providing an attractive system for comparing changes in gene expression and genome structure among relatively 'young' and 'established' or 'ancient' polyploids. Powerful reverse and forward genetics in Arabidopsis offer an exceptional means by which regulatory mechanisms of gene and genome duplication may be revealed. Moreover, the Arabidopsis genome is completely sequenced; both coding and non-coding sequences are available. We have developed spotted oligo-gene and chromosome microarrays using the complete Arabidopsis genome sequence. The oligo-gene microarray consists of ~26 000 70-mer oligonucleotides that are designed from all annotated genes in Arabidopsis, and the chromosome microarray contains 1 kb genomic tiling fragments amplified from a chromosomal region or the complete sequence of chromosome 4. We have demonstrated the utility of microarrays for genome-wide analysis of changes in gene expression, genome organization and chromatin structure in Arabidopsis polyploids and related species.