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

The Epigenomic Features and Potential Functions of PEG- and PDS-Favorable DNA G-Quadruplexes in Rice

A G-quadruplex (G4) is a typical non-B DNA structure and involved in various DNA-templated events in eukaryotic genomes. PEG and PDS chemicals have been widely applied for promoting the folding of in vivo or in vitro G4s. However, how PEG and PDS preferentially affect a subset of G4 formation genome-wide is still largely unknown. We here conducted a BG4-based IP-seq in vitro under K++PEG or K++PDS conditions in the rice genome. We found that PEG-favored IP-G4s+ have distinct sequence features, distinct genomic distributions and distinct associations with TEGs, non-TEGs and subtypes of TEs compared to PDS-favored ones. Strikingly, PEG-specific IP-G4s+ are associated with euchromatin with less enrichment levels of DNA methylation but with more enriched active histone marks, while PDS-specific IP-G4s+ are associated with heterochromatin with higher enrichment levels of DNA methylation and repressive marks. Moreover, we found that genes with PEG-specific IP-G4s+ are more expressed than those with PDS-specific IP-G4s+, suggesting that PEG/PDS-specific IP-G4s+ alone or coordinating with epigenetic marks are involved in the regulation of the differential expression of related genes, therefore functioning in distinct biological processes. Thus, our study provides new insights into differential impacts of PEG and PDS on G4 formation, thereby advancing our understanding of G4 biology.

Beyond the Primary Structure of Nucleic Acids: Potential Roles of Epigenetics and Noncanonical Structures in the Regulations of Plant Growth and Stress Responses

Epigenetics deals with changes in gene expression that are not caused by modifications in the primary sequence of nucleic acids. These changes beyond primary structures of nucleic acids not only include DNA/RNA methylation, but also other reversible conversions, together with histone modifications or RNA interference. In addition, under particular conditions (such as specific ion concentrations or protein-induced stabilization), the right-handed double-stranded DNA helix (B-DNA) can form noncanonical structures commonly described as "non-B DNA" structures. These structures comprise, for example, cruciforms, i-motifs, triplexes, and G-quadruplexes. Their formation often leads to significant differences in replication and transcription rates. Noncanonical RNA structures have also been documented to play important roles in translation regulation and the biology of noncoding RNAs. In human and animal studies, the frequency and dynamics of noncanonical DNA and RNA structures are intensively investigated, especially in the field of cancer research and neurodegenerative diseases. In contrast, noncanonical DNA and RNA structures in plants have been on the fringes of interest for a long time and only a few studies deal with their formation, regulation, and physiological importance for plant stress responses. Herein, we present a review focused on the main fields of epigenetics in plants and their possible roles in stress responses and signaling, with special attention dedicated to noncanonical DNA and RNA structures.

Exploring the Relationship between G-Quadruplex Nucleic Acids and Plants: From Plant G-Quadruplex Function to Phytochemical G4 Ligands with Pharmaceutic Potential

G-quadruplex (G4) oligonucleotides are higher-order DNA and RNA secondary structures of enormous relevance due to their implication in several biological processes and pathological states in different organisms. Strategies aiming at modulating human G4 structures and their interrelated functions are first-line approaches in modern research aiming at finding new potential anticancer treatments or G4-based aptamers for various biomedical and biotechnological applications. Plants offer a cornucopia of phytocompounds that, in many cases, are effective in binding and modulating the thermal stability of G4s and, on the other hand, contain almost unexplored G4 motifs in their genome that could inspire new biotechnological strategies. Herein, we describe some G4 structures found in plants, summarizing the existing knowledge of their functions and biological role. Moreover, we review some of the most promising G4 ligands isolated from vegetal sources and report on the known relationships between such phytochemicals and G4-mediated biological processes that make them potential leads in the pharmaceutical sector.

The Newly Sequenced Genome of Pisum sativum Is Replete with Potential G-Quadruplex-Forming Sequences-Implications for Evolution and Biological Regulation

