Boosting AthaMap Database Content with Data from Protein Binding Microarrays

Machine Learning Inference of Gene Regulatory Networks in Developing Mimulus Seeds

The angiosperm seed represents a critical evolutionary breakthrough that has been shown to propel the reproductive success and radiation of flowering plants. Seeds promote the rapid diversification of angiosperms by establishing postzygotic reproductive barriers, such as hybrid seed inviability. While prezygotic barriers to reproduction tend to be transient, postzygotic barriers are often permanent and therefore can play a pivotal role in facilitating speciation. This property of the angiosperm seed is exemplified in the Mimulus genus. In order to further the understanding of the gene regulatory mechanisms important in the Mimulus seed, we performed gene regulatory network (GRN) inference analysis by using time-series RNA-seq data from developing hybrid seeds from a viable cross between Mimulus guttatus and Mimulus pardalis. GRN inference has the capacity to identify active regulatory mechanisms in a sample and highlight genes of potential biological importance. In our case, GRN inference also provided the opportunity to uncover active regulatory relationships and generate a reference set of putative gene regulations. We deployed two GRN inference algorithms-RTP-STAR and KBoost-on three different subsets of our transcriptomic dataset. While the two algorithms yielded GRNs with different regulations and topologies when working with the same data subset, there was still significant overlap in the specific gene regulations they inferred, and they both identified potential novel regulatory mechanisms that warrant further investigation.

SnRK2.10 kinase differentially modulates expression of hub WRKY transcription factors genes under salinity and oxidative stress in Arabidopsis thaliana

In nature, all living organisms must continuously sense their surroundings and react to the occurring changes. In the cell, the information about these changes is transmitted to all cellular compartments, including the nucleus, by multiple phosphorylation cascades. Sucrose Non-Fermenting 1 Related Protein Kinases (SnRK2s) are plant-specific enzymes widely distributed across the plant kingdom and key players controlling abscisic acid (ABA)-dependent and ABA-independent signaling pathways in the plant response to osmotic stress and salinity. The main deleterious effects of salinity comprise water deficiency stress, disturbances in ion balance, and the accompanying appearance of oxidative stress. The reactive oxygen species (ROS) generated at the early stages of salt stress are involved in triggering intracellular signaling required for the fast stress response and modulation of gene expression. Here we established in Arabidopsis thaliana that salt stress or induction of ROS accumulation by treatment of plants with H2O2 or methyl viologen (MV) induces the expression of several genes encoding transcription factors (TFs) from the WRKY DNA-Binding Protein (WRKY) family. Their induction by salinity was dependent on SnRK2.10, an ABA non-activated kinase, as it was strongly reduced in snrk2.10 mutants. The effect of ROS was clearly dependent on their source. Following the H2O2 treatment, SnRK2.10 was activated in wild-type (wt) plants and the induction of the WRKY TFs expression was only moderate and was enhanced in snrk2.10 lines. In contrast, MV did not activate SnRK2.10 and the WRKY induction was very strong and was similar in wt and snrk2.10 plants. A bioinformatic analysis indicated that the WRKY33, WRKY40, WRKY46, and WRKY75 transcription factors have a similar target range comprising numerous stress-responsive protein kinases. Our results indicate that the stress-related functioning of SnRK2.10 is fine-tuned by the source and intracellular distribution of ROS and the co-occurrence of other stress factors.

Design and Analysis of Native Photorespiration Gene Motifs of Promoter Untranslated Region Combinations Under Short Term Abiotic Stress Conditions

