Regulatory network identification by genetical genomics: signaling downstream of the Arabidopsis receptor-like kinase ERECTA

Plasticity of maternal environment-dependent expression-QTLs of tomato seeds

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

BREVIPEDICELLUS and ERECTA control the expression of AtPRX17 to prevent Arabidopsis callus browning

Efficient in vitro callus generation is required for tissue culture propagation, a process that allows for plant regeneration and transgenic breeding for desired phenotypes. Identifying genes and regulatory elements that prevent impaired callus growth and callus browning is essential for the development of in vitro callus systems. Here, we show that the BREVIPEDICELLUS and ERECTA pathways in Arabidopsis calli converge to prevent callus browning, and positively regulate the expression of the isoperoxidase gene AtPRX17 in rapidly growing calli. Loss-of-function mutations in both BREVIPEDICELLUS and ERECTA resulted in markedly increased callus browning. Transgenic lines expressing 35S pro::AtPRX17 in the bp-5 er105 double mutant background fully rescued this phenotypic abnormality. Using in vivo (chromatin immunoprecipitation-PCR and transient expression) and in vitro (electrophoretic mobility shift assays) assays, we observed that the BREVIPEDICELLUS protein binds directly to the upstream sequence of AtPRX17 to promote its transcription during callus growth. ERECTA is a ubiquitous factor required for cell proliferation and growth. We show that ERECTA positively regulates the expression of the transcription factor WRKY6, which directly binds to a separate site on the AtPRX17 promoter, further increasing its expression. Our data reveal an important molecular mechanism involved in the regulation of peroxidase isozyme expression to reduce Arabidopsis callus browning.

Network Analysis Prioritizes DEWAX and ICE1 as the Candidate Genes for Major eQTL Hotspots in Seed Germination of Arabidopsis thaliana

Seed germination is characterized by a constant change of gene expression across different time points. These changes are related to specific processes, which eventually determine the onset of seed germination. To get a better understanding on the regulation of gene expression during seed germination, we performed a quantitative trait locus mapping of gene expression (eQTL) at four important seed germination stages (primary dormant, after-ripened, six-hour after imbibition, and radicle protrusion stage) using Arabidopsis thaliana Bay x Sha recombinant inbred lines (RILs). The mapping displayed the distinctness of the eQTL landscape for each stage. We found several eQTL hotspots across stages associated with the regulation of expression of a large number of genes. Interestingly, an eQTL hotspot on chromosome five collocates with hotspots for phenotypic and metabolic QTL in the same population. Finally, we constructed a gene co-expression network to prioritize the regulatory genes for two major eQTL hotspots. The network analysis prioritizes transcription factors DEWAX and ICE1 as the most likely regulatory genes for the hotspot. Together, we have revealed that the genetic regulation of gene expression is dynamic along the course of seed germination.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

WormQTL2: an interactive platform for systems genetics in Caenorhabditis elegans

Quantitative genetics provides the tools for linking polymorphic loci to trait variation. Linkage analysis of gene expression is an established and widely applied method, leading to the identification of expression quantitative trait loci (eQTLs). (e)QTL detection facilitates the identification and understanding of the underlying molecular components and pathways, yet (e)QTL data access and mining often is a bottleneck. Here, we present WormQTL2, a database and platform for comparative investigations and meta-analyses of published (e)QTL data sets in the model nematode worm C. elegans. WormQTL2 integrates six eQTL studies spanning 11 conditions as well as over 1000 traits from 32 studies and allows experimental results to be compared, reused and extended upon to guide further experiments and conduct systems-genetic analyses. For example, one can easily screen a locus for specific cis-eQTLs that could be linked to variation in other traits, detect gene-by-environment interactions by comparing eQTLs under different conditions, or find correlations between QTL profiles of classical traits and gene expression. WormQTL2 makes data on natural variation in C. elegans and the identified QTLs interactively accessible, allowing studies beyond the original publications. Database URL: www.bioinformatics.nl/WormQTL2/.

Genome-wide cis-regulatory signatures for modulation of agronomic traits as exemplified by drought yield index (DYI) in chickpea

