Berberine enhances defects in the establishment of leaf polarity in asymmetric leaves1 and asymmetric leaves2 of Arabidopsis thaliana

A Common Molecular Signature Indicates the Pre-Meristematic State of Plant Calli

In response to different degrees of mechanical injury, certain plant cells re-enter the division cycle to provide cells for tissue replenishment, tissue rejoining, de novo organ formation, and/or wound healing. The intermediate tissue formed by the dividing cells is called a callus. Callus formation can also be induced artificially in vitro by wounding and/or hormone (auxin and cytokinin) treatments. The callus tissue can be maintained in culture, providing starting material for de novo organ or embryo regeneration and thus serving as the basis for many plant biotechnology applications. Due to the biotechnological importance of callus cultures and the scientific interest in the developmental flexibility of somatic plant cells, the initial molecular steps of callus formation have been studied in detail. It was revealed that callus initiation can follow various ways, depending on the organ from which it develops and the inducer, but they converge on a seemingly identical tissue. It is not known, however, if callus is indeed a special tissue with a defined gene expression signature, whether it is a malformed meristem, or a mass of so-called "undifferentiated" cells, as is mostly believed. In this paper, I review the various mechanisms of plant regeneration that may converge on callus initiation. I discuss the role of plant hormones in the detour of callus formation from normal development. Finally, I compare various Arabidopsis gene expression datasets obtained a few days, two weeks, or several years after callus induction and identify 21 genes, including genes of key transcription factors controlling cell division and differentiation in meristematic regions, which were upregulated in all investigated callus samples. I summarize the information available on all 21 genes that point to the pre-meristematic nature of callus tissues underlying their wide regeneration potential.

Arabidopsis ASYMMETRIC LEAVES2 (AS2): roles in plant morphogenesis, cell division, and pathogenesis

The ASYMMETRIC LEAVES2 (AS2) gene in Arabidopsis thaliana is responsible for the development of flat, symmetric, and extended leaf laminae and their vein systems. AS2 protein is a member of the plant-specific AS2/LOB protein family, which includes 42 members comprising the conserved amino-terminal domain referred to as the AS2/LOB domain, and the variable carboxyl-terminal region. Among the members, AS2 has been most intensively investigated on both genetic and molecular levels. AS2 forms a complex with the myb protein AS1, and is involved in epigenetic repression of the abaxial genes ETTIN/AUXIN RESPONSE FACTOR3 (ETT/ARF3), ARF4, and class 1 KNOX homeobox genes. The repressed expression of these genes by AS2 is markedly enhanced by the cooperative action of various modifier genes, some of which encode nucleolar proteins. Further downstream, progression of the cell division cycle in the developing organs is stimulated; meristematic states are suppressed in determinate leaf primordia; and the extension of leaf primordia is induced. AS2 binds the specific sequence in exon 1 of ETT/ARF3 and maintains methylated CpGs in several exons of ETT/ARF3. AS2 forms bodies (designated as AS2 bodies) at nucleolar peripheries. AS2 bodies partially overlap chromocenters, including inactive 45S ribosomal DNA repeats, suggesting the presence of molecular and functional links among AS2, the 45S rDNAs, and the nucleolus to exert the repressive regulation of ETT/ARF3. The AS2/LOB domain is characterized by three subdomains, the zinc finger (ZF) motif, the internally conserved-glycine containing (ICG) region, and the leucine-zipper-like (LZL) region. Each of these subdomains is essential for the formation of AS2 bodies. ICG to LZL are required for nuclear localization, but ZF is not. LZL intrinsically has the potential to be exported to the cytoplasm. In addition to its nuclear function, it has been reported that AS2 plays a positive role in geminivirus infection: its protein BV1 stimulates the expression of AS2 and recruits AS2 to the cytoplasm, which enhances virus infectivity by suppression of cytoplasmic post transcriptional gene silencing.

Cytotoxicity of the natural herbicidal chemical, berberine, on Nicotiana tabacum Bright yellow-2 cells