G-quadruplexes (G4s) have been long considered rare and physiologically unimportant in vitro curiosities, but recent methodological advances have proved their presence and functions in vivo. Moreover, in addition to their functional relevance in bacteria and animals, including humans, their importance has been recently demonstrated in evolutionarily distinct plant species. In this study, we analyzed the genome of Pisum sativum (garden pea, or the so-called green pea), a unique member of the Fabaceae family. Our results showed that this genome contained putative G4 sequences (PQSs). Interestingly, these PQSs were located nonrandomly in the nuclear genome. We also found PQSs in mitochondrial (mt) and chloroplast (cp) DNA, and we experimentally confirmed G4 formation for sequences found in these two organelles. The frequency of PQSs for nuclear DNA was 0.42 PQSs per thousand base pairs (kbp), in the same range as for cpDNA (0.53/kbp), but significantly lower than what was found for mitochondrial DNA (1.58/kbp). In the nuclear genome, PQSs were mainly associated with regulatory regions, including 5'UTRs, and upstream of the rRNA region. In contrast to genomic DNA, PQSs were located around RNA genes in cpDNA and mtDNA. Interestingly, PQSs were also associated with specific transposable elements such as TIR and LTR and around them, pointing to their role in their spreading in nuclear DNA. The nonrandom localization of PQSs uncovered their evolutionary and functional significance in the Pisum sativum genome.

Epigenomic features of DNA G-quadruplexes and their roles in regulating rice gene transcription

A DNA G-quadruplex (G4) is a non-canonical four-stranded nucleic acid structure involved in many biological processes in mammals. The current knowledge on plant DNA G4s, however, is limited; whether and how DNA G4s impact gene expression in plants is still largely unknown. Here, we applied a protocol referred to as BG4-DNA-IP-seq followed by a comprehensive characterization of DNA G4s in rice (Oryza sativa L.); we next integrated dG4s (experimentally detectable G4s) with existing omics data and found that dG4s exhibited differential DNA methylation between transposable element (TE) and non-TE genes. dG4 regions displayed genic-dependent enrichment of epigenomic signatures; finally, we showed that these sites displayed a positive association with expression of DNA G4-containing genes when located at promoters, and a negative association when located in the gene body, suggesting localization-dependent promotional/repressive roles of DNA G4s in regulating gene transcription. This study reveals interrelations between DNA G4s and epigenomic signatures, as well as implicates DNA G4s in modulating gene transcription in rice. Our study provides valuable resources for the functional characterization or bioengineering of some of key DNA G4s in rice.

Genome-wide discovery of G-quadruplexes in barley

G-quadruplexes (G4s) are four-stranded nucleic acid structures with closely spaced guanine bases forming square planar G-quartets. Aberrant formation of G4 structures has been associated with genomic instability. However, most plant species are lacking comprehensive studies of G4 motifs. In this study, genome-wide identification of G4 motifs in barley was performed, followed by a comparison of genomic distribution and molecular functions to other monocot species, such as wheat, maize, and rice. Similar to the reports on human and some plants like wheat, G4 motifs peaked around the 5′ untranslated region (5′ UTR), the first coding domain sequence, and the first intron start sites on antisense strands. Our comparative analyses in human, Arabidopsis, maize, rice, and sorghum demonstrated that the peak points could be erroneously merged into a single peak when large window sizes are used. We also showed that the G4 distributions around genic regions are relatively similar in the species studied, except in the case of Arabidopsis. G4 containing genes in monocots showed conserved molecular functions for transcription initiation and hydrolase activity. Additionally, we provided examples of imperfect G4 motifs.

Genome-Wide Discovery of G-Quadruplexes in Wheat: Distribution and Putative Functional Roles

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world's most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

Excess sterols disrupt plant cellular activity by inducing stress-responsive gene expression

