Genome-wide analysis of gene expression during early Arabidopsis flower development

Detailed information about stage-specific changes in gene expression is crucial for the understanding of the gene regulatory networks underlying development. Here, we describe the global gene expression dynamics during early flower development, a key process in the life cycle of a plant, during which floral patterning and the specification of floral organs is established. We used a novel floral induction system in Arabidopsis, which allows the isolation of a large number of synchronized floral buds, in conjunction with whole-genome microarray analysis to identify genes with differential expression at distinct stages of flower development. We found that the onset of flower formation is characterized by a massive downregulation of genes in incipient floral primordia, which is followed by a predominance of gene activation during the differentiation of floral organs. Among the genes we identified as differentially expressed in the experiment, we detected a significant enrichment of closely related members of gene families. The expression profiles of these related genes were often highly correlated, indicating similar temporal expression patterns. Moreover, we found that the majority of these genes is specifically up-regulated during certain developmental stages. Because co-expressed members of gene families in Arabidopsis frequently act in a redundant manner, these results suggest a high degree of functional redundancy during early flower development, but also that its extent may vary in a stage-specific manner.     

Multiple GUS expression patterns of a single Arabidopsis gene

Ten independent transposant lines with gene or enhancer traps (ET) inserted into the same gene (At2g01170) were identified in Arabidopsis thaliana . Transposon insertions were confirmed for each line. Only three of five ET lines and only one of the five gene trap (GT) lines displayed uidA (GUS) staining. The GUS (β-glucuronidase) expression patterns of the ET lines were different in all three lines. In the GT line, the GUS expression was restricted to the vascular tissue under all conditions examined. The variation in ET GUS expression suggests that each ET was controlled by different enhancer elements or the different elements of the trapped locus may give rise to different GUS expression patterns. Of five GT lines, three have the GUS gene in the same orientation as the At2g01170 open reading frame, yet only one yielded GUS staining. Regardless of the insertion construct, only those transposants with an insertion at the 3' end of the gene yielded GUS staining. Some transposants displayed a longer root phenotype in the presence of kanamycin that was also observed in 3' insertion sites in At2g01170. Taken together, these data show that insertions in the 5' end of the gene disrupted expression and emphasise the complexity encountered with ET and GT constructs to characterise the expression patterns of genes of interest based solely on GUS expression patterns.     