Berberine is a naturally occurring plant secondary metabolite with allelopathic and cytotoxic properties. We investigated the cytotoxic effects of berberine against tobacco Bright Yellow-2 (BY-2) cells. We found that berberine inhibited tobacco BY-2 cell growth in a concentration-dependent manner and the potency of berberine was comparable to the traditional herbicide glyphosate. Meanwhile, the in vivo test revealed that the herbicidal activity of berberine was also comparable to that of glyphosate. Further mechanism studies for the cytotoxicity demonstrated that berberine at concentrations of 40 μg/mL and 80 μg/mL induced cell death by causing mitochondrial membrane depolarization, irregular nuclei and chromatin condensation but not leading to significant apoptosis. Ultra-structure analysis through transmission electron microscopy (TEM) indicated that treatments with 40 μg/mL berberine for 2 d or 80 μg/mL berberine for 1 d induced cell damage, including nuclei morphological changes, cytoplasm and mitochondria degradation and cell wall fibrosis. Treatment at higher concentration of 80 μg/mL for longer period of 2 d, induced plasmolysis and severe damage to basic cell structure, which are indicative of explosive necrosis. Our results suggest that berberine was cytotoxic to tobacco BY-2 cells by inducing necrosis but not apoptosis. Cell wall, mitochondria, nuclei and chromatin were damaged and may be possible primary targets. Therefore, the herbicidal activity of berberine appears to be a complex process associated with multiple mechanisms of action.

A genetic link between epigenetic repressor AS1-AS2 and DNA replication factors in establishment of adaxial-abaxial leaf polarity of Arabidopsis

Balanced development of adaxial and abaxial domains in leaf primordia is critical for the formation of flat symmetric leaf lamina. Arabidopsis ASYMMETRIC LEAVES1 (AS1) and AS2 proteins form a complex (AS1-AS2), which acts as key regulators for the adaxial development by the direct repression of expression of the abaxial gene ETTIN/AUXIN RESPONSE FACTOR3 (ETT/ARF3). Many modifier mutations have been identified, which enhance the defect of as1 and as2 mutations to generate abaxialized filamentous leaves without adaxial traits, suggesting that the development of the adaxial domain is achieved by cooperative repression by AS1-AS2 and the wild-type proteins corresponding to the modifiers. Mutations of several genes for DNA replication-related chromatin remodeling factors such as Chromatin Assembly Factor-1 (CAF-1) have been also identified as modifiers. It is still unknown, however, whether mutations in genes involved in DNA replication themselves might act as modifiers. Here we report that as1 and as2 mutants grown in the presence of hydroxyurea, a known inhibitor of DNA replication, form abaxialized filamentous leaves in a concentration-dependent manner. We further show that a mutation of the INCURVATA2 (ICU2) gene, which encodes the putative catalytic subunit of DNA polymerase α, and a mutation of the Replication Factor C Subunit3 (RFC3) gene, which encodes a protein used in replication as a clamp loader, act as modifiers. In addition, as2-1 icu2-1 double mutants showed increased mRNA levels of the genes for leaf abaxialization. These results suggest a tight link between DNA replication and the function of AS1-AS2 in the development of flat leaves.

The Allelochemical MDCA Inhibits Lignification and Affects Auxin Homeostasis

