Transcriptomic Analysis of Flowering Time Genes in Cultivated Chickpea and Wild Cicer.

Chickpea (Cicer arietinum L.) is a major grain legume and a good source of plant-based protein. However, comprehensive knowledge of flowering time control in Cicer is lacking. In this study, we acquire high-throughput transcriptome sequencing data and analyze changes in gene expression during floral transition in the early flowering cultivar ICCV 96029, later flowering C. arietinum accessions, and two wild species, C. reticulatum and C. echinospermum. We identify Cicer orthologs of A. thaliana flowering time genes and analyze differential expression of 278 genes between four species/accessions, three tissue types, and two conditions. Our results show that the differences in gene expression between ICCV 96029 and other cultivated chickpea accessions are vernalization-dependent. In addition, we highlight the role of FTa3, an ortholog of FLOWERING LOCUS T in Arabidopsis, in the vernalization response of cultivated chickpea. A common set of differentially expressed genes was found for all comparisons between wild species and cultivars. The direction of expression change for different copies of the FT-INTERACTING PROTEIN 1 gene was variable in different comparisons, which suggests complex mechanisms of FT protein transport. Our study makes a contribution to the understanding of flowering time control in Cicer, and can provide genetic strategies to further improve this important agronomic trait.

Modeling the Flowering Activation Motif during Vernalization in Legumes: A Case Study of M. trancatula.

In many plant species, flowering is promoted by the cold treatment or vernalization. The mechanism of vernalization-induced flowering has been extensively studied in Arabidopsis but remains largely unknown in legumes. The orthologs of the FLC gene, a major regulator of vernalization response in Arabidopsis, are absent or non-functional in the vernalization-sensitive legume species. Nevertheless, the legume integrator genes FT and SOC1 are involved in the transition of the vernalization signal to meristem identity genes, including PIM (AP1 ortholog). However, the regulatory contribution of these genes to PIM activation in legumes remains elusive. Here, we presented the theoretical and data-driven analyses of a feed-forward regulatory motif that includes a vernalization-responsive FT gene and several SOC1 genes, which independently activate PIM and thereby mediate floral transition. Our theoretical model showed that the multiple regulatory branches in this regulatory motif facilitated the elimination of no-sense signals and amplified useful signals from the upstream regulator. We further developed and analyzed four data-driven models of PIM activation in Medicago trancatula in vernalized and non-vernalized conditions in wild-type and fta1-1 mutants. The model with FTa1 providing both direct activation and indirect activation via three intermediate activators, SOC1a, SOC1b, and SOC1c, resulted in the most relevant PIM dynamics. In this model, the difference between regulatory inputs of SOC1 genes was nonessential. As a result, in the M. trancatula model, the cumulative action of SOC1a, SOC1b, and SOC1c was favored. Overall, in this study, we first presented the in silico analysis of vernalization-induced flowering in legumes. The considered vernalization network motif can be supplemented with additional regulatory branches as new experimental data become available.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration.

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

Mechanisms of Vernalization-Induced Flowering in Legumes.

Vernalization is the requirement for exposure to low temperatures to trigger flowering. The best knowledge about the mechanisms of vernalization response has been accumulated for Arabidopsis and cereals. In Arabidopsis thaliana, vernalization involves an epigenetic silencing of the MADS-box gene FLOWERING LOCUS C (FLC), which is a flowering repressor. FLC silencing releases the expression of the main flowering inductor FLOWERING LOCUS T (FT), resulting in a floral transition. Remarkably, no FLC homologues have been identified in the vernalization-responsive legumes, and the mechanisms of cold-mediated transition to flowering in these species remain elusive. Nevertheless, legume FT genes have been shown to retain the function of the main vernalization signal integrators. Unlike Arabidopsis, legumes have three subclades of FT genes, which demonstrate distinct patterns of regulation with respect to environmental cues and tissue specificity. This implies complex mechanisms of vernalization signal propagation in the flowering network, that remain largely elusive. Here, for the first time, we summarize the available information on the genetic basis of cold-induced flowering in legumes with a special focus on the role of FT genes.