Sterols are important lipid constituents of cellular membranes in plants and other organisms. Sterol homeostasis is under strict regulation in plants because excess sterols negatively impact plant growth. HIGH STEROL ESTER 1 (HISE1) functions as a negative regulator of sterol accumulation. If sterol production exceeds a certain threshold, excess sterols are detoxified via conversion to sterol esters by PHOSPHOLIPID STEROL ACYL TRANSFERASE 1 (PSAT1). We previously reported that the Arabidopsis thaliana double mutant hise1-3 psat1-2 shows 1.5-fold higher sterol content than the wild type and consequently a severe growth defect. However, the specific defects caused by excess sterol accumulation in plants remain unknown. In this study, we investigated the effects of excess sterols on plants by analyzing the phenotypes and transcriptomes of the hise1-3 psat1-2 double mutant. Transcriptomic analysis revealed that 435 genes were up-regulated in hise1-3 psat1-2 leaves compared with wild-type leaves. Gene ontology (GO) enrichment analysis revealed that abiotic and biotic stress-responsive genes including RESPONSIVE TO DESICCATION 29B/LOW-TEMPERATURE-INDUCED 65 (RD29B/LTI65) and COLD-REGULATED 15A (COR15A) were up-regulated in hise1-3 psat1-2 leaves compared with wild-type leaves. Expression levels of senescence-related genes were also much higher in hise1-3 psat1-2 leaves than in wild-type leaves. hise1-3 psat1-2 leaves showed early senescence, suggesting that excess sterols induce senescence of leaves. In the absence of sucrose, hise1-3 psat1-2 exhibited defects in seedling growth and root elongation. Together, our data suggest that excess sterol accumulation disrupts cellular activities of vegetative organs including leaves and roots, resulting in multiple damages to plants.

Tight DNA-protein complexes isolated from barley seedlings are rich in potential guanine quadruplex sequences

Background

The concept of chromatin domains attached to the nuclear matrix is being revisited, with nucleus described as a set of topologically associating domains. The significance of the tightly bound to DNA proteins (TBP), a protein group that remains attached to DNA after its deproteinization should be also revisited, as the existence of these interactions is in good agreement with the concept of the topologically associating domain. The work aimed to characterize the DNA component of TBP isolated from barley seedlings.

Methods

The tight DNA-protein complexes from the first leaves, coleoptiles, and roots of barley seedlings were isolated by purification with chromatography on nitrocellulose or exhaustive digestion of DNA with DNase I. Cloning and transformation were performed using pMOSBBlue Blunt Ended Cloning Kit. Inserts were amplified by PCR, and sequencing was performed on the MegaBace 1000 Sequencing System. The BLAST search was performed using sequence databases at NCBI, CR-EST, and TREP and Ensembl Plants databases. Comparison to MAR/SAR sequences was performed using http://smartdb.bioinf.med.uni-goettingen.de/cgi-bin/SMARtDB/smar.cgi database. The prediction of G quadruplexes (GQ) was performed with the aid of R-studio library pqsfinder. CD spectra were recorded on a Chirascan CS/3D spectrometer.

Results

Although the barley genome is AT-rich (43% of GC pairs), most DNA fragments associated with TBP were GC-rich (up to 70% in some fractions). Both fractionation procedures yielded a high proportion of CT-motif sequences presented predominantly by the 16-bp CC(TCTCCC)₂ TC fragment present in clones derived from the TBP-bound DNA and absent in free DNA. BLAST analysis revealed alignment with different barley repeats. Some clones, however, aligned with both nuclear and chloroplast structural genes. Alignments with MAR/SAR motifs were very few. The analysis produced by the pqsfinder program revealed numerous potential quadruplex-forming sites in the TBP-bound sequences. A set of oligonucleotides containing sites of possible GQs were designed and ordered. Three of them represented the minus strand of the CT-repeat. Two were derived from sequences of two clones of nitrocellulose retained fraction from leaves and contained GC-rich motifs different from the CT motif. Circular dichroism spectroscopy revealed profound changes in spectra when oligonucleotides were incubated with 100 mM KCl. There was either an increase of positive band in the area of 260 nm or the formation of a positive band at 290 nm. In the former case, changes are typical for parallel G-quadruplexes and, in the latter, 3 + 1 structures.

Discussion

The G-quadruplexes anchor proteins are probably involved in the maintenance of the topologically associated domain structure.

DNA G-quadruplexes: functional significance in plant and farm animal science

G-quadruplexes (G4s) are non-canonical structures that can be formed in DNA and RNA sequences which carry four short runs of guanines. They are distributed in the whole genome but are enriched in gene promoter regions, gene UTRs and chromosome telomeres. The whole array of their functional roles is not fully explored yet but there is solid evidence supporting their implication in a number of processes like regulation of transcription, replication and telomere organization, among others. During the last decade, there is an increased research interest for G4s that has resulted in a better understanding of their role in several physiological and pathological conditions. On the other hand, these structures are poorly studied in plant species and animals of agricultural interest. Here, we summarize the current methods that are used for studying G4s, we review the studies concerning plants and farm animals and we discuss the advantages of a more thorough inclusion of G4s research in the agricultural sciences.