The phenylpropanoid 3,4-(methylenedioxy)cinnamic acid (MDCA) is a plant-derived compound first extracted from roots of Asparagus officinalis and further characterized as an allelochemical. Later on, MDCA was identified as an efficient inhibitor of 4-COUMARATE-CoA LIGASE (4CL), a key enzyme of the general phenylpropanoid pathway. By blocking 4CL, MDCA affects the biosynthesis of many important metabolites, which might explain its phytotoxicity. To decipher the molecular basis of the allelochemical activity of MDCA, we evaluated the effect of this compound on Arabidopsis thaliana seedlings. Metabolic profiling revealed that MDCA is converted in planta into piperonylic acid (PA), an inhibitor of CINNAMATE-4-HYDROXYLASE (C4H), the enzyme directly upstream of 4CL. The inhibition of C4H was also reflected in the phenolic profile of MDCA-treated plants. Treatment of in vitro grown plants resulted in an inhibition of primary root growth and a proliferation of lateral and adventitious roots. These observed growth defects were not the consequence of lignin perturbation, but rather the result of disturbing auxin homeostasis. Based on DII-VENUS quantification and direct measurement of cellular auxin transport, we concluded that MDCA disturbs auxin gradients by interfering with auxin efflux. In addition, mass spectrometry was used to show that MDCA triggers auxin biosynthesis, conjugation, and catabolism. A similar shift in auxin homeostasis was found in the c4h mutant ref3-2, indicating that MDCA triggers a cross talk between the phenylpropanoid and auxin biosynthetic pathways independent from the observed auxin efflux inhibition. Altogether, our data provide, to our knowledge, a novel molecular explanation for the phytotoxic properties of MDCA.

Analysis of gene expression profiles of soft tissue sarcoma using a combination of knowledge-based filtering with integration of multiple statistics

The diagnosis and treatment of soft tissue sarcomas (STS) have been difficult. Of the diverse histological subtypes, undifferentiated pleomorphic sarcoma (UPS) is particularly difficult to diagnose accurately, and its classification per se is still controversial. Recent advances in genomic technologies provide an excellent way to address such problems. However, it is often difficult, if not impossible, to identify definitive disease-associated genes using genome-wide analysis alone, primarily because of multiple testing problems. In the present study, we analyzed microarray data from 88 STS patients using a combination method that used knowledge-based filtering and a simulation based on the integration of multiple statistics to reduce multiple testing problems. We identified 25 genes, including hypoxia-related genes (e.g., MIF, SCD1, P4HA1, ENO1, and STAT1) and cell cycle- and DNA repair-related genes (e.g., TACC3, PRDX1, PRKDC, and H2AFY). These genes showed significant differential expression among histological subtypes, including UPS, and showed associations with overall survival. STAT1 showed a strong association with overall survival in UPS patients (logrank p = 1.84 × 10(-6) and adjusted p value 2.99 × 10(-3) after the permutation test). According to the literature, the 25 genes selected are useful not only as markers of differential diagnosis but also as prognostic/predictive markers and/or therapeutic targets for STS. Our combination method can identify genes that are potential prognostic/predictive factors and/or therapeutic targets in STS and possibly in other cancers. These disease-associated genes deserve further preclinical and clinical validation.

Application of a combination of a knowledge-based algorithm and 2-stage screening to hypothesis-free genomic data on irinotecan-treated patients for identification of a candidate single nucleotide polymorphism related to an adverse effect

Interindividual variation in a drug response among patients is known to cause serious problems in medicine. Genomic information has been proposed as the basis for “personalized” health care. The genome-wide association study (GWAS) is a powerful technique for examining single nucleotide polymorphisms (SNPs) and their relationship with drug response variation; however, when using only GWAS, it often happens that no useful SNPs are identified due to multiple testing problems. Therefore, in a previous study, we proposed a combined method consisting of a knowledge-based algorithm, 2 stages of screening, and a permutation test for identifying SNPs. In the present study, we applied this method to a pharmacogenomics study where 109,365 SNPs were genotyped using Illumina Human-1 BeadChip in 168 cancer patients treated with irinotecan chemotherapy. We identified the SNP rs9351963 in potassium voltage-gated channel subfamily KQT member 5 (KCNQ5) as a candidate factor related to incidence of irinotecan-induced diarrhea. The p value for rs9351963 was 3.31×10⁻⁵ in Fisher's exact test and 0.0289 in the permutation test (when multiple testing problems were corrected). Additionally, rs9351963 was clearly superior to the clinical parameters and the model involving rs9351963 showed sensitivity of 77.8% and specificity of 57.6% in the evaluation by means of logistic regression. Recent studies showed that KCNQ4 and KCNQ5 genes encode members of the M channel expressed in gastrointestinal smooth muscle and suggested that these genes are associated with irritable bowel syndrome and similar peristalsis diseases. These results suggest that rs9351963 in KCNQ5 is a possible predictive factor of incidence of diarrhea in cancer patients treated with irinotecan chemotherapy and for selecting chemotherapy regimens, such as irinotecan alone or a combination of irinotecan with a KCNQ5 opener. Nonetheless, clinical importance of rs9351963 should be further elucidated.