Quantitative traits are rarely controlled by a single gene, thereby making multi-gene transformation an indispensable component of modern synthetic biology approaches. However, the shortage of unique gene regulatory elements (GREs) for the robust simultaneous expression of multiple nuclear transgenes is a major bottleneck that impedes the engineering of complex pathways in plants. In this study, we compared the transcriptional efficacies of a comprehensive list of well-documented promoter and untranslated region (UTR) sequences side by side. The strength of GREs was examined by a dual-luciferase assay in conjunction with transient expression in tobacco. In addition, we created suites of new GREs with higher transcriptional efficacies by combining the best performing promoter-UTR sequences. We also tested the impact of elevated temperature and high irradiance on the effectiveness of these GREs. While constitutive promoters ensure robust expression of transgenes, they lack spatiotemporal regulations exhibited by native promoters. Here, we present a proof-of-principle study on the characterization of synthetic promoters based on cis-regulatory elements of three key photorespiratory genes. This conserved biochemical process normally increases under elevated temperature, low CO2, and high irradiance stress conditions and results in ∼25% loss in fixed CO2. To select stress-responsive cis-regulatory elements involved in photorespiration, we analyzed promoters of two chloroplast transporters (AtPLGG1 and AtBASS6) and a key plastidial enzyme, AtPGLP using PlantPAN3.0 and AthaMap. Our results suggest that these motifs play a critical role for PLGG1, BASS6, and PGLP in mediating response to elevated temperature and high-intensity light stress. These findings will not only enable the advancement of metabolic and genetic engineering of photorespiration but will also be instrumental in related synthetic biology approaches.

Statistical estimates of multiple transcription factors binding in the model plant genomes based on ChIP-seq data

The development of high-throughput genomic sequencing coupled with chromatin immunoprecipitation technologies allows studying the binding sites of the protein transcription factors (TF) in the genome scale. The growth of data volume on the experimentally determined binding sites raises qualitatively new problems for the analysis of gene expression regulation, prediction of transcription factors target genes, and regulatory gene networks reconstruction. Genome regulation remains an insufficiently studied though plants have complex molecular regulatory mechanisms of gene expression and response to environmental stresses. It is important to develop new software tools for the analysis of the TF binding sites location and their clustering in the plant genomes, visualization, and the following statistical estimates. This study presents application of the analysis of multiple TF binding profiles in three evolutionarily distant model plant organisms. The construction and analysis of non-random ChIP-seq binding clusters of the different TFs in mammalian embryonic stem cells were discussed earlier using similar bioinformatics approaches. Such clusters of TF binding sites may indicate the gene regulatory regions, enhancers and gene transcription regulatory hubs. It can be used for analysis of the gene promoters as well as a background for transcription networks reconstruction. We discuss the statistical estimates of the TF binding sites clusters in the model plant genomes. The distributions of the number of different TFs per binding cluster follow same power law distribution for all the genomes studied. The binding clusters in Arabidopsis thaliana genome were discussed here in detail.

NIN-like protein7 and PROTEOLYSIS6 functional interaction enhances tolerance to sucrose, ABA, and submergence

Nitrate (NO3) assimilation and signaling regulate plant growth through the relevant function of the transcription factor NIN-like Protein7 (NLP7). NO3 is also the main source for plants to produce nitric oxide (NO), which regulates growth and stress responses. NO-mediated regulation requires efficient sensing via the PROTEOLYSIS6 (PRT6)-mediated proteasome-triggered degradation of group VII of ethylene response transcription factors through the Cys/Arg N-degron pathway. The convergence of NO3 signaling and N-degron proteolysis on NO-mediated regulation remains largely unknown. Here, we investigated the functional interaction between NLP7 and PRT6 using Arabidopsis (Arabidopsis thaliana) double prt6 nlp7 mutant plants as well as complementation lines overexpressing NLP7 in different mutant genetic backgrounds. prt6 nlp7 mutant plants displayed several potentiated prt6 characteristic phenotypes, including slower vegetative growth, increased NO content, and diminished tolerance to abiotic stresses such as high-sucrose concentration, abscisic acid, and hypoxia-reoxygenation. Although NLP7 has an N-terminus that could be targeted by the N-degron proteolytic pathway, it was not a PRT6 substrate. The potential PRT6- and NO-regulated nucleocytoplasmic translocation of NLP7, which is likely modulated by posttranslational modifications, is proposed to act as a regulatory loop to control NO homeostasis and action.