Sequence Dynamics of Pre-mRNA G-Quadruplexes in Plants

Intramolecular G-quadruplexes (G4s) are secondary structures that may form within G-rich stretches of nucleic acids. Although their presence has been associated with genomic instability and mutagenicity, recent reports suggest their involvement in regulation of diverse cellular events, including transcription and translation. The majority of data regarding G4s stems from mammalian and yeast studies, leaving the plant G4s almost unexplored. Using the publicly available Arabidopsis thaliana and Oryza sativa WGS data, we examined the single nucleotide variability of sequences predicted to form G4s (pG4s) structures. We focused our analysis on protein coding transcripts and compared the results to well-characterized Homo sapiens data. We demonstrate that the overall high variability of pG4s is not uniform and differs between gene structural elements. Specifically, plant AUG-containing pG4s, located within 5'UTR/CDS junctions, are abundant and appear not to be affected by a higher frequency of sequence change, indicating their functional relevance. Furthermore, we show that substitutions lowering the probability of G4s' formation are preferred over neutral or stabilizing modifications.

Case studies on potential G-quadruplex-forming sequences from the bacterial orders Deinococcales and Thermales derived from a survey of published genomes

Genomes provide a platform for storage of chemical information that must be stable under the context in which an organism thrives. The 2'-deoxyguanosine (G) nucleotide has the potential to provide additional chemical information beyond its Watson-Crick base-pairing capacity. Sequences with four or more runs of three G nucleotides each are potential G-quadruplex forming sequences (PQSs) that can adopt G-quadruplex folds. Herein, we analyzed sequenced genomes from the NCBI database to determine the PQS densities of the genome sequences. First, we found organisms with large genomes, including humans, alligators, and maize, have similar densities of PQSs (~300 PQSs/Mbp), and the genomes are significantly enriched in PQSs with more than four G tracks. Analysis of microorganism genomes found a greater diversity of PQS densities. In general, PQS densities positively tracked with the GC% of the genome. Exceptions to this observation were the genomes from thermophiles that had many more PQSs than expected by random chance. Analysis of the location of these PQSs in annotated genomes from the order Thermales showed these G-rich sequences to be randomly distributed; in contrast, in the order Deinococcales the PQSs were enriched and biased around transcription start sites of genes. Four representative PQSs, two each from the Thermales and Deinococcales, were studied by biophysical methods to establish the ability of them to fold to G-quadruplexes. The experiments found the two PQSs in the Thermales did not adopt G-quadruplex folds, while the two most common in the Deinococcales adopted stable parallel-stranded G-quadruplexes. The findings lead to a hypothesis that thermophilic organisms are enriched with PQSs as an unavoidable consequence to stabilize thermally their genomes to live at high temperature; in contrast, the genomes from stress-resistant bacteria found in the Deinococcales may utilize PQSs for gene regulatory purposes.

Review: Plant G-quadruplex (G4) motifs in DNA and RNA; abundant, intriguing sequences of unknown function

DNA sequences capable of forming G-quadruplex (G4) structures can be predicted and mapped in plant genomes using computerized pattern search programs. Non-telomeric G4 motifs have recently been found to number in the thousands across many plant species and enriched around gene promoters, prompting speculation that they may represent a newly uncovered and ubiquitous family of cis-acting elements. Comparative analysis shows that monocots exhibit five to ten times higher G4 motif density than eudicots, but the significance of this difference has not been determined. The vast scale and complexity of G4 functions, actual or theoretical, are reviewed in relation to the multiple modes of action and myriad genetic functions for which G4s have been implicated in DNA and RNA. Future experimental strategies and opportunities include identifying plant G4-interactomes, resolving the structures of G4s with and without their binding partners, and defining molecular mechanisms through reporter gene, genetic, or genome editing approaches. Given the global importance of plants for food, clothing, medicine, and energy, together with the potential role of G4 motifs as a widely conserved set of DNA sequences that could coordinate gene regulation, future plant G4 research holds great potential for use in plant improvement strategies.