Regulation of compound leaf development by PHANTASTICA in Medicago truncatula

Plant leaves, simple or compound, initiate as peg-like structures from the peripheral zone of the shoot apical meristem, which requires class I KNOTTED-LIKE HOMEOBOXI (KNOXI) transcription factors to maintain its activity. The MYB domain protein encoded by the ASYMMETRIC LEAVES1/ROUGH SHEATH2/PHANTASTICA (ARP) gene, together with other factors, excludes KNOXI gene expression from incipient leaf primordia to initiate leaves and specify leaf adaxial identity. However, the regulatory relationship between ARP and KNOXI is more complex in compound-leafed species. Here, we investigated the role of ARP and KNOXI genes in compound leaf development in Medicago truncatula. We show that the M. truncatula phantastica mutant exhibited severe compound leaf defects, including curling and deep serration of leaf margins, shortened petioles, increased rachises, petioles acquiring motor organ characteristics, and ectopic development of petiolules. On the other hand, the M. truncatula brevipedicellus mutant did not exhibit visible compound leaf defects. Our analyses show that the altered petiole development requires ectopic expression of ELONGATED PETIOLULE1, which encodes a lateral organ boundary domain protein, and that the distal margin serration requires the auxin efflux protein M. truncatula PIN-FORMED10 in the M. truncatula phantastica mutant.

Macrophage migration inhibitory factor and stearoyl-CoA desaturase 1: potential prognostic markers for soft tissue sarcomas based on bioinformatics analyses

The diagnosis and treatment of soft tissue sarcomas (STSs) has been particularly difficult, because STSs are a group of highly heterogeneous tumors in terms of histopathology, histological grade, and primary site. Recent advances in genome technologies have provided an excellent opportunity to determine the complete biological characteristics of neoplastic tissues, resulting in improved diagnosis, treatment selection, and investigation of therapeutic targets. We had previously developed a novel bioinformatics method for marker gene selection and applied this method to gene expression data from STS patients. This previous analysis revealed that the extracted gene combination of macrophage migration inhibitory factor (MIF) and stearoyl-CoA desaturase 1 (SCD1) is an effective diagnostic marker to discriminate between subtypes of STSs with highly different outcomes. In the present study, we hypothesize that the combination of MIF and SCD1 is also a prognostic marker for the overall outcome of STSs. To prove this hypothesis, we first analyzed microarray data from 88 STS patients and their outcomes. Our results show that the survival rates for MIF- and SCD1-positive groups were lower than those for negative groups, and the p values of the log-rank test are 0.0146 and 0.00606, respectively. In addition, survival rates are more significantly different (p = 0.000116) between groups that are double-positive and double-negative for MIF and SCD1. Furthermore, in vitro cell growth inhibition experiments by MIF and SCD1 inhibitors support the hypothesis. These results suggest that the gene set is useful as a prognostic marker associated with tumor progression.

Identification of a candidate single-nucleotide polymorphism related to chemotherapeutic response through a combination of knowledge-based algorithm and hypothesis-free genomic data