The Plastid-Localized AtFtsHi3 Pseudo-Protease of Arabidopsis thaliana Has an Impact on Plant Growth and Drought Tolerance

While drought severely affects plant growth and crop production, the molecular mechanisms of the drought response of plants remain unclear. In this study, we demonstrated for the first time the effect of the pseudo-protease AtFtsHi3 of Arabidopsis thaliana on overall plant growth and in drought tolerance. An AtFTSHi3 knock-down mutant [ftshi3-1(kd)] displayed a pale-green phenotype with lower photosynthetic efficiency and Darwinian fitness compared to wild type (Wt). An observed delay in seed germination of ftshi3-1(kd) was attributed to overaccumulation of abscisic acid (ABA); ftshi3-1(kd) seedlings showed partial sensitivity to exogenous ABA. Being exposed to similar severity of soil drying, ftshi3-1(kd) was drought-tolerant up to 20 days after the last irrigation, while wild type plants wilted after 12 days. Leaves of ftshi3-1(kd) contained reduced stomata size, density, and a smaller stomatic aperture. During drought stress, ftshi3-1(kd) showed lowered stomatal conductance, increased intrinsic water-use efficiency (WUEi), and slower stress acclimation. Expression levels of ABA-responsive genes were higher in leaves of ftshi3-1(kd) than Wt; DREB1A, but not DREB2A, was significantly upregulated during drought. However, although ftshi3-1(kd) displayed a drought-tolerant phenotype in aboveground tissue, the root-associated bacterial community responded to drought.

Negative modulation of SA signaling components by the capsid protein of tobacco mosaic virus is required for viral long-distance movement

An important aspect of plant-virus interaction is the way viruses dynamically move over long distances and how plant immunity modulates viral systemic movement. Salicylic acid (SA), a well-characterized hormone responsible for immune responses against virus, is activated through different transcription factors including TGA and WRKY. In tobamoviruses, evidence suggests that capsid protein (CP) is required for long-distance movement, although its precise role has not been fully characterized yet. Previously, we showed that the CP of Tobacco Mosaic Virus (TMV)-Cg negatively modulates the SA-mediated defense. In this study, we analyzed the impact of SA-defense mechanism on the long-distance transport of a truncated version of TMV (TMV ∆CP virus) that cannot move to systemic tissues. The study showed that the negative modulation of NPR1 and TGA10 factors allows the long-distance transport of TMV ∆CP virus. Moreover, we observed that the stabilization of DELLA proteins promotes TMV ∆CP systemic movement. We also characterized a group of genes, part of a network modulated by CP, involved in TMV ∆CP long-distance transport. Altogether, our results indicate that CP-mediated downregulation of SA signaling pathway is required for the virus systemic movement, and this role of CP may be linked to its ability to stabilize DELLA proteins.

Signaling Toward Reactive Oxygen Species-Scavenging Enzymes in Plants

Reactive oxygen species (ROS) are signaling molecules essential for plant responses to abiotic and biotic stimuli as well as for multiple developmental processes. They are produced as byproducts of aerobic metabolism and are affected by adverse environmental conditions. The ROS content is controlled on the side of their production but also by scavenging machinery. Antioxidant enzymes represent a major ROS-scavenging force and are crucial for stress tolerance in plants. Enzymatic antioxidant defense occurs as a series of redox reactions for ROS elimination. Therefore, the deregulation of the antioxidant machinery may lead to the overaccumulation of ROS in plants, with negative consequences both in terms of plant development and resistance to environmental challenges. The transcriptional activation of antioxidant enzymes accompanies the long-term exposure of plants to unfavorable environmental conditions. Fast ROS production requires the immediate mobilization of the antioxidant defense system, which may occur via retrograde signaling, redox-based modifications, and the phosphorylation of ROS detoxifying enzymes. This review aimed to summarize the current knowledge on signaling processes regulating the enzymatic antioxidant capacity of plants.

