Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators.

Gene regulation is a dynamic process in which transcription factors (TFs) play an important role in controlling spatiotemporal gene expression. To enhance our global understanding of regulatory interactions in Arabidopsis thaliana, different regulatory input networks capturing complementary information about DNA motifs, open chromatin, TF-binding and expression-based regulatory interactions were combined using a supervised learning approach, resulting in an integrated gene regulatory network (iGRN) covering 1,491 TFs and 31,393 target genes (1.7 million interactions). This iGRN outperforms the different input networks to predict known regulatory interactions and has a similar performance to recover functional interactions compared to state-of-the-art experimental methods. The iGRN correctly inferred known functions for 681 TFs and predicted new gene functions for hundreds of unknown TFs. For regulators predicted to be involved in reactive oxygen species (ROS) stress regulation, we confirmed in total 75% of TFs with a function in ROS and/or physiological stress responses. This includes 13 ROS regulators, previously not connected to any ROS or stress function, that were experimentally validated in our ROS-specific phenotypic assays of loss- or gain-of-function lines. In conclusion, the presented iGRN offers a high-quality starting point to enhance our understanding of gene regulation in plants by integrating different experimental data types.

Understanding genetic control of root system architecture in soybean: Insights into the genetic basis of lateral root number.

Developing crops with better root systems is a promising strategy to ensure productivity in both optimum and stress environments. Root system architectural traits in 397 soybean accessions were characterized and a high-density single nucleotide polymorphisms (SNPs)-based genome-wide association study was performed to identify the underlying genes associated with root structure. SNPs associated with root architectural traits specific to landraces and elite germplasm pools were detected. Four loci were detected in landraces for lateral root number (LRN) and distribution of root thickness in diameter Class I with a major locus on chromosome 16. This major loci was detected in the coding region of unknown protein, and subsequent analyses demonstrated that root traits are affected with mutated haplotypes of the gene. In elite germplasm pool, 3 significant SNPs in alanine-glyoxalate aminotransferase, Leucine-Rich Repeat receptor/No apical meristem, and unknown functional genes were found to govern multiple traits including root surface area and volume. However, no major loci were detected for LRN in elite germplasm. Nucleotide diversity analysis found evidence of selective sweeps around the landraces LRN gene. Soybean accessions with minor and mutated allelic variants of LRN gene were found to perform better in both water-limited and optimal field conditions.

KIRA1 and ORESARA1 terminate flower receptivity by promoting cell death in the stigma of Arabidopsis.

Flowers have a species-specific functional life span that determines the time window in which pollination, fertilization and seed set can occur. The stigma tissue plays a key role in flower receptivity by intercepting pollen and initiating pollen tube growth toward the ovary. In this article, we show that a developmentally controlled cell death programme terminates the functional life span of stigma cells in Arabidopsis. We identified the leaf senescence regulator ORESARA1 (also known as ANAC092) and the previously uncharacterized KIRA1 (also known as ANAC074) as partially redundant transcription factors that modulate stigma longevity by controlling the expression of programmed cell death-associated genes. KIRA1 expression is sufficient to induce cell death and terminate floral receptivity, whereas lack of both KIRA1 and ORESARA1 substantially increases stigma life span. Surprisingly, the extension of stigma longevity is accompanied by only a moderate extension of flower receptivity, suggesting that additional processes participate in the control of the flower's receptive life span.

TF2Network: predicting transcription factor regulators and gene regulatory networks in Arabidopsis using publicly available binding site information.

A gene regulatory network (GRN) is a collection of regulatory interactions between transcription factors (TFs) and their target genes. GRNs control different biological processes and have been instrumental to understand the organization and complexity of gene regulation. Although various experimental methods have been used to map GRNs in Arabidopsis thaliana, their limited throughput combined with the large number of TFs makes that for many genes our knowledge about regulating TFs is incomplete. We introduce TF2Network, a tool that exploits the vast amount of TF binding site information and enables the delineation of GRNs by detecting potential regulators for a set of co-expressed or functionally related genes. Validation using two experimental benchmarks reveals that TF2Network predicts the correct regulator in 75-92% of the test sets. Furthermore, our tool is robust to noise in the input gene sets, has a low false discovery rate, and shows a better performance to recover correct regulators compared to other plant tools. TF2Network is accessible through a web interface where GRNs are interactively visualized and annotated with various types of experimental functional information. TF2Network was used to perform systematic functional and regulatory gene annotations, identifying new TFs involved in circadian rhythm and stress response.