The Genetic Landscape of Fiber Flax.

Genetic diversity in a breeding program is essential to overcome modern-day environmental challenges faced by humanity and produce robust, resilient crop cultivars with improved agronomic characteristics, as well as to trace crop domestication history. Flax (Linum usitatissimum), one of the first crops domesticated by mankind, has been traditionally cultivated for fiber as well as for medicinal purposes and as a nutritional product. The origins of fiber flax are hidden in the mists of time and can be hypothetically traced back to either the Indo-Afghan region or Fertile Crescent. To shed new light on fiber flax genetic diversity and breeding history, in this study, we presented a comprehensive analysis of the core collection of flax (306 accessions) of different morphotypes and geographic origins maintained by the Russian Federal Research Center for Bast Fiber Crops. We observed significant population differentiation between oilseed and fiber morphotypes, as well as mapped genomic regions affected by recent breeding efforts. We also sought to unravel the origins of kryazhs, Russian heritage landraces, and their genetic relatedness to modern fiber flax cultivars. For the first time, our results provide strong genetic evidence in favor of the hypothesis on kryazh's mixed origin from both the Indo-Afghan diversity center and Fertile Crescent. Finally, we showed predominant contribution from Russian landraces and kryazhs into the ancestry of modern fiber flax varieties. Taken together, these findings may have practical implications on the development of new improved flax varieties with desirable traits that give farmers greater choice in crop management and meet the aspirations of breeders.

Application of Domain- and Genotype-Specific Models to Infer Post-Transcriptional Regulation of Segmentation Gene Expression in Drosophila.

Unlike transcriptional regulation, the post-transcriptional mechanisms underlying zygotic segmentation gene expression in early Drosophila embryo have been insufficiently investigated. Condition-specific post-transcriptional regulation plays an important role in the development of many organisms. Our recent study revealed the domain- and genotype-specific differences between mRNA and the protein expression of Drosophila hb, gt, and eve genes in cleavage cycle 14A. Here, we use this dataset and the dynamic mathematical model to recapitulate protein expression from the corresponding mRNA patterns. The condition-specific nonuniformity in parameter values is further interpreted in terms of possible post-transcriptional modifications. For hb expression in wild-type embryos, our results predict the position-specific differences in protein production. The protein synthesis rate parameter is significantly higher in hb anterior domain compared to the posterior domain. The parameter sets describing Gt protein dynamics in wild-type embryos and Kr mutants are genotype-specific. The spatial discrepancy between gt mRNA and protein posterior expression in Kr mutants is well reproduced by the whole axis model, thus rejecting the involvement of post-transcriptional mechanisms. Our models fail to describe the full dynamics of eve expression, presumably due to its complex shape and the variable time delays between mRNA and protein patterns, which likely require a more complex model. Overall, our modeling approach enables the prediction of regulatory scenarios underlying the condition-specific differences between mRNA and protein expression in early embryo.

Interplay between sex determination cascade and major signaling pathways during Drosophila eye development: Perspectives for future research.

Understanding molecular mechanisms of sexually dimorphic organ growth is a fundamental problem of developmental biology. Recent quantitative studies showed that the Drosophila compound eye is a convenient model to study the determination of the final organ size. In Drosophila, females have larger eyes than males and this is evident even after correction for the larger body size. Moreover, female eyes include more ommatidia (photosensitive units) than male eyes and this difference is specified at the third larval instar in the eye primordia called eye imaginal discs. This may result in different visual capabilities between the two sexes and have behavioral consequences. Despite growing evidence on the genetic bases of eye size variation between different Drosophila species and strains, mechanisms responsible for within-species sexual dimorphism still remain elusive. Here, we discuss a presumptive crosstalk between the sex determination cascade and major signaling pathways during dimorphic eye development. Male- and female-specific isoforms of Doublesex (Dsx) protein are known to control sex-specific differentiation in the somatic tissues. However, no data on Dsx function during eye disc growth and patterning are currently available. Remarkably, Sex lethal (Sxl), the sex determination switch protein, was shown to directly affect Hedgehog (Hh) and Notch (N) signaling in the Drosophila wing disc. The similarity of signaling pathways involved in the wing and eye disc growth suggests that Sxl might be integrated into regulation of eye development. Dsx role in the eye disc requires further investigation. We discuss currently available data on sex-biased gene expression in the Drosophila eye and highlight perspectives for future studies.