Unusual DNA-binding properties of the Arabidopsis thaliana WRKY50 transcription factor at target gene promoters

WRKY50 from A. thaliana requires WT-boxes at target gene promoters for activation and binding. Based on the genome-wide prediction of WRKY50 target genes and the similarity of a WRKY50 binding site to WT-boxes in microbe-associated molecular pattern (MAMP)-responsive cis-regulatory modules (CRM), four WT-box containing CRMs from the promoter region of three WRKY50 target genes were investigated for their interaction with WRKY50. These target genes are DJ1E, WRKY30 and ATBBE4. Two of the four CRMs, one from DJ1E and one from WRKY30, were able to activate reporter gene expression in the presence of WRKY50. Activation requires the WT-boxes GGACTTTT, GGACTTTG from DJ1E and GGACTTTC from WRKY30. WRKY50 does not activate a second CRM from WRKY30 and the CRM from ATBBE4, both containing the WT-box TGACTTTT. In vitro gel-shift assays demonstrate WT-box-specific binding of the WRKY50 DNA-binding domain to all four CRMs. This work shows a high flexibility of WRKY50 binding site recognition beyond the classic W-box TTGACC/T.

A strong NF-κB p65 responsive cis-regulatory sequence from Arabidopsis thaliana interacts with WRKY40

Transcription factors from mammals and plants, which play a role in innate immunity, interact with the same microbe-associated molecular pattern (MAMP)-responsive sequences from Arabidopsis thaliana. The interaction of mouse NF-κB p65 with MAMP-responsive sequences containing the core motif GACTTT of the WT-box was investigated. This revealed one sequence, derived from the promoter of the A. thaliana gene At1g76960, a gene with unknown function, to activate NF-κB p65 dependent reporter gene expression in plant cells very strongly. A bioinformatic analysis predicts three putative NF-κB p65 binding sites in this sequence and all three sites are required for reporter gene activation and binding. The sequence is one of the weakest MAMP-responsive sequences previously isolated, but the introduction of a GCC-box increases its MAMP responsivity in combination with upstream WT-box sequences. Although a bioinformatic analysis of the unmutated cis-sequence only predicts NF-κB p65 binding, A. thaliana WRKY40 was selected in a yeast one-hybrid screen. WRKY40, which is a transcriptional repressor, requires the sequence TTTTCTA for direct binding. This sequence is similar to the WK-box TTTTCCAC, previously shown to interact with tobacco NtWRKY12. In summary, this work supports the similarity in binding site recognition between NF-κB and WRKY factors.

Arabidopsis mlo3 mutant plants exhibit spontaneous callose deposition and signs of early leaf senescence

Arabidopsis thaliana mlo3 mutant plants are not affected in pathogen infection phenotypes but-reminiscent of mlo2 mutant plants-exhibit spontaneous callose deposition and signs of early leaf senescence. The family of Mildew resistance Locus O (MLO) proteins is best known for its profound effect on the outcome of powdery mildew infections: when the appropriate MLO protein is absent, the plant is fully resistant to otherwise virulent powdery mildew fungi. However, most members of the MLO protein family remain functionally unexplored. Here, we investigate Arabidopsis thaliana MLO3, the closest relative of AtMLO2, AtMLO6 and AtMLO12, which are the Arabidopsis MLO genes implicated in the powdery mildew interaction. The co-expression network of AtMLO3 suggests association of the gene with plant defense-related processes such as salicylic acid homeostasis. Our extensive analysis shows that mlo3 mutants are unaffected regarding their infection phenotype upon challenge with the powdery mildew fungi Golovinomyces orontii and Erysiphe pisi, the oomycete Hyaloperonospora arabidopsidis, and the bacterial pathogen Pseudomonas syringae (the latter both in terms of basal and systemic acquired resistance), indicating that the protein does not play a major role in the response to any of these pathogens. However, mlo3 genotypes display spontaneous callose deposition as well as signs of early senescence in 6- or 7-week-old rosette leaves in the absence of any pathogen challenge, a phenotype that is reminiscent of mlo2 mutant plants. We hypothesize that de-regulated callose deposition in mlo3 genotypes might be the result of a subtle transient aberration of salicylic acid-jasmonic acid homeostasis during development.