The Transcription Factor MYB29 Is a Regulator of ALTERNATIVE OXIDASE1a.

Plants sense and integrate a variety of signals from the environment through different interacting signal transduction pathways that involve hormones and signaling molecules. Using ALTERNATIVE OXIDASE1a (AOX1a) gene expression as a model system of retrograde or stress signaling between mitochondria and the nucleus, MYB DOMAIN PROTEIN29 (MYB29) was identified as a negative regulator (regulator of alternative oxidase1a 7 [rao7] mutant) in a genetic screen of Arabidopsis (Arabidopsis thaliana). rao7/myb29 mutants have increased levels of AOX1a transcript and protein compared to wild type after induction with antimycin A. A variety of genes previously associated with the mitochondrial stress response also display enhanced transcript abundance, indicating that RAO7/MYB29 negatively regulates mitochondrial stress responses in general. Meta-analysis of hormone-responsive marker genes and identification of downstream transcription factor networks revealed that MYB29 functions in the complex interplay of ethylene, jasmonic acid, salicylic acid, and reactive oxygen species signaling by regulating the expression of various ETHYLENE RESPONSE FACTOR and WRKY transcription factors. Despite an enhanced induction of mitochondrial stress response genes, rao7/myb29 mutants displayed an increased sensitivity to combined moderate light and drought stress. These results uncover interactions between mitochondrial retrograde signaling and the regulation of glucosinolate biosynthesis, both regulated by RAO7/MYB29. This common regulator can explain why perturbation of the mitochondrial function leads to transcriptomic responses overlapping with responses to biotic stress.

Physical Interaction of T Cells with Dendritic Cells Is Not Required for the Immunomodulatory Effects of the Edible Mushroom Agaricus subrufescens.

Mushrooms are well known for their immunomodulating capacities. However, little is known about how mushroom-stimulated dendritic cells (DCs) affect T cells. Therefore, we investigated the effect of mushroom compounds derived from seven edible mushroom species on DCs, their fate in DCs, and the effect of the mushroom-stimulated DCs on T cells. Each mushroom species stimulated DCs in a different manner as was revealed from the DC's cytokine response. Assessing DC maturation revealed that only one mushroom species, Agaricus subrufescens, induced complete DC maturation. The other six mushroom species upregulated MHC-II and CD86 expression, but did not significantly affect the expression of CD40 and CD11c. Nevertheless, mushroom compounds of all investigated mushroom species are endocytosed by DCs. Endocytosis is most likely mediated by C-type lectin receptors (CLRs) because CLR binding is Ca2+ dependent, and EGTA reduces TNF-α secretion with more than 90%. Laminarin partly inhibited TNF-α secretion indicating that the CLR dectin-1, among other CLRs, is involved in binding mushroom compounds. Stimulated DCs were shown to stimulate T cells; however, physical contact of DCs and T cells is not required. Because CLRs seem to play a prominent role in DC stimulation, mushrooms may function as a carbohydrate containing adjuvant to be used in conjunction with anti-fungal vaccines.

A Collection of Conserved Noncoding Sequences to Study Gene Regulation in Flowering Plants.