Inter-individual variations in drug responses among patients are known to cause serious problems in medicine. Genome-wide association study (GWAS) is powerful for examining single-nucleotide polymorphisms (SNPs) and their relationships with drug response variations. However, no significant SNP has been identified using GWAS due to multiple testing problems. Therefore, we propose a combination method consisting of knowledge-based algorithm, two stages of screening, and permutation test for identifying SNPs in the present study. We applied this method to a genome-wide pharmacogenomics study for which 109,365 SNPs had been genotyped using Illumina Human-1 BeadChip for 119 gastric cancer patients treated with fluoropyrimidine. We identified rs2293347 in epidermal growth factor receptor (EGFR) is as a candidate SNP related to chemotherapeutic response. The p value for the rs2293347 was 2.19 × 10(-5) for Fisher's exact test, and the p value was 0.00360 for the permutation test (multiple testing problems are corrected). Additionally, rs2293347 was clearly superior to clinical parameters and showed a sensitivity value of 55.0% and specificity value of 94.4% in the evaluation by using multiple regression models. Recent studies have shown that combination chemotherapy of fluoropyrimidine and EGFR-targeting agents is effective for gastric cancer patients highly expressing EGFR. These results suggest that rs2293347 is a potential predictive factor for selecting chemotherapies, such as fluoropyrimidine alone or combination chemotherapies.

Dual regulation of ETTIN (ARF3) gene expression by AS1-AS2, which maintains the DNA methylation level, is involved in stabilization of leaf adaxial-abaxial partitioning in Arabidopsis

Leaf primordia are generated at the periphery of the shoot apex, developing into flat symmetric organs with adaxial-abaxial polarity, in which the indeterminate state is repressed. Despite the crucial role of the ASYMMETRIC LEAVES1 (AS1)-AS2 nuclear-protein complex in leaf adaxial-abaxial polarity specification, information on mechanisms controlling their downstream genes has remained elusive. We systematically analyzed transcripts by microarray and chromatin immunoprecipitation assays and performed genetic rescue of as1 and as2 phenotypic abnormalities, which identified a new target gene, ETTIN (ETT)/AUXIN RESPONSE FACTOR3 (ARF3), which encodes an abaxial factor acting downstream of the AS1-AS2 complex. While the AS1-AS2 complex represses ETT by direct binding of AS1 to the ETT promoter, it also indirectly activates miR390- and RDR6-dependent post-transcriptional gene silencing to negatively regulate both ETT and ARF4 activities. Furthermore, AS1-AS2 maintains the status of DNA methylation in the ETT coding region. In agreement, filamentous leaves formed in as1 and as2 plants treated with a DNA methylation inhibitor were rescued by loss of ETT and ARF4 activities. We suggest that negative transcriptional, post-transcriptional and epigenetic regulation of the ARFs by AS1-AS2 is important for stabilizing early leaf partitioning into abaxial and adaxial domains.

Abaxial Greening Phenotype in Hybrid Aspen

The typical angiosperm leaf, as in Arabidopsis, is bifacial consisting of top (adaxial) and bottom (abaxial) surfaces readily distinguishable by the underlying cell type (palisade and spongy mesophyll, respectively). Species of the genus Populus have leaves that are either conventionally bifacial or isobilateral. Isobilateral leaves have palisade mesophyll on the top and bottom of the leaf, making the two sides virtually indistinguishable at the macroscopic level. In poplars this has been termed the “abaxial greening” phenotype. Previous work has implicated ASYMMETRIC LEAVES1 (AS1) as an essential determinant of palisade mesophyll development. This gene, as well as other genes (84 in all) putatively involved in setting the dorsiventral axis of leaves, were investigated in two Populus species: black cottonwood (Populus trichocarpa) and hybrid aspen (P. tremula x tremuloides), representative of each leaf type (bifacial and isobilateral, respectively). Poplar orthologs of AS1 have significantly higher expression in aspen leaf blade and lower in the petiole, suggestive of a potential role in the isobilateral leaf phenotype consistent with the previously observed phenotypes. Furthermore, an ABERRANT TESTA SHAPE (ATS) ortholog has significantly lower expression in aspen leaf tissue, also suggesting a possible contribution of this gene to abaxial greening.

Meta-analyses of microarrays of Arabidopsis asymmetric leaves1 (as1), as2 and their modifying mutants reveal a critical role for the ETT pathway in stabilization of adaxial-abaxial patterning and cell division during leaf development