Quantitative analysis reveals genotype- and domain- specific differences between mRNA and protein expression of segmentation genes in Drosophila.

In many biological systems gene expression at mRNA and protein levels is not identical. Rigorous comparison of such differences on a spatio-temporal scale is still not feasible by high-throughput transcriptomic and proteomic analyses of early embryo development. Here, we characterize differences between mRNA and protein expression of Drosophila segmentation genes at the level of individual gene expression domains. We obtained quantitative imaging data on expression of gap genes gt and hb and pair-rule gene eve for Drosophila wild type embryos, Kr null mutants and Kr+/Kr- heterozygotes. To compare mRNA and protein expression we use several criteria including difference in amplitude and positions of expression domains, pattern shape and positional variability. For a number of gene expression domains we show examples where protein expression does not repeat mRNA expression even after a temporal delay. We calculated time delays between eve pattern formation at the level of mRNA and protein for wild type embryos, Kr mutants and Kr+/Kr- heterozygotes. We detect that in wild type embryos, the amplitudes of eve stripes 3 and 7 do not differ significantly at the level of mRNA, however, stripe 3 is higher than stripe 7 at the protein level. We further show that hb mRNA and protein expression in both anterior and posterior domains significantly differs at specific time points. The formation of hb PS4 stripe at the mRNA level proceeds five times faster than at the level of protein. With regard to spatial expression, we show that the offset between posterior gt mRNA and protein domains is much larger in Kr mutants than in wild type embryos and heterozygotes. Finally, we analyze differences in positional variability of eve stripe 7 expression in Kr mutants and Kr+/Kr- heterozygotes at the level of mRNA and protein. These results enable further perspectives to uncover mechanisms underlying discrepancies between mRNA and protein expression in early embryo.

Dynamic maternal gradients and morphogenetic networks in Drosophila early embryo.

Here, we review the latest publications on dynamics and interpretation of morphogen gradients in Drosophila early embryo. The instructive cues provided by these gradients are further interpreted by target genes to enable correct cell fate specification. Moreover, recent studies point on the dynamic and active input from gradients themselves. Latest research has demonstrated that the decay of maternal gradients affects the positional dynamics of target gap gene expression. Study of temporal interpretation of the Bicoid (Bcd) morphogen signal revealed that the Bcd-dependent cell fate specification proceeds sequentially from the posterior towards the anterior. The 'morphogenetic network' of maternal and zygotic regulators functions to spatially organize expression of Bcd-dependent target genes. The analysis of molecular mechanisms of Bcd interaction with target enhancers showed that Bcd gradient induces different chromatin states depending on its concentration. All these findings significantly update the model of morphogen gradient interpretation and provide new perspectives to uncover mechanisms of pattern formation in early embryogenesis.

Preparation of Drosophila embryos for quantitative imaging of gene expression.

Quantitative measurements derived using sophisticated microscopy techniques are essential for understanding the basic principles that control the behavior of biological systems. We have developed a five-step data pipeline to extract quantitative data on segmentation gene expression from confocal images of gene expression patterns in Drosophila. This protocol describes the preparation of Drosophila embryos for imaging by confocal microscopy. Embryos are collected at the appropriate developmental stage and fixed. They are then stained with both primary antibodies and secondary antibodies conjugated with fluorophores to reveal the segmentation gene expression patterns.

Quantitative imaging of gene expression in Drosophila embryos.