Transcription factors (TFs) regulate gene expression by binding cis-regulatory elements, of which the identification remains an ongoing challenge owing to the prevalence of large numbers of nonfunctional TF binding sites. Powerful comparative genomics methods, such as phylogenetic footprinting, can be used for the detection of conserved noncoding sequences (CNSs), which are functionally constrained and can greatly help in reducing the number of false-positive elements. In this study, we applied a phylogenetic footprinting approach for the identification of CNSs in 10 dicot plants, yielding 1,032,291 CNSs associated with 243,187 genes. To annotate CNSs with TF binding sites, we made use of binding site information for 642 TFs originating from 35 TF families in Arabidopsis (Arabidopsis thaliana). In three species, the identified CNSs were evaluated using TF chromatin immunoprecipitation sequencing data, resulting in significant overlap for the majority of data sets. To identify ultraconserved CNSs, we included genomes of additional plant families and identified 715 binding sites for 501 genes conserved in dicots, monocots, mosses, and green algae. Additionally, we found that genes that are part of conserved mini-regulons have a higher coherence in their expression profile than other divergent gene pairs. All identified CNSs were integrated in the PLAZA 3.0 Dicots comparative genomics platform (http://bioinformatics.psb.ugent.be/plaza/versions/plaza_v3_dicots/) together with new functionalities facilitating the exploration of conserved cis-regulatory elements and their associated genes. The availability of this data set in a user-friendly platform enables the exploration of functional noncoding DNA to study gene regulation in a variety of plant species, including crops.

Functional Analysis of the Arabidopsis TETRASPANIN Gene Family in Plant Growth and Development.

TETRASPANIN (TET) genes encode conserved integral membrane proteins that are known in animals to function in cellular communication during gamete fusion, immunity reaction, and pathogen recognition. In plants, functional information is limited to one of the 17 members of the Arabidopsis (Arabidopsis thaliana) TET gene family and to expression data in reproductive stages. Here, the promoter activity of all 17 Arabidopsis TET genes was investigated by pAtTET::NUCLEAR LOCALIZATION SIGNAL-GREEN FLUORESCENT PROTEIN/ß-GLUCURONIDASE reporter lines throughout the life cycle, which predicted functional divergence in the paralogous genes per clade. However, partial overlap was observed for many TET genes across the clades, correlating with few phenotypes in single mutants and, therefore, requiring double mutant combinations for functional investigation. Mutational analysis showed a role for TET13 in primary root growth and lateral root development and redundant roles for TET5 and TET6 in leaf and root growth through negative regulation of cell proliferation. Strikingly, a number of TET genes were expressed in embryonic and seedling progenitor cells and remained expressed until the differentiation state in the mature plant, suggesting a dynamic function over developmental stages. The cis-regulatory elements together with transcription factor-binding data provided molecular insight into the sites, conditions, and perturbations that affect TET gene expression and positioned the TET genes in different molecular pathways; the data represent a hypothesis-generating resource for further functional analyses.

BLSSpeller: exhaustive comparative discovery of conserved cis-regulatory elements.

MOTIVATION: The accurate discovery and annotation of regulatory elements remains a challenging problem. The growing number of sequenced genomes creates new opportunities for comparative approaches to motif discovery. Putative binding sites are then considered to be functional if they are conserved in orthologous promoter sequences of multiple related species. Existing methods for comparative motif discovery usually rely on pregenerated multiple sequence alignments, which are difficult to obtain for more diverged species such as plants. As a consequence, misaligned regulatory elements often remain undetected. RESULTS: We present a novel algorithm that supports both alignment-free and alignment-based motif discovery in the promoter sequences of related species. Putative motifs are exhaustively enumerated as words over the IUPAC alphabet and screened for conservation using the branch length score. Additionally, a confidence score is established in a genome-wide fashion. In order to take advantage of a cloud computing infrastructure, the MapReduce programming model is adopted. The method is applied to four monocotyledon plant species and it is shown that high-scoring motifs are significantly enriched for open chromatin regions in Oryza sativa and for transcription factor binding sites inferred through protein-binding microarrays in O.sativa and Zea mays. Furthermore, the method is shown to recover experimentally profiled ga2ox1-like KN1 binding sites in Z.mays. AVAILABILITY AND IMPLEMENTATION: BLSSpeller was written in Java. Source code and manual are available at http://bioinformatics.intec.ugent.be/blsspeller CONTACT: Klaas.Vandepoele@psb.vib-ugent.be or jan.fostier@intec.ugent.be. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Selection for Improved Energy Use Efficiency and Drought Tolerance in Canola Results in Distinct Transcriptome and Epigenome Changes.

To increase both the yield potential and stability of crops, integrated breeding strategies are used that have mostly a direct genetic basis, but the utility of epigenetics to improve complex traits is unclear. A better understanding of the status of the epigenome and its contribution to agronomic performance would help in developing approaches to incorporate the epigenetic component of complex traits into breeding programs. Starting from isogenic canola (Brassica napus) lines, epilines were generated by selecting, repeatedly for three generations, for increased energy use efficiency and drought tolerance. These epilines had an enhanced energy use efficiency, drought tolerance, and nitrogen use efficiency. Transcriptome analysis of the epilines and a line selected for its energy use efficiency solely revealed common differentially expressed genes related to the onset of stress tolerance-regulating signaling events. Genes related to responses to salt, osmotic, abscisic acid, and drought treatments were specifically differentially expressed in the drought-tolerant epilines. The status of the epigenome, scored as differential trimethylation of lysine-4 of histone 3, further supported the phenotype by targeting drought-responsive genes and facilitating the transcription of the differentially expressed genes. From these results, we conclude that the canola epigenome can be shaped by selection to increase energy use efficiency and stress tolerance. Hence, these findings warrant the further development of strategies to incorporate epigenetics into breeding.

Discrimination within epitope specific antibody populations against Classical swine fever virus is a new means of differentiating infection from vaccination.

Serological differentiation between infection and vaccination depends on the detection of pathogen specific antibodies for an epitope that is modified or lacking in a vaccine. Here we describe a new assay principle that is based on differences in the binding properties of epitope specific antibodies. C-DIVA is a potent Classical swine fever vaccine candidate that differs from the parental C-strain life attenuated vaccine in the highly immunogenic TAVSPTTLR epitope by the deletion of two and the mutation of one amino acid (TAGSΔΔTLR). We show that C-DIVA vaccination elicits antibodies with high affinity for both the TAGSΔΔTLR and TAVSPTTLR epitope, whereas infection elicits only TAVSPTTLR specific antibodies. Differentiation is achieved with a double competition assay with negative selection for antibodies with affinity for the TAGSΔΔTLR epitope followed by positive selection for antibodies with affinity for the TAVSPTTLR epitope. Our findings add a new strategy for the development of marker vaccines and their accompanying discrimination assays and offer an alternative to the devastating stamping out policy for Classical swine fever.

Assessing the immunomodulatory potential of high-molecular-weight extracts from mushrooms; an assay based on THP-1 macrophages.

BACKGROUND: Food is a potential source of immunomodulating compounds that may be used to steer immune responses towards a desired status such as reducing inflammatory disorders. However, to identify and characterize such bioactive compounds, biologically relevant and standardized assays are required. Macrophages play an important role in immunomodulation and are suited for developing cell-based assays. An assay was developed based on macrophages, in a homogeneous differentiation state, using the human monocytic cell line THP-1 previously used to assess immunomodulatory properties of low-molecular-weight allergens, hormones, dietary supplements and therapeutic drugs. RESULTS: Zymosan and mushroom polysaccharide extracts lead to a heterogeneous differentiation state of THP-1 monocytes, and these cells secrete low levels of cytokines upon stimulation. Differentiation into macrophages using a low concentration of phorbol 12-myristate 13-acetate improved responsiveness. Elevated levels of cytokines were secreted by cells in a homogenous differentiation state. In addition, it was determined that the assay performs best when using cells at a concentration of (2.5-5) × 10(5) cells mL(-1). CONCLUSION: An assay was developed suitable to distinguish the immunomodulatory properties of food compounds in a reproducible manner. It was evaluated using eight mushroom species by measuring the secretion of relevant cytokines TNF-α, IL-1ß, IL-6 and IL-10.

A functional and evolutionary perspective on transcription factor binding in Arabidopsis thaliana.

Understanding the mechanisms underlying gene regulation is paramount to comprehend the translation from genotype to phenotype. The two are connected by gene expression, and it is generally thought that variation in transcription factor (TF) function is an important determinant of phenotypic evolution. We analyzed publicly available genome-wide chromatin immunoprecipitation experiments for 27 TFs in Arabidopsis thaliana and constructed an experimental network containing 46,619 regulatory interactions and 15,188 target genes. We identified hub targets and highly occupied target (HOT) regions, which are enriched for genes involved in development, stimulus responses, signaling, and gene regulatory processes in the currently profiled network. We provide several lines of evidence that TF binding at plant HOT regions is functional, in contrast to that in animals, and not merely the result of accessible chromatin. HOT regions harbor specific DNA motifs, are enriched for differentially expressed genes, and are often conserved across crucifers and dicots, even though they are not under higher levels of purifying selection than non-HOT regions. Distal bound regions are under purifying selection as well and are enriched for a chromatin state showing regulation by the Polycomb repressive complex. Gene expression complexity is positively correlated with the total number of bound TFs, revealing insights in the regulatory code for genes with different expression breadths. The integration of noncanonical and canonical DNA motif information yields new hypotheses on cobinding and tethering between specific TFs involved in flowering and light regulation.

Inference of transcriptional networks in Arabidopsis through conserved noncoding sequence analysis.

Transcriptional regulation plays an important role in establishing gene expression profiles during development or in response to (a)biotic stimuli. Transcription factor binding sites (TFBSs) are the functional elements that determine transcriptional activity, and the identification of individual TFBS in genome sequences is a major goal to inferring regulatory networks. We have developed a phylogenetic footprinting approach for the identification of conserved noncoding sequences (CNSs) across 12 dicot plants. Whereas both alignment and non-alignment-based techniques were applied to identify functional motifs in a multispecies context, our method accounts for incomplete motif conservation as well as high sequence divergence between related species. We identified 69,361 footprints associated with 17,895 genes. Through the integration of known TFBS obtained from the literature and experimental studies, we used the CNSs to compile a gene regulatory network in Arabidopsis thaliana containing 40,758 interactions, of which two-thirds act through binding events located in DNase I hypersensitive sites. This network shows significant enrichment toward in vivo targets of known regulators, and its overall quality was confirmed using five different biological validation metrics. Finally, through the integration of detailed expression and function information, we demonstrate how static CNSs can be converted into condition-dependent regulatory networks, offering opportunities for regulatory gene annotation.

A DNA-binding-site landscape and regulatory network analysis for NAC transcription factors in Arabidopsis thaliana.

Target gene identification for transcription factors is a prerequisite for the systems wide understanding of organismal behaviour. NAM-ATAF1/2-CUC2 (NAC) transcription factors are amongst the largest transcription factor families in plants, yet limited data exist from unbiased approaches to resolve the DNA-binding preferences of individual members. Here, we present a TF-target gene identification workflow based on the integration of novel protein binding microarray data with gene expression and multi-species promoter sequence conservation to identify the DNA-binding specificities and the gene regulatory networks of 12 NAC transcription factors. Our data offer specific single-base resolution fingerprints for most TFs studied and indicate that NAC DNA-binding specificities might be predicted from their DNA-binding domain's sequence. The developed methodology, including the application of complementary functional genomics filters, makes it possible to translate, for each TF, protein binding microarray data into a set of high-quality target genes. With this approach, we confirm NAC target genes reported from independent in vivo analyses. We emphasize that candidate target gene sets together with the workflow associated with functional modules offer a strong resource to unravel the regulatory potential of NAC genes and that this workflow could be used to study other families of transcription factors.

3D domain swapping causes extensive multimerisation of human interleukin-10 when expressed in planta.

Heterologous expression platforms of biopharmaceutical proteins have been significantly improved over the last decade. Further improvement can be established by examining the intrinsic properties of proteins. Interleukin-10 (IL-10) is an anti-inflammatory cytokine with a short half-life that plays an important role in re-establishing immune homeostasis. This homodimeric protein of 36 kDa has significant therapeutic potential to treat inflammatory and autoimmune diseases. In this study we show that the major production bottleneck of human IL-10 is not protein instability as previously suggested, but extensive multimerisation due to its intrinsic 3D domain swapping characteristic. Extensive multimerisation of human IL-10 could be visualised as granules in planta. On the other hand, mouse IL-10 hardly multimerised, which could be largely attributed to its glycosylation. By introducing a short glycine-serine-linker between the fourth and fifth alpha helix of human IL-10 a stable monomeric form of IL-10 (hIL-10(mono)) was created that no longer multimerised and increased yield up to 20-fold. However, hIL-10(mono) no longer had the ability to reduce pro-inflammatory cytokine secretion from lipopolysaccharide-stimulated macrophages. Forcing dimerisation restored biological activity. This was achieved by fusing human IL-10(mono) to the C-terminal end of constant domains 2 and 3 of human immunoglobulin A (Fcα), a natural dimer. Stable dimeric forms of IL-10, like Fcα-IL-10, may not only be a better format for improved production, but also a more suitable format for medical applications.