Genome-wide discovery of G-quadruplex forming sequences and their functional relevance in plants

DNA, in addition to the canonical B-form, can acquire a variety of alternate structures, such as G-quadruplexes. These structures have been implicated in several cellular processes in animals. In this study, we identified different types of G-quadruplex forming sequences (GQSes) in 15 sequenced plants and analyzed their distribution in various genomic features, including gene body, coding, intergenic and promoter regions. G2-type GQSes were most abundant in all the plant species analyzed. A strong association of G3-type GQSes with intergenic, promoter and intronic regions was found. However, G2-type GQSes were enriched in genic, CDS, exonic and untranslated regions. Further, we identified GQSes present in the conserved genes among monocots and dicots. The genes involved in development, cell growth and size, transmembrane transporter, and regulation of gene expression were found to be significantly enriched. In the promoter region, we detected strong co-occurrence of Telobox, ERF, MYB, RAV1B and E2F motifs with GQSes. Further, we validated the structure formation of several plant GQSes, demonstrated their effect on stalling in-vitro replication and revealed their interaction with plant nuclear proteins. Our data provide insights into the prevalence of GQSes in plants, establish their association with different genomic features and functional relevance.

Construction of possible integrated predictive index based on EGFR and ANXA3 polymorphisms for chemotherapy response in fluoropyrimidine-treated Japanese gastric cancer patients using a bioinformatic method

Background

Variability in drug response between individual patients is a serious concern in medicine. To identify single-nucleotide polymorphisms (SNPs) related to drug response variability, many genome-wide association studies have been conducted.

Methods

We previously applied a knowledge-based bioinformatic approach to a pharmacogenomics study in which 119 fluoropyrimidine-treated gastric cancer patients were genotyped at 109,365 SNPs using the Illumina Human-1 BeadChip. We identified the SNP rs2293347 in the human epidermal growth factor receptor (EGFR) gene as a novel genetic factor related to chemotherapeutic response. In the present study, we reanalyzed these hypothesis-free genomic data using extended knowledge.

Results

We identified rs2867461 in annexin A3 (ANXA3) gene as another candidate. Using logistic regression, we confirmed that the performance of the rs2867461 + rs2293347 model was superior to those of the single factor models. Furthermore, we propose a novel integrated predictive index (iEA) based on these two polymorphisms in EGFR and ANXA3. The p value for iEA was 1.47 × 10⁻⁸ by Fisher’s exact test. Recent studies showed that the mutations in EGFR is associated with high expression of dihydropyrimidine dehydrogenase, which is an inactivating and rate-limiting enzyme for fluoropyrimidine, and suggested that the combination of chemotherapy with fluoropyrimidine and EGFR-targeting agents is effective against EGFR-overexpressing gastric tumors, while ANXA3 overexpression confers resistance to tyrosine kinase inhibitors targeting the EGFR pathway.

Conclusions

These results suggest that the iEA index or a combination of polymorphisms in EGFR and ANXA3 may serve as predictive factors of drug response, and therefore could be useful for optimal selection of chemotherapy regimens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1721-z) contains supplementary material, which is available to authorized users.

Strand-Biased Formation of G-Quadruplexes in DNA Duplexes Transcribed with T7 RNA Polymerase

G-quadruplex-forming sequences are enriched near transcription start sites (TSSs) in animal genes. They readily form G-quadruplexes in transcription, which in turn regulate transcription. Therefore, the control of G-quadruplex formation is important for their functionality. It is now shown that G-quadruplexes form efficiently on the non-template, but hardly on the template DNA strand in the downstream vicinity of TSSs in DNA duplexes when they are transcribed by the T7 RNA polymerase (RNAP). Structural analysis reveals that the T7 RNAP causes distortion in a DNA duplex both inside and in front of the enzyme. This structural distortion leads to strand-biased G-quadruplex formation when a G-quadruplex-forming sequence is partially fed into the T7 RNAP to a position about seven nucleotides away from the front of RNA synthesis. Based on these facts, we propose a model for the strand-biased formation of G-quadruplexes in transcribed DNA duplexes.

Flower bud transcriptome analysis of Sapium sebiferum (Linn.) Roxb. and primary investigation of drought induced flowering: pathway construction and G-quadruplex prediction based on transcriptome