Quantitative measurements derived using sophisticated microscopy techniques are essential for understanding the basic principles that control the behavior of biological systems. Here we describe a data pipeline developed to extract quantitative data on segmentation gene expression from confocal images of gene expression patterns in Drosophila. The pipeline consists of image segmentation, background removal, temporal characterization of an embryo, data registration, and data averaging. This pipeline has been successfully applied to obtain quantitative gene expression data at cellular resolution in space and at 6.5-min resolution in time. It has also enabled the construction of a spatiotemporal atlas of segmentation gene expression. We describe the software used to construct a workflow for extracting quantitative data on segmentation gene expression and the BREReA package, which implements the methods for background removal and registration of segmentation gene expression patterns.

Lack of tailless leads to an increase in expression variability in Drosophila embryos.

Developmental processes are robust, or canalised: dynamic patterns of gene expression across space and time are regulated reliably and precisely in the presence of genetic and environmental perturbations. It remains unclear whether canalisation relies on specific regulatory factors (such as heat-shock proteins), or whether it is based on more general redundancy and distributed robustness at the network level. The latter explanation implies that mutations in many regulatory factors should exhibit loss of canalisation. Here, we present a quantitative characterisation of segmentation gene expression patterns in mutants of the terminal gap gene tailless (tll) in Drosophila melanogaster. Our analysis provides new insights into the dynamic mechanisms underlying gap gene regulation, and reveals significantly increased variability of gene expression in the mutant compared to the wild-type background. We show that both position and timing of posterior segmentation gene expression domains vary strongly from embryo-to-embryo in tll mutants. This variability must be caused by a vulnerability in the regulatory system which is hidden or buffered in the wild-type, but becomes uncovered by the deletion of tll. Our analysis provides evidence that loss of canalisation in mutants could be more widespread than previously thought.

Quantitative dynamics and increased variability of segmentation gene expression in the Drosophila Krüppel and knirps mutants.

Here we characterize the response of the Drosophila segmentation system to mutations in two gap genes, Kr and kni, in the form of single or double homozygotes and single heterozygotes. Segmentation gene expression in these genotypes was quantitatively monitored with cellular resolution in space and 6.5 to 13min resolution in time. As is the case with wild type, we found that gene expression domains in the posterior portion of the embryo shift to the anterior over time. In certain cases, such as the gt posterior domain in Kr mutants, the shifts are significantly larger than is seen in wild type embryos. We also investigated the effects of Kr and kni on the variability of gene expression. Mutations often produce variable phenotypes, and it is well known that the cuticular phenotype of Kr mutants is variable. We sought to understand the molecular basis of this effect. We find that throughout cycle 14A the relative levels of eve and ftz expression in stripes 2 and 3 are variable among individual embryos. Moreover, in Kr and kni mutants, unlike wild type, the variability in positioning of the posterior Hb domain and eve stripe 7 is not decreased or filtered with time. The posterior Gt domain in Kr mutants is highly variable at early times, but this variability decreases when this domain shifts in the anterior direction to the position of the neighboring Kni domain. In contrast to these findings, positional variability throughout the embryo does not decrease over time in double Kr;kni mutants. In heterozygotes the early expression patterns of segmentation genes resemble patterns seen in homozygous mutants but by the onset of gastrulation they become similar to the wild type patterns. Finally, we note that gene expression levels are reduced in Kr and kni mutant embryos and have a tendency to decrease over time. This is a surprising result in view of the role that mutual repression is thought to play in the gap gene system.

Modeling of gap gene expression in Drosophila Kruppel mutants.

The segmentation gene network in Drosophila embryo solves the fundamental problem of embryonic patterning: how to establish a periodic pattern of gene expression, which determines both the positions and the identities of body segments. The gap gene network constitutes the first zygotic regulatory tier in this process. Here we have applied the systems-level approach to investigate the regulatory effect of gap gene Kruppel (Kr) on segmentation gene expression. We acquired a large dataset on the expression of gap genes in Kr null mutants and demonstrated that the expression levels of these genes are significantly reduced in the second half of cycle 14A. To explain this novel biological result we applied the gene circuit method which extracts regulatory information from spatial gene expression data. Previous attempts to use this formalism to correctly and quantitatively reproduce gap gene expression in mutants for a trunk gap gene failed, therefore here we constructed a revised model and showed that it correctly reproduces the expression patterns of gap genes in Kr null mutants. We found that the remarkable alteration of gap gene expression patterns in Kr mutants can be explained by the dynamic decrease of activating effect of Cad on a target gene and exclusion of Kr gene from the complex network of gap gene interactions, that makes it possible for other interactions, in particular, between hb and gt, to come into effect. The successful modeling of the quantitative aspects of gap gene expression in mutant for the trunk gap gene Kr is a significant achievement of this work. This result also clearly indicates that the oversimplified representation of transcriptional regulation in the previous models is one of the reasons for unsuccessful attempts of mutant simulations.