It is necessary to use algorithms to analyze gene expression data from DNA microarrays, such as in clustering and machine learning. Previously, we developed the knowledge-based fuzzy adaptive resonance theory (KB-FuzzyART), a clustering algorithm suitable for analyzing gene expression data, to find clues for identifying gene networks. Leaf primordia form around the shoot apical meristem (SAM), which consists of indeterminate stem cells. Upon initiation of leaf development, adaxial-abaxial patterning is crucial for lateral expansion, via cellular proliferation, and the formation of flat symmetric leaves. Many regulatory genes that specify such patterning have been identified. Analysis by the KB-FuzzyART and subsequent molecular and genetic analyses previously showed that ASYMMETRIC LEAVES1 (AS1) and AS2 repress the expression of some abaxial-determinant genes, such as AUXIN RESPONSE FACTOR3 (ARF3)/ETTIN (ETT) and ARF4, which are responsible for defects in leaf adaxial-abaxial polarity in as1 and as2. In the present study, genetic analysis revealed that ARF3/ETT and ARF4 were regulated by modifier genes, BOBBER1 (BOB1) and ELONGATA3 (ELO3), together with AS1-AS2. We analyzed expression arrays with as2 elo3 and as2 bob1, and extracted genes downstream of ARF3/ETT by using KB-FuzzyART and molecular analyses. The results showed that expression of Kip-related protein (KRP) (for inhibitors of cyclin-dependent protein kinases) and Isopentenyltransferase (IPT) (for biosynthesis of cytokinin) genes were controlled by AS1-AS2 through ARF3/ETT and ARF4 functions, which suggests that the AS1-AS2-ETT pathway plays a critical role in controlling the cell division cycle and the biosynthesis of cytokinin around SAM to stabilize leaf development in Arabidopsis thaliana.

Leaf development

Leaves are the most important organs for plants. Without leaves, plants cannot capture light energy or synthesize organic compounds via photosynthesis. Without leaves, plants would be unable perceive diverse environmental conditions, particularly those relating to light quality/quantity. Without leaves, plants would not be able to flower because all floral organs are modified leaves. Arabidopsis thaliana is a good model system for analyzing mechanisms of eudicotyledonous, simple-leaf development. The first section of this review provides a brief history of studies on development in Arabidopsis leaves. This history largely coincides with a general history of advancement in understanding of the genetic mechanisms operating during simple-leaf development in angiosperms. In the second section, I outline events in Arabidopsis leaf development, with emphasis on genetic controls. Current knowledge of six important components in these developmental events is summarized in detail, followed by concluding remarks and perspectives.

Discovery of novel rules for G-quadruplex-forming sequences in plants by using bioinformatics methods

The G-quadruplex is one of the most frequently studied secondary DNA structures and consists of 4 guanine residues that interact through Watson-Crick and Hoogsteen pairing. The G-quadruplex formation is thought to be a molecular switch for gene expression. Genome-wide analyses of G-quadruplexes have been published for many species; however, only one genome-wide analysis of G-quadruplexes in plants has been reported. Here, we propose a new approach involving a two-step procedure for identifying G-quadruplex-forming sequences (potential G4 DNA motif regions: G4MRs) and classifying positional relationships between G4MRs and genes. By using this approach, we exhaustively searched for G4MRs in the whole genomes of 8 species: Arabidopsis thaliana, Oryza sativa subsp. japonica, Populus trichocarpa, Vitis vinifera, Homo sapiens, Danio rerio, Drosophila melanogaster, and Caenorhabditis elegans. We classified genes on the basis of their positional relationships to their proximal G4MRs. We identified novel rules for G4MRs in plants, such as G4MR-enrichment in the template strands at transcription start sites (TSSs). Next, we focused on the template strands of TSSs and conducted gene ontology (GO) analysis of genes proximal to G4MRs. We identified GO terms such as chloroplast and nucleosome (or histone) in O. sativa. Although these terms were strongly associated in O. sativa, weak associations were identified in other plants. These results will be helpful for elucidating the functional roles of G4 DNA.