Developing functional molecular tags from the cis-regulatory sequence components of genes is vital for their deployment in efficient genetic dissection of complex quantitative traits in crop plants including chickpea. The current study identified 431,194 conserved non-coding SNP (CNSNP) from the cis-regulatory element regions of genes which were annotated on a chickpea genome. These genome-wide CNSNP marker resources are made publicly accessible through a user-friendly web-database ( http://www.cnsnpcicarbase.com ). The CNSNP-based quantitative trait loci (QTL) and expression QTL (eQTL) mapping and genome-wide association study (GWAS) were further integrated with global gene expression landscapes, molecular haplotyping, and DNA-protein interaction study in the association panel and recombinant inbred lines (RIL) mapping population to decode complex genetic architecture of one of the vital seed yield trait under drought stress, drought yield index (DYI), in chickpea. This delineated two constituted natural haplotypes and alleles from a histone H3 protein-coding gene and its transcriptional regulator NAC transcription factor (TF) harboring the major QTLs and trans-acting eQTL governing DYI in chickpea. The effect of CNSNPs in TF-binding cis-element of a histone H3 gene in altering the binding affinity and transcriptional activity of NAC TF based on chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR) assay was evident. The CNSNP-led promising molecular tags scanned will essentially have functional significance to decode transcriptional gene regulatory function and thus can drive translational genomic analysis in chickpea.

Genome-wide Analysis of Transcriptional Variability in a Large Maize-Teosinte Population

Gene expression regulation plays an important role in controlling plant phenotypes and adaptation. Here, we report a comprehensive assessment of gene expression variation through the transcriptome analyses of a large maize-teosinte experimental population. Genome-wide mapping identified 25 660 expression quantitative trait loci (eQTL) for 17 311 genes, capturing an unprecedented range of expression variation. We found that local eQTL were more frequently mapped to adjacent genes, displaying a mode of expression piggybacking, which consequently created co-regulated gene clusters. Genes within the co-regulated gene clusters tend to have relevant functions and shared chromatin modifications. Distant eQTL formed 125 significant distant eQTL hotspots with their targets significantly enriched in specific functional categories. By integrating different sources of information, we identified putative trans- regulators for a variety of metabolic pathways. We demonstrated that the bHLH transcription factor R1 and hexokinase HEX9 might act as crucial regulators for flavonoid biosynthesis and glycolysis, respectively. Moreover, we showed that domestication or improvement has significantly affected global gene expression, with many genes targeted by selection. Of particular interest, the Bx genes for benzoxazinoid biosynthesis may have undergone coordinated cis-regulatory divergence between maize and teosinte, and a transposon insertion that inactivates Bx12 was under strong selection as maize spread into temperate environments with a distinct herbivore community.

Construction of a High-Density Genetic Map from RNA-Seq Data for an Arabidopsis Bay-0 × Shahdara RIL Population

High-density genetic maps are essential for high resolution mapping of quantitative traits. Here, we present a new genetic map for an Arabidopsis Bayreuth × Shahdara recombinant inbred line (RIL) population, built on RNA-seq data. RNA-seq analysis on 160 RILs of this population identified 30,049 single-nucleotide polymorphisms (SNPs) covering the whole genome. Based on a 100-kbp window SNP binning method, 1059 bin-markers were identified, physically anchored on the genome. The total length of the RNA-seq genetic map spans 471.70 centimorgans (cM) with an average marker distance of 0.45 cM and a maximum marker distance of 4.81 cM. This high resolution genotyping revealed new recombination breakpoints in the population. To highlight the advantages of such high-density map, we compared it to two publicly available genetic maps for the same population, comprising 69 PCR-based markers and 497 gene expression markers derived from microarray data, respectively. In this study, we show that SNP markers can effectively be derived from RNA-seq data. The new RNA-seq map closes many existing gaps in marker coverage, saturating the previously available genetic maps. Quantitative trait locus (QTL) analysis for published phenotypes using the available genetic maps showed increased QTL mapping resolution and reduced QTL confidence interval using the RNA-seq map. The new high-density map is a valuable resource that facilitates the identification of candidate genes and map-based cloning approaches.

Quantitative trait loci from identification to exploitation for crop improvement

Advancement in the field of genetics and genomics after the discovery of Mendel's laws of inheritance has led to map the genes controlling qualitative and quantitative traits in crop plant species. Mapping of genomic regions controlling the variation of quantitatively inherited traits has become routine after the advent of different types of molecular markers. Recently, the next generation sequencing methods have accelerated the research on QTL analysis. These efforts have led to the identification of more closely linked molecular markers with gene/QTLs and also identified markers even within gene/QTL controlling the trait of interest. Efforts have also been made towards cloning gene/QTLs or identification of potential candidate genes responsible for a trait. Further new concepts like crop QTLome and QTL prioritization have accelerated precise application of QTLs for genetic improvement of complex traits. In the past years, efforts have also been made in exploitation of a number of QTL for improving grain yield or other agronomic traits in various crops through markers assisted selection leading to cultivation of these improved varieties at farmers' field. In present article, we reviewed QTLs from their identification to exploitation in plant breeding programs and also reviewed that how improved cultivars developed through introgression of QTLs have improved the yield productivity in many crops.

Contribution of trans regulatory eQTL to cryptic genetic variation in C. elegans

Background

Cryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population.

Results

We measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL.

Conclusions

These results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.

Electronic supplementary material

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

AraQTL - workbench and archive for systems genetics in Arabidopsis thaliana

Genetical genomics studies uncover genome-wide genetic interactions between genes and their transcriptional regulators. High-throughput measurement of gene expression in recombinant inbred line populations has enabled investigation of the genetic architecture of variation in gene expression. This has the potential to enrich our understanding of the molecular mechanisms affected by and underlying natural variation. Moreover, it contributes to the systems biology of natural variation, as a substantial number of experiments have resulted in a valuable amount of interconnectable phenotypic, molecular and genotypic data. A number of genetical genomics studies have been published for Arabidopsis thaliana, uncovering many expression quantitative trait loci (eQTLs). However, these complex data are not easily accessible to the plant research community, leaving most of the valuable genetic interactions unexplored as cross-analysis of these studies is a major effort. We address this problem with AraQTL (http://www.bioinformatics.nl/Ara QTL/), an easily accessible workbench and database for comparative analysis and meta-analysis of all published Arabidopsis eQTL datasets. AraQTL provides a workbench for comparing, re-using and extending upon the results of these experiments. For example, one can easily screen a physical region for specific local eQTLs that could harbour candidate genes for phenotypic QTLs, or detect gene-by-environment interactions by comparing eQTLs under different conditions.

Genotypes, Networks, Phenotypes: Moving Toward Plant Systems Genetics

One of the central goals in biology is to understand how and how much of the phenotype of an organism is encoded in its genome. Although many genes that are crucial for organismal processes have been identified, much less is known about the genetic bases underlying quantitative phenotypic differences in natural populations. We discuss the fundamental gap between the large body of knowledge generated over the past decades by experimental genetics in the laboratory and what is needed to understand the genotype-to-phenotype problem on a broader scale. We argue that systems genetics, a combination of systems biology and the study of natural variation using quantitative genetics, will help to address this problem. We present major advances in these two mostly disconnected areas that have increased our understanding of the developmental processes of flowering time control and root growth. We conclude by illustrating and discussing the efforts that have been made toward systems genetics specifically in plants.

Natural Genetic Variation Differentially Affects the Proteome and Transcriptome in Caenorhabditis elegans

Natural genetic variation is the raw material of evolution and influences disease development and progression. An important question is how this genetic variation translates into variation in protein abundance. To analyze the effects of the genetic background on gene and protein expression in the nematode Caenorhabditis elegans, we quantitatively compared the two genetically highly divergent wild-type strains N2 and CB4856. Gene expression was analyzed by microarray assays, and proteins were quantified using stable isotope labeling by amino acids in cell culture. Among all transcribed genes, we found 1,532 genes to be differentially transcribed between the two wild types. Of the total 3,238 quantified proteins, 129 proteins were significantly differentially expressed between N2 and CB4856. The differentially expressed proteins were enriched for genes that function in insulin-signaling and stress-response pathways, underlining strong divergence of these pathways in nematodes. The protein abundance of the two wild-type strains correlates more strongly than protein abundance versus transcript abundance within each wild type. Our findings indicate that in C. elegans only a fraction of the changes in protein abundance can be explained by the changes in mRNA abundance. These findings corroborate with the observations made across species.

Learning from Co-expression Networks: Possibilities and Challenges

Plants are fascinating and complex organisms. A comprehensive understanding of the organization, function and evolution of plant genes is essential to disentangle important biological processes and to advance crop engineering and breeding strategies. The ultimate aim in deciphering complex biological processes is the discovery of causal genes and regulatory mechanisms controlling these processes. The recent surge of omics data has opened the door to a system-wide understanding of the flow of biological information underlying complex traits. However, dealing with the corresponding large data sets represents a challenging endeavor that calls for the development of powerful bioinformatics methods. A popular approach is the construction and analysis of gene networks. Such networks are often used for genome-wide representation of the complex functional organization of biological systems. Network based on similarity in gene expression are called (gene) co-expression networks. One of the major application of gene co-expression networks is the functional annotation of unknown genes. Constructing co-expression networks is generally straightforward. In contrast, the resulting network of connected genes can become very complex, which limits its biological interpretation. Several strategies can be employed to enhance the interpretation of the networks. A strategy in coherence with the biological question addressed needs to be established to infer reliable networks. Additional benefits can be gained from network-based strategies using prior knowledge and data integration to further enhance the elucidation of gene regulatory relationships. As a result, biological networks provide many more applications beyond the simple visualization of co-expressed genes. In this study we review the different approaches for co-expression network inference in plants. We analyse integrative genomics strategies used in recent studies that successfully identified candidate genes taking advantage of gene co-expression networks. Additionally, we discuss promising bioinformatics approaches that predict networks for specific purposes.

On predicting regulatory genes by analysis of functional networks in C. elegans

Background

Connectivity networks, which reflect multiple interactions between genes and proteins, possess not only a descriptive but also a predictive value, as new connections can be extrapolated and tested by means of computational analysis. Integration of different types of connectivity data (such as co-expression and genetic interactions) in one network has proven to benefit ‘guilt by association’ analysis. However predictive values of connectives of different types, that had their specific functional meaning and topological characteristics were not obvious, and have been addressed in this analysis.

Methods

eQTL data for 3 experimental C.elegans age groups were retrieved from WormQTL. WormNet has been used to obtain pair-wise gene interactions. The Shortest Path Function (SPF) has been adopted for statistical validation of the co-expressed gene clusters and for computational prediction of their potential gene expression regulators from a network context. A new SPF-based algorithm has been applied to genetic interactions sub-networks adjacent to the clusters of co-expressed genes for ranking the most likely gene expression regulators causal to eQTLs.

Results

We have demonstrated that known co-expression and genetic interactions between C. elegans genes can be complementary in predicting gene expression regulators. Several algorithms were compared in respect to their predictive potential in different network connectivity contexts. We found that genes associated with eQTLs are highly clustered in a C. elegans co-expression sub-network, and their adjacent genetic interactions provide the optimal functional connectivity environment for application of the new SPF-based algorithm. It was successfully tested in the reverse-prediction analysis on groups of genes with known regulators and applied to co-expressed genes and experimentally observed expression quantitative trait loci (eQTLs).

Conclusions

This analysis demonstrates differences in topology and connectivity of co-expression and genetic interactions sub-networks in WormNet. The modularity of less continuous genetic interaction network does not correspond to modularity of the dense network comprised by gene co-expression interactions. However the genetic interaction network can be used much more efficiently with the SPF method in prediction of potential regulators of gene expression. The developed method can be used for validation of functional significance of suggested eQTLs and a discovery of new regulatory modules.

Genome-wide conserved non-coding microsatellite (CNMS) marker-based integrative genetical genomics for quantitative dissection of seed weight in chickpea

Phylogenetic footprinting identified 666 genome-wide paralogous and orthologous CNMS (conserved non-coding microsatellite) markers from 5'-untranslated and regulatory regions (URRs) of 603 protein-coding chickpea genes. The (CT)n and (GA)n CNMS carrying CTRMCAMV35S and GAGA8BKN3 regulatory elements, respectively, are abundant in the chickpea genome. The mapped genic CNMS markers with robust amplification efficiencies (94.7%) detected higher intraspecific polymorphic potential (37.6%) among genotypes, implying their immense utility in chickpea breeding and genetic analyses. Seventeen differentially expressed CNMS marker-associated genes showing strong preferential and seed tissue/developmental stage-specific expression in contrasting genotypes were selected to narrow down the gene targets underlying seed weight quantitative trait loci (QTLs)/eQTLs (expression QTLs) through integrative genetical genomics. The integration of transcript profiling with seed weight QTL/eQTL mapping, molecular haplotyping, and association analyses identified potential molecular tags (GAGA8BKN3 and RAV1AAT regulatory elements and alleles/haplotypes) in the LOB-domain-containing protein- and KANADI protein-encoding transcription factor genes controlling the cis-regulated expression for seed weight in the chickpea. This emphasizes the potential of CNMS marker-based integrative genetical genomics for the quantitative genetic dissection of complex seed weight in chickpea.

Molecular phenotyping of plant single cell-types enhances forward genetic analyses

Recent advances in the isolation of single cell-types in plants provides an opportunity to conduct detailed analyses of their molecular characteristics at high resolution. This kind of cell-type specific molecular phenotyping is likely to enhance forward genetics studies to dissect the effect of mutations and thereby aid gene function assignment. Recent experimental results support this view, demonstrating that different cell-types exhibit substantial variation in transcript, protein, and metabolite accumulation and these molecular phenotypes are often sensitive to genetic and environmental alterations. The use of single cell-type molecular phenotyping approach to define plant gene function is most amenable to cell-types with well-characterized molecular tools and isolation protocols.

Highly polygenic variation in environmental perception determines dauer larvae formation in growing populations of Caenorhabditis elegans

Background

Determining how complex traits are genetically controlled is a requirement if we are to predict how they evolve and how they might respond to selection. This requires understanding how distinct, and often more simple, life history traits interact and change in response to environmental conditions. In order to begin addressing such issues, we have been analyzing the formation of the developmentally arrested dauer larvae of Caenorhabditis elegans under different conditions.

Results

We find that 18 of 22 previously identified quantitative trait loci (QTLs) affecting dauer larvae formation in growing populations, assayed by determining the number of dauer larvae present at food patch exhaustion, can be recovered under various environmental conditions. We also show that food patch size affects both the ability to detect QTLs and estimates of effect size, and demonstrate that an allele of nath-10 affects dauer larvae formation in growing populations. To investigate the component traits that affect dauer larvae formation in growing populations we map, using the same introgression lines, QTLs that affect dauer larvae formation in response to defined amounts of pheromone. This identifies 36 QTLs, again demonstrating the highly polygenic nature of the genetic variation underlying dauer larvae formation.

Conclusions

These data indicate that QTLs affecting the number of dauer larvae at food exhaustion in growing populations of C. elegans are highly reproducible, and that nearly all can be explained by variation affecting dauer larvae formation in response to defined amounts of pheromone. This suggests that most variation in dauer larvae formation in growing populations is a consequence of variation in the perception of the food and pheromone environment (i.e. chemosensory variation) and in the integration of these cues.

Regulatory network rewiring for secondary metabolism in Arabidopsis thaliana under various conditions

BACKGROUND

Plant secondary metabolites are critical to various biological processes. However, the regulations of these metabolites are complex because of regulatory rewiring or crosstalk. To unveil how regulatory behaviors on secondary metabolism reshape biological processes, we constructed and analyzed a dynamic regulatory network of secondary metabolic pathways in Arabidopsis.

RESULTS

The dynamic regulatory network was constructed through integrating co-expressed gene pairs and regulatory interactions. Regulatory interactions were either predicted by conserved transcription factor binding sites (TFBSs) or proved by experiments. We found that integrating two data (co-expression and predicted regulatory interactions) enhanced the number of highly confident regulatory interactions by over 10% compared with using single data. The dynamic changes of regulatory network systematically manifested regulatory rewiring to explain the mechanism of regulation, such as in terpenoids metabolism, the regulatory crosstalk of RAV1 (AT1G13260) and ATHB1 (AT3G01470) on HMG1 (hydroxymethylglutaryl-CoA reductase, AT1G76490); and regulation of RAV1 on epoxysqualene biosynthesis and sterol biosynthesis. Besides, we investigated regulatory rewiring with expression, network topology and upstream signaling pathways. Regulatory rewiring was revealed by the variability of genes' expression: pathway genes and transcription factors (TFs) were significantly differentially expressed under different conditions (such as terpenoids biosynthetic genes in tissue experiments and E2F/DP family members in genotype experiments). Both network topology and signaling pathways supported regulatory rewiring. For example, we discovered correlation among the numbers of pathway genes, TFs and network topology: one-gene pathways (such as δ-carotene biosynthesis) were regulated by a fewer TFs, and were not critical to metabolic network because of their low degrees in topology. Upstream signaling pathways of 50 TFs were identified to comprehend the underlying mechanism of TFs' regulatory rewiring.

CONCLUSION

Overall, this dynamic regulatory network largely improves the understanding of perplexed regulatory rewiring in secondary metabolism in Arabidopsis.

Loss-of-function of β-catenin bar-1 slows development and activates the Wnt pathway in Caenorhabditis elegans

C. elegans is extensively used to study the Wnt-pathway and most of the core-signalling components are known. Four β-catenins are important gene expression regulators in Wnt-signalling. One of these, bar-1, is part of the canonical Wnt-pathway. Together with Wnt effector pop-1, bar-1 forms a transcription activation complex which regulates the transcription of downstream genes. The effects of bar-1 loss-of-function mutations on many phenotypes have been studied well. However, the effects on global gene expression are unknown. Here we report the effects of a loss-of-function mutation bar-1(ga80). By analysing the transcriptome and developmental phenotyping we show that bar-1(ga80) impairs developmental timing. This developmental difference confounds the comparison of the gene expression profile between the mutant and the reference strain. When corrected for this difference it was possible to identify genes that were directly affected by the bar-1 mutation. We show that the Wnt-pathway itself is activated, as well as transcription factors elt-3, pqm-1, mdl-1 and pha-4 and their associated genes. The outcomes imply that this response compensates for the loss of functional bar-1. Altogether we show that bar-1 loss-of function leads to delayed development possibly caused by an induction of a stress response, reflected by daf-16 activated genes.

An expression quantitative trait loci-guided co-expression analysis for constructing regulatory network using a rice recombinant inbred line population

The ability to reveal the regulatory architecture of genes at the whole-genome level by constructing a regulatory network is critical for understanding the biological processes and developmental programmes of organisms. Here, we conducted an eQTL-guided function-related co-expression analysis to identify the putative regulators and construct gene regulatory network. We performed an eQTL analysis of 210 recombinant inbred lines (RILs) derived from a cross between two indica rice lines, Zhenshan 97 and Minghui 63, the parents of an elite hybrid, using data obtained by hybridizing RNA samples of flag leaves at the heading stage with Affymetrix whole-genome arrays. Making use of an ultrahigh-density single-nucleotide polymorphism bin map constructed by population sequencing, 13 647 eQTLs for 10 725 e-traits were detected, comprising 5079 cis-eQTLs (37.2%) and 8568 trans-eQTLs (62.8%). The analysis revealed 138 trans-eQTLs hotspots, each of which apparently regulates the expression variations of many genes. Co-expression analysis of functionally related genes within the framework of regulator-target relationships outlined by the eQTLs led to the identification of putative regulators in the system. The usefulness of the strategy was demonstrated with the genes known to be involved in flowering. We also applied this strategy to the analysis of QTLs for yield traits, which also suggested likely candidate genes. eQTL-guided co-expression analysis may provide a promising solution for outlining a framework for the complex regulatory network of an organism.

Diverse roles of ERECTA family genes in plant development

Multiple receptor-like kinases (RLKs) enable intercellular communication that coordinates growth and development of plant tissues. ERECTA family receptors (ERfs) are an ancient family of leucine-rich repeat RLKs that in Arabidopsis consists of three genes: ERECTA, ERL1, and ERL2. ERfs sense secreted cysteine-rich peptides from the EPF/EPFL family and transmit the signal through a MAP kinase cascade. This review discusses the functions of ERfs in stomata development, in regulation of longitudinal growth of aboveground organs, during reproductive development, and in the shoot apical meristem. In addition the role of ERECTA in plant responses to biotic and abiotic factors is examined. Elena D. Shpak (Corresponding author).

WormQTLHD--a web database for linking human disease to natural variation data in C. elegans

Interactions between proteins are highly conserved across species. As a result, the molecular basis of multiple diseases affecting humans can be studied in model organisms that offer many alternative experimental opportunities. One such organism—Caenorhabditis elegans—has been used to produce much molecular quantitative genetics and systems biology data over the past decade. We present WormQTL^HD (Human Disease), a database that quantitatively and systematically links expression Quantitative Trait Loci (eQTL) findings in C. elegans to gene–disease associations in man. WormQTL^HD, available online at http://www.wormqtl-hd.org, is a user-friendly set of tools to reveal functionally coherent, evolutionary conserved gene networks. These can be used to predict novel gene-to-gene associations and the functions of genes underlying the disease of interest. We created a new database that links C. elegans eQTL data sets to human diseases (34 337 gene–disease associations from OMIM, DGA, GWAS Central and NHGRI GWAS Catalogue) based on overlapping sets of orthologous genes associated to phenotypes in these two species. We utilized QTL results, high-throughput molecular phenotypes, classical phenotypes and genotype data covering different developmental stages and environments from WormQTL database. All software is available as open source, built on MOLGENIS and xQTL workbench.

Genetic mapping of variation in dauer larvae development in growing populations of Caenorhabditis elegans

In the nematode Caenorhabditis elegans, the appropriate induction of dauer larvae development within growing populations is likely to be a primary determinant of genotypic fitness. The underlying genetic architecture of natural genetic variation in dauer formation has, however, not been thoroughly investigated. Here, we report extensive natural genetic variation in dauer larvae development within growing populations across multiple wild isolates. Moreover, bin mapping of introgression lines (ILs) derived from the genetically divergent isolates N2 and CB4856 reveals 10 quantitative trait loci (QTLs) affecting dauer formation. Comparison of individual ILs to N2 identifies an additional eight QTLs, and sequential IL analysis reveals six more QTLs. Our results also show that a behavioural, laboratory-derived, mutation controlled by the neuropeptide Y receptor homolog npr-1 can affect dauer larvae development in growing populations. These findings illustrate the complex genetic architecture of variation in dauer larvae formation in C. elegans and may help to understand how the control of variation in dauer larvae development has evolved.

Tissue-specific profiling reveals transcriptome alterations in Arabidopsis mutants lacking morphological phenotypes

Traditional genetic analysis relies on mutants with observable phenotypes. Mutants lacking visible abnormalities may nevertheless exhibit molecular differences useful for defining gene function. To examine this, we analyzed tissue-specific transcript profiles from Arabidopsis thaliana transcription factor gene mutants with known roles in root epidermis development, but lacking a single-gene mutant phenotype due to genetic redundancy. We discovered substantial transcriptional changes in each mutant, preferentially affecting root epidermal genes in a manner consistent with the known double mutant effects. Furthermore, comparing transcript profiles of single and double mutants, we observed remarkable variation in the sensitivity of target genes to the loss of one or both paralogous genes, including preferential effects on specific branches of the epidermal gene network, likely reflecting the pathways of paralog subfunctionalization during evolution. In addition, we analyzed the root epidermal transcriptome of the transparent testa glabra2 mutant to clarify its role in the network. These findings provide insight into the molecular basis of genetic redundancy and duplicate gene diversification at the level of a specific gene regulatory network, and they demonstrate the usefulness of tissue-specific transcript profiling to define gene function in mutants lacking informative visible changes in phenotype.

Constitutively elevated salicylic acid levels alter photosynthesis and oxidative state but not growth in transgenic populus

Salicylic acid (SA) has long been implicated in plant responses to oxidative stress. SA overproduction in Arabidopsis thaliana leads to dwarfism, making in planta assessment of SA effects difficult in this model system. We report that transgenic Populus tremula × alba expressing a bacterial SA synthase hyperaccumulated SA and SA conjugates without negative growth consequences. In the absence of stress, endogenously elevated SA elicited widespread metabolic and transcriptional changes that resembled those of wild-type plants exposed to oxidative stress-promoting heat treatments. Potential signaling and oxidative stress markers azelaic and gluconic acids as well as antioxidant chlorogenic acids were strongly coregulated with SA, while soluble sugars and other phenylpropanoids were inversely correlated. Photosynthetic responses to heat were attenuated in SA-overproducing plants. Network analysis identified potential drivers of SA-mediated transcriptome rewiring, including receptor-like kinases and WRKY transcription factors. Orthologs of Arabidopsis SA signaling components NON-EXPRESSOR OF PATHOGENESIS-RELATED GENES1 and thioredoxins were not represented. However, all members of the expanded Populus nucleoredoxin-1 family exhibited increased expression and increased network connectivity in SA-overproducing Populus, suggesting a previously undescribed role in SA-mediated redox regulation. The SA response in Populus involved a reprogramming of carbon uptake and partitioning during stress that is compatible with constitutive chemical defense and sustained growth, contrasting with the SA response in Arabidopsis, which is transient and compromises growth if sustained.

Temperature-dependent shade avoidance involves the receptor-like kinase ERECTA

Plants detect the presence of neighbouring vegetation by monitoring changes in the ratio of red (R) to far-red (FR) wavelengths (R:FR) in ambient light. Reductions in R:FR are perceived by the phytochrome family of plant photoreceptors and initiate a suite of developmental responses termed the shade avoidance syndrome. These include increased elongation growth of stems and petioles, enabling plants to overtop competing vegetation. The majority of shade avoidance experiments are performed at standard laboratory growing temperatures (>20°C). In these conditions, elongation responses to low R:FR are often accompanied by reductions in leaf development and accumulation of plant biomass. Here we investigated shade avoidance responses at a cooler temperature (16°C). In these conditions, Arabidopsis thaliana displays considerable low R:FR-mediated increases in leaf area, with reduced low R:FR-mediated petiole elongation and leaf hyponasty responses. In Landsberg erecta, these strikingly different shade avoidance phenotypes are accompanied by increased leaf thickness, increased biomass and an altered metabolite profile. At 16°C, low R:FR treatment results in the accumulation of soluble sugars and metabolites associated with cold acclimation. Analyses of natural genetic variation in shade avoidance responses at 16°C have revealed a regulatory role for the receptor-like kinase ERECTA.

Natural variation identifies multiple loci controlling petal shape and size in Arabidopsis thaliana

Natural variation in organ morphologies can have adaptive significance and contribute to speciation. However, the underlying allelic differences responsible for variation in organ size and shape remain poorly understood. We have utilized natural phenotypic variation in three Arabidopsis thaliana ecotypes to examine the genetic basis for quantitative variation in petal length, width, area, and shape. We identified 23 loci responsible for such variation, many of which appear to correspond to genes not previously implicated in controlling organ morphology. These analyses also demonstrated that allelic differences at distinct loci can independently affect petal length, width, area or shape, suggesting that these traits behave as independent modules. We also showed that ERECTA (ER), encoding a leucine-rich repeat (LRR) receptor-like serine-threonine kinase, is a major effect locus determining petal shape. Allelic variation at the ER locus was associated with differences in petal cell proliferation and concomitant effects on petal shape. ER has been previously shown to be required for regulating cell division and expansion in other contexts; the ER receptor-like kinase functioning to also control organ-specific proliferation patterns suggests that allelic variation in common signaling components may nonetheless have been a key factor in morphological diversification.

Genetical Genomics Reveals Large Scale Genotype-By-Environment Interactions in Arabidopsis thaliana

One of the major goals of quantitative genetics is to unravel the complex interactions between molecular genetic factors and the environment. The effects of these genotype-by-environment interactions also affect and cause variation in gene expression. The regulatory loci responsible for this variation can be found by genetical genomics that involves the mapping of quantitative trait loci (QTLs) for gene expression traits also called expression-QTL (eQTLs). Most genetical genomics experiments published so far, are performed in a single environment and hence do not allow investigation of the role of genotype-by-environment interactions. Furthermore, most studies have been done in a steady state environment leading to acclimated expression patterns. However a response to the environment or change therein can be highly plastic and possibly lead to more and larger differences between genotypes. Here we present a genetical genomics study on 120 Arabidopsis thaliana, Landsberg erecta × Cape Verde Islands, recombinant inbred lines (RILs) in active response to the environment by treating them with 3 h of shade. The results of this experiment are compared to a previous study on seedlings of the same RILs from a steady state environment. The combination of two highly different conditions but exactly the same RILs with a fixed genetic variation showed the large role of genotype-by-environment interactions on gene expression levels. We found environment-dependent hotspots of transcript regulation. The major hotspot was confirmed by the expression profile of a near isogenic line. Our combined analysis leads us to propose CSN5A, a COP9 signalosome component, as a candidate regulator for the gene expression response to shade.

Sensing the environment: key roles of membrane-localized kinases in plant perception and response to abiotic stress

Adverse environmental conditions have negative effects on plant growth and development. Receptor proteins on the plasma membrane sense various environmental stimuli and transduce them to downstream intra- and intercellular signalling networks. Receptor-like kinases (RLKs) play important roles in perceiving the extracellular ligands and activating the downstream pathway via phosphorylation of intracellular serine/threonine kinase domains. The Arabidopsis genome possesses >600 RLK-encoding genes, some of which are implicated in the perception of environmental signals during the life cycle of the sessile plants. Histidine kinases are also membrane-localized kinases and perceive osmotic stress and plant hormones. In this review, we focus on the RLKs and histidine kinases that play a role in plant response to abiotic stresses. We summarize our recent understanding of their specific roles in stress responses and absicisic acid (ABA) regulation. Elucidation of the functions of these kinases in the osmotic stress response will provide a better understanding of stress-sensing mechanisms in plants and help to identify potential candidate genes for genetic engineering of improved stress-tolerant crops.

Aging Uncouples Heritability and Expression-QTL in Caenorhabditis elegans

The number and distribution of gene expression QTL (eQTL) represent the genetic architecture of many complex traits, including common human diseases. We previously reported that the heritable eQTL patterns are highly dynamic with age in an N2 × CB4856 recombinant inbred population of the nematode Caenorhabditis elegans. In particular, we showed that the number of eQTL decreased with age. Here, we investigated the reason for this decrease by combining gene expression profiles at three ages in the wild types N2 and CB4856 with the reported expression profiles of the RIL population. We determined heritability and transgression (when gene expression levels in the RILs are more extreme than the parents) and investigated their relation with eQTL changes with age. Transgressive segregation was widespread but depended on physiological age. The percentage of genes with an eQTL increased with a higher heritability in young worms. However, for old worms this percentage hardly increased. Using a single marker approach, we found that almost 20% of genes with heritability >0.9 had an eQTL in developing worms. Surprisingly, only 10% was found in old worms. Using a multimarker approach, this percentage increased to almost 30% for both age groups. Comparison of the single marker to a multiple marker eQTL mapping indicated that heritable regulation of gene expression becomes more polygenic in aging worms due to multiple loci and possible epistatic interactions. We conclude that linkage studies should account for the relation between increased polygenic regulation and diminished effects at older ages.

Regulation of inflorescence architecture by intertissue layer ligand-receptor communication between endodermis and phloem

Multicellular organisms achieve final body shape and size by coordinating cell proliferation, expansion, and differentiation. Loss of function in the Arabidopsis ERECTA (ER) receptor-kinase gene confers characteristic compact inflorescence architecture, but its underlying signaling pathways remain unknown. Here we report that the expression of ER in the phloem is sufficient to rescue compact er inflorescences. We further identified two Epidermal Patterning Factor-like (EPFL) secreted peptide genes, EPFL4 and EPFL6/CHALLAH (CHAL), as redundant, upstream components of ER-mediated inflorescence growth. The expression of EPFL4 or EPFL6 in the endodermis, a layer adjacent to phloem, is sufficient to rescue the er-like inflorescence of epfl4 epfl6 plants. EPFL4 and EPFL6 physically associate with ER in planta. Finally, transcriptome analysis of er and epfl4 epfl6 revealed a potential downstream component as well as a role for plant hormones in EPFL4/6- and ER-mediated inflorescence growth. Our results suggest that intercell layer communication between the endodermis and phloem mediated by peptide ligands and a receptor kinase coordinates proper inflorescence architecture in Arabidopsis.

Lessons from eQTL mapping studies: non-coding regions and their role behind natural phenotypic variation in plants

Even if considerable progress has been achieved towards the understanding of natural variation in plant systems, the contribution of transcript abundance variation to phenotypic diversity remains unappreciated. Over the last decade, efforts to characterise the genome-wide expression variation in natural accessions, structured populations and hybrids have improved our knowledge of the contribution of non-coding polymorphisms to gene expression regulation. Moreover, new studies are helping to unravel the role of expression polymorphisms and their orchestrated performance. Recent advances involving classical linkage analysis, GWAS and improved eQTL mapping strategies will provide a greater resolution to determine the genetic variants shaping the broad diversity in plant systems.

Organ specificity and transcriptional control of metabolic routes revealed by expression QTL profiling of source--sink tissues in a segregating potato population

BACKGROUND

With the completion of genome sequences belonging to some of the major crop plants, new challenges arise to utilize this data for crop improvement and increased food security. The field of genetical genomics has the potential to identify genes displaying heritable differential expression associated to important phenotypic traits. Here we describe the identification of expression QTLs (eQTLs) in two different potato tissues of a segregating potato population and query the potato genome sequence to differentiate between cis- and trans-acting eQTLs in relation to gene subfunctionalization.

RESULTS

Leaf and tuber samples were analysed and screened for the presence of conserved and tissue dependent eQTLs. Expression QTLs present in both tissues are predominantly cis-acting whilst for tissue specific QTLs, the percentage of trans-acting QTLs increases. Tissue dependent eQTLs were assigned to functional classes and visualized in metabolic pathways. We identified a potential regulatory network on chromosome 10 involving genes crucial for maintaining circadian rhythms and controlling clock output genes. In addition, we show that the type of genetic material screened and sampling strategy applied, can have a high impact on the output of genetical genomics studies.

CONCLUSIONS

Identification of tissue dependent regulatory networks based on mapped differential expression not only gives us insight in tissue dependent gene subfunctionalization but brings new insights into key biological processes and delivers targets for future haplotyping and genetic marker development.

Developmental genetics and new sequencing technologies: the rise of nonmodel organisms

Much of developmental biology in the past decades has been driven by forward genetic studies in a few model organisms. We review recent work with relatives of these species, motivated by a desire to understand the evolutionary and ecological context for morphological innovation. Unfortunately, despite a number of shining examples, progress in nonmodel systems has often been slow. The current revolution in DNA sequencing has, however, enormous potential in extending the reach of genetics. We discuss how developmental biology will benefit from these advances, particularly by increasing the universe of study species.

Arabidopsis ERECTA-family receptor kinases mediate morphological alterations stimulated by activation of NB-LRR-type UNI proteins

Shoot apical meristems (SAMs), which maintain stem cells at the tips of stems, and axillary meristems (AMs), which arise at leaf axils for branch formation, play significant roles in the establishment of plant architecture. Previously, we showed that, in Arabidopsis thaliana, activation of NB-LRR (nucleotide-binding site-leucine-rich repeat)-type UNI proteins affects plant morphology through modulation of the regulation of meristems. However, information about genes involved in the processes was still lacking. Here, we report that ERECTA (ER) receptor kinase family members cooperatively mediate the morphological alterations that are stimulated by activation of UNI proteins. uni-1D is a gain-of-function mutation in the UNI gene and uni-1D mutants exhibit early termination of inflorescence stem growth and also formation of extra AMs at leaf axils. The former defect involves modulation of the SAM activity and is suppressed by er mutation. Though the AM phenotype is not affected by a single er mutation, it is suppressed by simultaneous mutations of ER-family members. It was previously shown that trans-zeatin (tZ)-type cytokinins were involved in the morphological phenotypes of uni-1D mutants and that expression of CYP735A2, which is essential for biosynthesis of tZ-type cytokinins, was modulated in uni-1D mutants. We show that this modulation of CYP735A2 expression requires activities of ER-family members. Moreover, the ER activity in UNI-expressing cells contributes to all morphological phenotypes of uni-1D mutants, suggesting that a cross-talk between ER-family-dependent and UNI-triggered signaling pathways plays a significant role in the morphological alterations observed in uni-1D mutants.