Comparative transcriptome meta-analysis of Arabidopsis thaliana under drought and cold stress

Multiple environmental stresses adversely affect plant growth and development. Plants under multiple stress condition trigger cascade of signals and show response unique to specific stress as well as shared responses, common to individual stresses. Here, we aim to identify common and unique genetic components during stress response mechanisms liable for cross-talk between stresses. Although drought and cold stress have been widely studied, insignificant information is available about how their combination affects plants. To that end, we performed meta-analysis and co-expression network comparison of drought and cold stress response in Arabidopsis thaliana by analyzing 390 microarray samples belonging to 29 microarray studies. We observed 6120 and 7079 DEGs (differentially expressed genes) under drought and cold stress respectively, using Rank Product methodology. Statistically, 28% (2890) DEGs were found to be common in both the stresses (i.e.; drought and cold stress) with most of them having similar expression pattern. Further, gene ontology-based enrichment analysis have identified shared biological processes and molecular mechanisms such as—‘photosynthesis’, ‘respiratory burst’, ‘response to hormone’, ‘signal transduction’, ‘metabolic process’, ‘response to water deprivation’, which were affected under cold and drought stress. Forty three transcription factor families were found to be expressed under both the stress conditions. Primarily, WRKY, NAC, MYB, AP2/ERF and bZIP transcription factor family genes were highly enriched in all genes sets and were found to regulate 56% of common genes expressed in drought and cold stress. Gene co-expression network analysis by WGCNA (weighted gene co-expression network analysis) revealed 21 and 16 highly inter-correlated gene modules with specific expression profiles under drought and cold stress respectively. Detection and analysis of gene modules shared between two stresses revealed the presence of four consensus gene modules.

Eukaryotic and prokaryotic promoter databases as valuable tools in exploring the regulation of gene transcription: a comprehensive overview

The complete exploration of the regulation of gene expression remains one of the top-priority goals for researchers. As the regulation is mainly controlled at the level of transcription by promoters, study on promoters and findings are of great importance. This review summarizes forty selected databases that centralize experimental and theoretical knowledge regarding the organization of promoters, interacting transcription factors (TFs) and microRNAs (miRNAs) in many eukaryotic and prokaryotic species. The presented databases offer researchers valuable support in elucidating the regulation of gene transcription.

Genome-wide chromatin mapping with size resolution reveals a dynamic sub-nucleosomal landscape in Arabidopsis