Mechanisms of developmental robustness.

We present a review of noise buffering mechanisms responsible for developmental robustness. We focus on functions of chaperone Hsp90, miRNA, and cross-regulation of gap genes in Drosophila. The noise buffering mechanisms associated with these functions represent specific examples of the developmental canalization, reducing the phenotypical variability in presence of either genetic or environmental perturbations. We demonstrate that robustness often appears as a function of a network of interacting elements and that the system level approach is needed to understand the mechanisms of noise filtering.

A regression system for estimation of errors introduced by confocal imaging into gene expression data in situ.

BACKGROUND: Accuracy of the data extracted from two-dimensional confocal images is limited due to experimental errors that arise in course of confocal scanning. The common way to reduce the noise in images is sequential scanning of the same specimen several times with the subsequent averaging of multiple frames. Attempts to increase the dynamical range of an image by setting too high values of microscope PMT parameters may cause clipping of single frames and introduce errors into the data extracted from the averaged images. For the estimation and correction of this kind of errors a method based on censoring technique (Myasnikova et al., 2009) is used. However, the method requires the availability of all the confocal scans along with the averaged image, which is normally not provided by the standard scanning procedure. RESULTS: To predict error size in the data extracted from the averaged image we developed a regression system. The system is trained on the learning sample composed of images obtained from three different microscopes at different combinations of PMT parameters, and for each image all the scans are saved. The system demonstrates high prediction accuracy and was applied for correction of errors in the data on segmentation gene expression in Drosophila blastoderm stored in the FlyEx database (http://urchin.spbcas.ru/flyex/, http://flyex.uchicago.edu/flyex/). The prediction method is realized as a software tool CorrectPattern freely available at http://urchin.spbcas.ru/asp/2011/emm/. CONCLUSIONS: We created a regression system and software to predict the magnitude of errors in the data obtained from a confocal image based on information about microscope parameters used for the image acquisition. An important advantage of the developed prediction system is the possibility to accurately correct the errors in data obtained from strongly clipped images, thereby allowing to obtain images of the higher dynamical range and thus to extract more detailed quantitative information from them.

Canalization of gene expression in the Drosophila blastoderm by gap gene cross regulation.

Developing embryos exhibit a robust capability to reduce phenotypic variations that occur naturally or as a result of experimental manipulation. This reduction in variation occurs by an epigenetic mechanism called canalization, a phenomenon which has resisted understanding because of a lack of necessary molecular data and of appropriate gene regulation models. In recent years, quantitative gene expression data have become available for the segment determination process in the Drosophila blastoderm, revealing a specific instance of canalization. These data show that the variation of the zygotic segmentation gene expression patterns is markedly reduced compared to earlier levels by the time gastrulation begins, and this variation is significantly lower than the variation of the maternal protein gradient Bicoid. We used a predictive dynamical model of gene regulation to study the effect of Bicoid variation on the downstream gap genes. The model correctly predicts the reduced variation of the gap gene expression patterns and allows the characterization of the canalizing mechanism. We show that the canalization is the result of specific regulatory interactions among the zygotic gap genes. We demonstrate the validity of this explanation by showing that variation is increased in embryos mutant for two gap genes, Krüppel and knirps, disproving competing proposals that canalization is due to an undiscovered morphogen, or that it does not take place at all. In an accompanying article in PLoS Computational Biology (doi:10.1371/journal.pcbi.1000303), we show that cross regulation between the gap genes causes their expression to approach dynamical attractors, reducing initial variation and providing a robust output. These results demonstrate that the Bicoid gradient is not sufficient to produce gap gene borders having the low variance observed, and instead this low variance is generated by gap gene cross regulation. More generally, we show that the complex multigenic phenomenon of canalization can be understood at a quantitative and predictive level by the application of a precise dynamical model.

GCPReg package for registration of the segmentation gene expression data in Drosophila.

In modern functional genomics registration techniques areused to construct reference gene expression patterns and createa spatiotemporal atlas of the expression of all the genes in anetwork. In this paper we present a software package calledGCPReg, which can be used to register the expression patterns ofsegmentation genes in the early Drosophila embryo. The key task,which this package performs, is the extraction of spatially localizedcharacteristic features of expression patterns. To facilitatethis task, we have developed an easy-to-use interactive graphicalinterface. We describe GCPReg usage and demonstrate how thispackage can be applied to register gene expression patterns inwild type and mutants. GCPReg has been designed to operate ona UNIX platform and is freely available via the Internet at http://urchin.spbcas.ru/downloads/GCPReg/GCPReg.htm.

Canalization of gene expression and domain shifts in the Drosophila blastoderm by dynamical attractors.

The variation in the expression patterns of the gap genes in the blastoderm of the fruit fly Drosophila melanogaster reduces over time as a result of cross regulation between these genes, a fact that we have demonstrated in an accompanying article in PLoS Biology (see Manu et al., doi:10.1371/journal.pbio.1000049). This biologically essential process is an example of the phenomenon known as canalization. It has been suggested that the developmental trajectory of a wild-type organism is inherently stable, and that canalization is a manifestation of this property. Although the role of gap genes in the canalization process was established by correctly predicting the response of the system to particular perturbations, the stability of the developmental trajectory remains to be investigated. For many years, it has been speculated that stability against perturbations during development can be described by dynamical systems having attracting sets that drive reductions of volume in phase space. In this paper, we show that both the reduction in variability of gap gene expression as well as shifts in the position of posterior gap gene domains are the result of the actions of attractors in the gap gene dynamical system. Two biologically distinct dynamical regions exist in the early embryo, separated by a bifurcation at 53% egg length. In the anterior region, reduction in variation occurs because of stability induced by point attractors, while in the posterior, the stability of the developmental trajectory arises from a one-dimensional attracting manifold. This manifold also controls a previously characterized anterior shift of posterior region gap domains. Our analysis shows that the complex phenomena of canalization and pattern formation in the Drosophila blastoderm can be understood in terms of the qualitative features of the dynamical system. The result confirms the idea that attractors are important for developmental stability and shows a richer variety of dynamical attractors in developmental systems than has been previously recognized.

Estimation of errors introduced by confocal imaging into the data on segmentation gene expression in Drosophila.

MOTIVATION: Currently the confocal scanning microscopy of fluorescently tagged molecules is extensively employed to acquire quantitative data on gene expression at cellular resolution. Following this approach, we generated a large dataset on the expression of segmentation genes in the Drosophila blastoderm, that is widely used in systems biology studies. As data accuracy is of critical importance for the success of studies in this field, we took a shot to evaluate possible errors introduced in the data by acquisition and processing methods. This article deals with errors introduced by confocal microscope. RESULTS: In confocal imaging, the inevitable photon noise is commonly reduced by the averaging of multiple frames. The averaging may introduce errors into the data, if single frames are clipped by microscope hardware. A method based on censoring technique is used to estimate and correct this type of errors. Additional source of errors is the quantification of blurred images. To estimate and correct these errors, the Richardson-Lucy deconvolution method was modified to provide the higher accuracy of data read off from blurred images of the Drosophila blastoderm. We have found that the sizes of errors introduced by confocal imaging make up approximately 5-7% of the mean intensity values and do not disguise the dynamic behavior and characteristic features of gene expression patterns. We also defined a range of microscope parameters for the acquisition of sufficiently accurate data. AVAILABILITY: http://urchin.spbcas.ru/downloads/step/step.htm