Genome-wide quantification of homeolog expression ratio revealed nonstochastic gene regulation in synthetic allopolyploid Arabidopsis

Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea (Arabidopsis petraea subsp. umbrosa), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica. The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated (r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.     

Genome-wide epigenetic perturbation jump-starts patterns of heritable variation found in nature

We extensively phenotyped 6000 Arabidopsis plants with experimentally perturbed DNA methylomes as well as a diverse panel of natural accessions in a common garden. We found that alterations in DNA methylation not only caused heritable phenotypic diversity but also produced heritability patterns closely resembling those of the natural accessions. Our findings indicate that epigenetically induced and naturally occurring variation in complex traits share part of their polygenic architecture and may offer complementary adaptation routes in ecological settings.     

Genome-wide analysis of gene expression profiles associated with cell cycle transitions in growing organs of Arabidopsis

Organ growth results from the progression of component cells through subsequent phases of proliferation and expansion before reaching maturity. We combined kinematic analysis, flowcytometry, and microarray analysis to characterize cell cycle regulation during the growth process of leaves 1 and 2 of Arabidopsis (Arabidopsis thaliana). Kinematic analysis showed that the epidermis proliferates until day 12; thereafter, cells expand until day 19 when leaves reach maturity. Flowcytometry revealed that endoreduplication occurs from the time cell division rates decline until the end of cell expansion. Analysis of 10 time points with a 6k-cDNA microarray showed that transitions between the growth stages were closely reflected in the mRNA expression data. Subsequent genome-wide microarray analysis on the three main stages allowed us to categorize known cell cycle genes into three major classes: constitutively expressed, proliferative, and inhibitory. Comparison with published expression data obtained from root zones corresponding to similar developmental stages and from synchronized cell cultures supported this categorization and enabled us to identify a high confidence set of 131 proliferation genes. Most of those had an M phase-dependent expression pattern and, in addition to many known cell cycle-related genes, there were at least 90 that were unknown or previously not associated with proliferation.     

Genome-wide analysis of lipoxygenase gene family in Arabidopsis and rice

The enzymes called lipoxygenases (LOXs) can dioxygenate unsaturated fatty acids, which leads to lipoperoxidation of biological membranes. This process causes synthesis of signaling molecules and also leads to changes in cellular metabolism. LOXs are known to be involved in apoptotic (programmed cell death) pathway, and biotic and abiotic stress responses in plants. Here, the members of LOX gene family in Arabidopsis and rice are identified. The Arabidopsis and rice genomes encode 6 and 14 LOX proteins, respectively, and interestingly, with more LOX genes in rice. The rice LOXs are validated based on protein alignment studies. This is the first report wherein LOXs are identified in rice which may allow better understanding the initiation, progression and effects of apoptosis, and responses to bitoic and abiotic stresses and signaling cascades in plants.Key words: apoptosis, biotic and abiotic stresses, genomics, jasmonic acid, lipidsLipoxygenases (linoleate:oxygen oxidoreductase, EC 1.13.11.-; LOXs) catalyze the conversion of polyunsaturated fatty acids (lipids) into conjugated hydroperoxides. This process is called hydroperoxidation of lipids. LOXs are monomeric, non-heme and non-sulfur, but iron-containing dioxygenases widely expressed in fungi, animal and plant cells, and are known to be absent in prokaryotes. However, a recent finding suggests the existence of LOX-related genomic sequences in bacteria but not in archaea.¹ The inflammatory conditions in mammals like bronchial asthama, psoriasis and arthritis are a result of LOXs reactions.² Further, several clinical conditions like HIV-1 infection,³ disease of kidneys due to the activation of 5-lipoxygenase,^4,5 aging of the brain due to neuronal 5-lipoxygenase⁶ and atherosclerosis⁷ are mediated by LOXs. In plants, LOXs are involved in response to biotic and abiotic stresses.⁸ They are involved in germination⁹ and also in traumatin and jasmonic acid biochemical pathways.^10,11 Studies on LOX in rice are conducted to develop novel strategies against insect pests¹² in response to wounding and insect attack,¹³ and on rice bran extracts as functional foods and dietary supplements for control of inflammation and joint health.¹⁴ In Arabidopsis, LOXs are studied in response to natural and stress-induced senescence,¹⁵ transition to flowering,¹⁶ regulation of lateral root development and defense response.¹⁷The arachidonic, linoleic and linolenic acids can act as substrates for different LOX isozymes. A hydroperoxy group is added at carbons 5, 12 or 15, when arachidonic acid is the substrate, and so the LOXs are designated as 5-, 12- or 15-lipoxygenases. Sequences are available in the database for plant lipoxygenases (EC:1.13.11.12), mammalian arachidonate 5-lipoxygenase (EC:1.13.11.34), mammalian arachidonate 12-lipoxygenase (EC:1.13.11.31) and mammalian erythroid cell-specific 15-lipoxygenase (EC:1.13.11.33). The prototype member for LOX family, LOX-1 of Glycine max L. (soybean) is a 15-lipoxygenase. The LOX isoforms of soybean (LOX-1, LOX-2, LOX-3a and LOX-3b) are the most characterized of plant LOXs.¹⁸ In addition, five vegetative LOXs (VLX-A, -B, -C, -D, -E) are detected in soybean leaves.¹⁹ The 3-dimensional structure of soybean LOX-1 has been determined.^20,21 LOX-1 was shown to be made of two domains, the N-terminal domain-I which forms a β-barrel of 146 residues, and a C-terminal domain-II of bundle of helices of 693 residues²¹ (Fig. 1). The iron atom was shown to be at the centre of domain-II bound by four coordinating ligands, of which three are histidine residues.²²Open in a separate windowFigure 1Three-dimensional structure of soybean lipoxygenase L-1. The domain I (N-terminal) and domain II (C-terminal) are indicated. The catalytic iron atom is embedded in domain II (PDB ID-1YGE).²¹This article describes identification of LOX genes in Arabidopsis and rice. The Arabidopsis genome encodes for six LOX proteins²³ (www.arabidopsis.org) (

小拟南芥MKK基因家族全基因组鉴定及进化和表达分析

丝裂原活化蛋白激酶激酶(mitogen-activated protein kinase kinase,MAPKK或MKK)是丝裂原活化蛋白激酶(mitogen-activatedproteinkinase,MAPK)级联的重要组成部分,在植物的生长发育和胁迫应答过程中发挥重要作用。目前,已在多种植物中鉴定了MKK基因家族,但在十字花科植物小拟南芥(Arabidopsis pumila)中MKK基因家族的系统鉴定与分析尚未见报道。为了探索小拟南芥MKK基因家族的进化和功能,本研究通过全基因组分析鉴定了小拟南芥中16个MKK基因,散布于小拟南芥的10条染色体上。基于系统发育分析和多重序列比对,将这些基因分为5个亚族:A亚族(5个)、B亚族(2个)、C亚族(4个)、D亚族(3个)和E亚族(2个)。分子进化和共线性分析表明小拟南芥中存在7对复制基因,分别是ApMKK1-1/1-2、ApMKK2-1/2-2、ApMKK3-1/3-2、ApMKK4-1/4-2、ApMKK5-1/5-2、ApMKK9-1/9-2和ApMKK10-1/10-2,其中ApMKK1-1/1-2在复制事件之后发生了加速进化。结合ApMKKs启动子区的顺式元件分布和ApMKKs在成熟叶片、茎、花和果实以及盐胁迫下的表达模式,结果发现复制基因的表达具有组织特异性和功能多样性。部分复制基因在组织中的表达模式存在差异,但在盐胁迫下的表达模式却基本相同。本研究结果为解析MKK介导的小拟南芥发育过程和非生物胁迫信号转导通路的复杂机制奠定了基础。     

拟南芥和水稻SET结构域基因家族全基因组鉴定、分类和表达

ET(Su(var), Enhancer of zeste (E(z)), and Trithorax)结构域基因家族是一组含有保守SET结构域的蛋白的统称, 它们参与蛋白甲基化, 影响染色体结构, 并且调控基因表达, 在植物发育中起着重要的作用。分析拟南芥和水稻中SET结构域基因家族进化关系, 对研究这一基因家族中各成员的功能有着重要的意义。我们系统地鉴定了47个拟南芥(Arabidopsis thaliana)和43个水稻(Orysa sativa japonica cultivar Nipponbare)的SET结构域基因, 染色体定位和基因复制分析表明SET结构域基因扩增是由片段复制和反转录引起的, 根据这些结构域差异和系统发育分析把拟南芥和水稻的SET结构域基因划分成5个亚家族。通过分析SET结构域基因家族在拟南芥和水稻各个发育阶段的表达谱, 发现SET结构域基因绝大部分至少在一个组织中表达; 大部分在花和花粉中高表达; 一些SET结构域基因在某些组织中有特异的表达模式, 表明与组织发育有密切的关系。在拟南芥和水稻中分别找到了4个差异表达基因。拟南芥4个差异基因都在花粉管高表达, 水稻4个差异基因有3个在雄性花蕊中高表达, 另一个在幼穗中高表达。     

Genome-wide crossover distribution in Arabidopsis thaliana meiosis reveals sex-specific patterns along chromosomes

In most species, crossovers (COs) are essential for the accurate segregation of homologous chromosomes at the first meiotic division. Their number and location are tightly regulated. Here, we report a detailed, genome-wide characterization of the rate and localization of COs in Arabidopsis thaliana, in male and female meiosis. We observed dramatic differences between male and female meiosis which included: (i) genetic map length; 575 cM versus 332 cM respectively; (ii) CO distribution patterns: male CO rates were very high at both ends of each chromosome, whereas female CO rates were very low; (iii) correlations between CO rates and various chromosome features: female CO rates correlated strongly and negatively with GC content and gene density but positively with transposable elements (TEs) density, whereas male CO rates correlated positively with the CpG ratio. However, except for CpG, the correlations could be explained by the unequal repartition of these sequences along the Arabidopsis chromosome. For both male and female meiosis, the number of COs per chromosome correlates with chromosome size expressed either in base pairs or as synaptonemal complex length. Finally, we show that interference modulates the CO distribution both in male and female meiosis.