All eukaryotic genomes are packaged as chromatin, with DNA interlaced with both regularly patterned nucleosomes and sub-nucleosomal-sized protein structures such as mobile and labile transcription factors (TF) and initiation complexes, together forming a dynamic chromatin landscape. Whilst details of nucleosome position in Arabidopsis have been previously analysed, there is less understanding of their relationship to more dynamic sub-nucleosomal particles (subNSPs) defined as protected regions shorter than the ~150bp typical of nucleosomes. The genome-wide profile of these subNSPs has not been previously analysed in plants and this study investigates the relationship of dynamic bound particles with transcriptional control. Here we combine differential micrococcal nuclease (MNase) digestion and a modified paired-end sequencing protocol to reveal the chromatin structure landscape of Arabidopsis cells across a wide particle size range. Linking this data to RNAseq expression analysis provides detailed insight into the relationship of identified DNA-bound particles with transcriptional activity. The use of differential digestion reveals sensitive positions, including a labile -1 nucleosome positioned upstream of the transcription start site (TSS) of active genes. We investigated the response of the chromatin landscape to changes in environmental conditions using light and dark growth, given the large transcriptional changes resulting from this simple alteration. The resulting shifts in the suites of expressed and repressed genes show little correspondence to changes in nucleosome positioning, but led to significant alterations in the profile of subNSPs upstream of TSS both globally and locally. We examined previously mapped positions for the TFs PIF3, PIF4 and CCA1, which regulate light responses, and found that changes in subNSPs co-localized with these binding sites. This small particle structure is detected only under low levels of MNase digestion and is lost on more complete digestion of chromatin to nucleosomes. We conclude that wide-spectrum analysis of the Arabidopsis genome by differential MNase digestion allows detection of sensitive features hereto obscured, and the comparisons between genome-wide subNSP profiles reveals dynamic changes in their distribution, particularly at distinct genomic locations (i.e. 5'UTRs). The method here employed allows insight into the complex influence of genetic and extrinsic factors in modifying the sub-nucleosomal landscape in association with transcriptional changes.

From experiment-driven database analyses to database-driven experiments in Arabidopsis thaliana transcription factor research

Experiment-driven database analysis is employed in forward genetics to predict the function of genes assocíated with a mutant phenotype. These analyses subsequently lead to database-driven experiments involving reverse genetics to verify functional predictions based on bioinformatic analyses. Genomic transcription factors (TFs) are key regulators of gene expression by binding to short regulatory sequences and by interacting with other TFs. Currently more than 2400 TFs are predicted for A. thaliana. As DNA-binding proteins they are particularly amenable to database-driven experiments, especially when their binding site specificities are known. Databases are available for predicting binding sites for specific TFs in regulatory sequences. Since most of these bioinformatically identified binding sites may not be functional, additional experiments for identifying the actual in vivo binding sites for TFs are required. Recently, large scale approaches were employed to determine binding sites for many A. thaliana TFs. With these approaches binding sites for 984 unique TFs were determined experimentally. An area deserving further research is proposed for interacting TFs. Most of the A. thaliana genes are under combinatorial control, and in vivo interacting TFs, similar to mammalian TFs, may bind to combinatorial elements in which the binding sites vary from those detected with the single TFs.

High-quality genome assembly of Capsella bursa-pastoris reveals asymmetry of regulatory elements at early stages of polyploid genome evolution

Polyploidization and subsequent sub- and neofunctionalization of duplicated genes represent a major mechanism of plant genome evolution. Capsella bursa-pastoris, a widespread ruderal plant, is a recent allotetraploid and, thus, is an ideal model organism for studying early changes following polyploidization. We constructed a high-quality assembly of C. bursa-pastoris genome and a transcriptome atlas covering a broad sample of organs and developmental stages (available online at http://travadb.org/browse/Species=Cbp). We demonstrate that expression of homeologs is mostly symmetric between subgenomes, and identify a set of homeolog pairs with discordant expression. Comparison of promoters within such pairs revealed emerging asymmetry of regulatory elements. Among them there are multiple binding sites for transcription factors controlling the regulation of photosynthesis and plant development by light (PIF3, HY5) and cold stress response (CBF). These results suggest that polyploidization in C. bursa-pastoris enhanced its plasticity of response to light and temperature, and allowed substantial expansion of its distribution range.

Computational inference of gene regulatory networks: Approaches, limitations and opportunities

Gene regulatory networks lie at the core of cell function control. In E. coli and S. cerevisiae, the study of gene regulatory networks has led to the discovery of regulatory mechanisms responsible for the control of cell growth, differentiation and responses to environmental stimuli. In plants, computational rendering of gene regulatory networks is gaining momentum, thanks to the recent availability of high-quality genomes and transcriptomes and development of computational network inference approaches. Here, we review current techniques, challenges and trends in gene regulatory network inference and highlight challenges and opportunities for plant science. We provide plant-specific application examples to guide researchers in selecting methodologies that suit their particular research questions. Given the interdisciplinary nature of gene regulatory network inference, we tried to cater to both biologists and computer scientists to help them engage in a dialogue about concepts and caveats in network inference. Specifically, we discuss problems and opportunities in heterogeneous data integration for eukaryotic organisms and common caveats to be considered during network model evaluation. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Bioinformatic Identification of Conserved Cis-Sequences in Coregulated Genes

Bioinformatics tools can be employed to identify conserved cis-sequences in sets of coregulated plant genes because more and more gene expression and genomic sequence data become available. Knowledge on the specific cis-sequences, their enrichment and arrangement within promoters, facilitates the design of functional synthetic plant promoters that are responsive to specific stresses. The present chapter illustrates an example for the bioinformatic identification of conserved Arabidopsis thaliana cis-sequences enriched in drought stress-responsive genes. This workflow can be applied for the identification of cis-sequences in any sets of coregulated genes. The workflow includes detailed protocols to determine sets of coregulated genes, to extract the corresponding promoter sequences, and how to install and run a software package to identify overrepresented motifs. Further bioinformatic analyses that can be performed with the results are discussed.

More than meets the eye: Emergent properties of transcription factors networks in Arabidopsis

Uncovering and mathematically modeling Transcription Factor Networks (TFNs) are the first steps in engineering plants with traits that are better equipped to respond to changing environments. Although several plant TFNs are well known, the framework for systematically modeling complex characteristics such as switch-like behavior, oscillations, and homeostasis that emerge from them remain elusive. This review highlights literature that provides, in part, experimental and computational techniques for characterizing TFNs. This review also outlines methodologies that have been used to mathematically model the dynamic characteristics of TFNs. We present several examples of TFNs in plants that are involved in developmental and stress response. In several cases, advanced algorithms capture or quantify emergent properties that serve as the basis for robustness and adaptability in plant responses. Increasing the use of mathematical approaches will shed new light on these regulatory properties that control plant growth and development, leading to mathematical models that predict plant behavior. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Purifying selection acts on coding and non-coding sequences of paralogous genes in Arabidopsis thaliana

Background

Whole-genome duplications in the ancestors of many diverse species provided the genetic material for evolutionary novelty. Several models explain the retention of paralogous genes. However, how these models are reflected in the evolution of coding and non-coding sequences of paralogous genes is unknown.

Results

Here, we analyzed the coding and non-coding sequences of paralogous genes in Arabidopsis thaliana and compared these sequences with those of orthologous genes in Arabidopsis lyrata. Paralogs with lower expression than their duplicate had more nonsynonymous substitutions, were more likely to fractionate, and exhibited less similar expression patterns with their orthologs in the other species. Also, lower-expressed genes had greater tissue specificity. Orthologous conserved non-coding sequences in the promoters, introns, and 3′ untranslated regions were less abundant at lower-expressed genes compared to their higher-expressed paralogs. A gene ontology (GO) term enrichment analysis showed that paralogs with similar expression levels were enriched in GO terms related to ribosomes, whereas paralogs with different expression levels were enriched in terms associated with stress responses.

Conclusions

Loss of conserved non-coding sequences in one gene of a paralogous gene pair correlates with reduced expression levels that are more tissue specific. Together with increased mutation rates in the coding sequences, this suggests that similar forces of purifying selection act on coding and non-coding sequences. We propose that coding and non-coding sequences evolve concurrently following gene duplication.

Electronic supplementary material

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

Identification of plant transcription factor target sequences

Regulation of gene expression depends on specific cis-regulatory sequences located in the gene promoter regions. These DNA sequences are recognized by transcription factors (TFs) in a sequence-specific manner, and their identification could help to elucidate the regulatory networks that underlie plant physiological responses to developmental programs or to environmental adaptation. Here we review recent advances in high throughput methodologies for the identification of plant TF binding sites. Several approaches offer a map of the TF binding locations in vivo and of the dynamics of the gene regulatory networks. As an alternative, high throughput in vitro methods provide comprehensive determination of the DNA sequences recognized by TFs. These advances are helping to decipher the regulatory lexicon and to elucidate transcriptional network hierarchies in plants in response to internal or external cues. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.