Sapium sebiferum (Linn.) Roxb. (Chinese Tallow Tree) is a perennial woody tree and its seeds are rich in oil which hold great potential for biodiesel production. Despite a traditional woody oil plant, our understanding on S. sebiferum genetics and molecular biology remains scant. In this study, the first comprehensive transcriptome of S. sebiferum flower has been generated by sequencing and de novo assembly. A total of 149,342 unigenes were generated from raw reads, of which 24,289 unigenes were successfully matched to public database. A total of 61 MADS box genes and putative pathways involved in S. sebiferum flower development have been identified. Abiotic stress response network was also constructed in this work, where 2,686 unigenes are involved in the pathway. As for lipid biosynthesis, 161 unigenes have been identified in fatty acid (FA) and triacylglycerol (TAG) biosynthesis. Besides, the G-Quadruplexes in RNA of S. sebiferum also have been predicted. An interesting finding is that the stress-induced flowering was observed in S. sebiferum for the first time. According to the results of semi-quantitative PCR, expression tendencies of flowering-related genes, GA1, AP2 and CRY2, accorded with stress-related genes, such as GRX50435 and PRXⅡ39562. This transcriptome provides functional genomic information for further research of S. sebiferum, especially for the genetic engineering to shorten the juvenile period and improve yield by regulating flower development. It also offers a useful database for the research of other Euphorbiaceae family plants.

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.

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.

Bioinformatic analysis reveals an evolutional selection for DNA:RNA hybrid G-quadruplex structures as putative transcription regulatory elements in warm-blooded animals

Recently, we reported the co-transcriptional formation of DNA:RNA hybrid G-quadruplex (HQ) structure by the non-template DNA strand and nascent RNA transcript, which in turn modulates transcription under both in vitro and in vivo conditions. Here we present bioinformatic analysis on putative HQ-forming sequences (PHQS) in the genomes of eukaryotic organisms. Starting from amphibian, PHQS motifs are concentrated in the immediate 1000-nt region downstream of transcription start sites, implying their potential role in transcription regulation. Moreover, their occurrence shows a strong bias toward the non-template versus the template strand. PHQS has become constitutional in genes in warm-blooded animals, and the magnitude of the strand bias correlates with the ability of PHQS to form HQ, suggesting a selection based on HQ formation. This strand bias is reversed in lower species, implying that the selection of PHQS/HQ depended on the living temperature of the organisms. In comparison with the putative intramolecular G-quadruplex-forming sequences (PQS), PHQS motifs are far more prevalent and abundant in the transcribed regions, making them the dominant candidates in the formation of G-quadruplexes in transcription. Collectively, these results suggest that the HQ structures are evolutionally selected to function in transcription and other transcription-mediated processes that involve guanine-rich non-template strand.

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.

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

Leaves develop as flat lateral organs from the indeterminate shoot apical meristem. The establishment of polarity along three-dimensional axes, proximal-distal, medial-lateral, and adaxial-abaxial axes, is crucial for the growth of normal leaves. The mutations of ASYMMETRIC LEAVES1 (AS1) and AS2 of Arabidopsis thaliana cause defects in repression of the indeterminate state and the establishment of axis formation in leaves. Although many mutations have been identified that enhance the adaxial-abaxial polarity defects of as1 and as2 mutants, the roles of the causative genes in leaf development are still unknown. In this study, we found that wild-type plants treated with berberine produced pointed leaves, which are often observed in the single mutants that enhance phenotypes of as1 and as2 mutants. The berberine-treated as1 and as2 mutants formed abaxialized filamentous leaves. Berberine, an isoquinoline alkaloid compound naturally produced in various plant sources, has a growth inhibitory effect on plants that do not produce berberine. We further showed that transcript levels of meristem-specific class 1 KNOX homeobox genes and abaxial determinant genes were increased in berberine-treated as1 and as2. Berberine treated plants carrying double mutations of AS2 and the large subunit ribosomal protein gene RPL5B showed more severe defects in polarity than did the as2 single mutant plants. We suggest that berberine inhibits (a) factor(s) that might be required for leaf adaxial cell differentiation through a pathway independent of AS1 and AS2. Multiple pathways might play important roles in the formation of flat symmetric leaves.