Spatiotemporal Gene Expression Patterns Research Articles

Controlling spatiotemporal pattern formation in a concentration gradient with a synthetic toggle switch.

The formation of spatiotemporal patterns of gene expression is frequently guided by gradients of diffusible signaling molecules. The toggle switch subnetwork, composed of two cross‐repressing transcription factors, is a common component of gene regulatory networks in charge of patterning, converting the continuous information provided by the gradient into discrete abutting stripes of gene expression. We present a synthetic biology framework to understand and characterize the spatiotemporal patterning properties of the toggle switch. To this end, we built a synthetic toggle switch controllable by diffusible molecules in Escherichia coli. We analyzed the patterning capabilities of the circuit by combining quantitative measurements with a mathematical reconstruction of the underlying dynamical system. The toggle switch can produce robust patterns with sharp boundaries, governed by bistability and hysteresis. We further demonstrate how the hysteresis, position, timing, and precision of the boundary can be controlled, highlighting the dynamical flexibility of the circuit.

Populations of genetic circuits are unable to find the fittest solution in a multilevel genotype-phenotype map.

The evolution of gene regulatory networks (GRNs) is of great relevance for both evolutionary and synthetic biology. Understanding the relationship between GRN structure and its function can allow us to understand the selective pressures that have shaped a given circuit. This is especially relevant when considering spatio-temporal expression patterns, where GRN models have been shown to be extremely robust and evolvable. However, previous models that studied GRN evolution did not include the evolution of protein and genetic elements that underlie GRN architecture. Here we use toyLIFE, a multilevel genotype-phenotype map, to show that not all GRNs are equally likely in genotype space and that evolution is biased to find the most common GRNs. toyLIFE rules create Boolean GRNs that, embedded in a one-dimensional tissue, develop a variety of spatio-temporal gene expression patterns. Populations of toyLIFE organisms choose the most common GRN out of a set of equally fit alternatives and, most importantly, fail to find a target pattern when it is very rare in genotype space. Indeed, we show that the probability of finding the fittest phenotype increases dramatically with its abundance in genotype space. This phenotypic bias represents a mechanism that can prevent the fixation in the population of the fittest phenotype, one that is inherent to the structure of genotype space and the genotype-phenotype map.

Transcriptome of Pectobacterium carotovorum subsp. carotovorum PccS1 infected in calla plants in vivo highlights a spatiotemporal expression pattern of genes related to virulence, adaptation, and host response.

Bacterial pathogens from the genus Pectobacterium cause soft rot in various plants, and result in important economic losses worldwide. We understand much about how these pathogens digest their hosts and protect themselves against plant defences, as well as some regulatory networks in these processes. However, the spatiotemporal expression of genome‐wide infection of Pectobacterium remains unclear, although researchers analysed this in some phytopathogens. In the present work, comparing the transcriptome profiles from cellular infection with growth in minimal and rich media, RNA‐Seq analyses revealed that the differentially expressed genes (log2‐fold ratio ≥ 1.0) in the cells of Pectobacterium carotovorum subsp. carotovorum PccS1 recovered at a series of time points after inoculation in the host in vivo covered approximately 50% of genes in the genome. Based on the dynamic expression changes in infection, the significantly differentially expressed genes (log2‐fold ratio ≥ 2.0) were classified into five types, and the main expression pattern of the genes for carbohydrate metabolism underlying the processes of infection was identified. The results are helpful to our understanding of the inducement of host plant and environmental adaption of Pectobacterium. In addition, our results demonstrate that maceration caused by PccS1 is due to the depression of callose deposition in the plant for resistance by the pathogenesis‐related genes and the superlytic ability of pectinolytic enzymes produced in PccS1, rather than the promotion of plant cell death elicited by the T3SS of bacteria as described in previous work.

A GDSL Lipase Is Required for Anther and Pollen Development.

Male fertility in angiosperms requires specific spatiotemporal gene expression patterns and cellular differentiation during male reproductive development ([Wilson and Zhang, 2009][1]; [Lin et al., 2017][2]). The formation of the male gametophyte begins inside the anther with the formation of a large

Image-Based Structural Modeling of the Early-Stage Zebrafish Embryo Brain

During embryogenesis, spatiotemporal gene expression patterns direct the segmentation of primary brain compartments that are progressively refined into functional subunits of the adult brain. At the midbrain-hindbrain boundary (MHB), early errors in gene expression result in structural malformations that impair, by an unknown mechanism, the subsequent functional development of the tectum, cerebellum, and tegmentum which help to coordinate visuomotor responses. A quantitative model of the dynamic relationship between gene expression and structure formation would shed further insight on the role of different genes in this process. However, time-series analysis of the 3-dimensional multiphoton microscopy images of early-stage zebrafish embryo brain primordia poses a distinct challenge since the brain initially lacks recognizable landmarks to use as reference. Additional difficulty arises due to noise and pixel intensity variations within and across image data sets. To address these difficulties, we have developed a computational method that assigns simple geometric elements such as “beads” and “bonds” in a hierarchically refined structural model of the zebrafish brain for comparing different phenotypes and developmental time points. Our measurement highlight morphological differences between the wild type and zebrafish fgf8a(acerebellar) mutant that does not form MHB. Our approach is applicable to building quantitative surface models for images of other systems.

Enhancers Facilitate the Birth of De Novo Genes and Gene Integration into Regulatory Networks.

Regulatory networks control the spatiotemporal gene expression patterns that give rise to and define the individual cell types of multicellular organisms. In eumetazoa, distal regulatory elements called enhancers play a key role in determining the structure of such networks, particularly the wiring diagram of “who regulates whom.” Mutations that affect enhancer activity can therefore rewire regulatory networks, potentially causing adaptive changes in gene expression. Here, we use whole-tissue and single-cell transcriptomic and chromatin accessibility data from mouse to show that enhancers play an additional role in the evolution of regulatory networks: They facilitate network growth by creating transcriptionally active regions of open chromatin that are conducive to de novo gene evolution. Specifically, our comparative transcriptomic analysis with three other mammalian species shows that young, mouse-specific intergenic open reading frames are preferentially located near enhancers, whereas older open reading frames are not. Mouse-specific intergenic open reading frames that are proximal to enhancers are more highly and stably transcribed than those that are not proximal to enhancers or promoters, and they are transcribed in a limited diversity of cellular contexts. Furthermore, we report several instances of mouse-specific intergenic open reading frames proximal to promoters showing evidence of being repurposed enhancers. We also show that open reading frames gradually acquire interactions with enhancers over macroevolutionary timescales, helping integrate genes—those that have arisen de novo or by other means—into existing regulatory networks. Taken together, our results highlight a dual role of enhancers in expanding and rewiring gene regulatory networks.

Evaluation of BMP-mediated patterning in a 3D mathematical model of the zebrafish blastula embryo.

Bone Morphogenetic Proteins (BMPs) play an important role in dorsal-ventral (DV) patterning of the early zebrafish embryo. BMP signaling is regulated by a network of extracellular and intracellular factors that impact the range and signaling of BMP ligands. Recent advances in understanding the mechanism of pattern formation support a source-sink mechanism, however it is not clear how the source-sink mechanism shapes patterns in 3D, nor how sensitive the pattern is to biophysical rates and boundary conditions along both the anteroposterior (AP) and DV axes of the embryo. We propose a new three-dimensional growing Partial Differential Equation (PDE)-based model to simulate the BMP patterning process during the blastula stage. This model provides a starting point to elucidate how different mechanisms and components work together in 3D to create and maintain the BMP gradient in the embryo. We also show how the 3D model fits the BMP signaling gradient data at multiple time points along both axes. Furthermore, sensitivity analysis of the model suggests that the spatiotemporal patterns of Chordin and BMP ligand gene expression are dominant drivers of shape in 3D and more work is needed to quantify the spatiotemporal profiles of gene and protein expression to further refine the models.

Spatiotemporal gene expression patterns reveal molecular relatedness between retinal laminae.

In several areas of the central nervous system, neurons are regionally organized into groups or layers that carry out specific activities. In this form of patterning, neurons of distinct types localize their cell bodies to just one or a few of the layers within a structure. However, little is known about whether diverse neuron types within a lamina share molecular features that coordinate their organization. To begin to identify such candidates, we used the laminated murine retina to screen 92 lacZ reporter lines available through the Knockout Mouse Project. Thirty-two of these displayed reporter expression in restricted subsets of inner retina neurons. We then identified the spatiotemporal expression patterns of these genes at key developmental stages. This uncovered several that were heavily enriched in development but reduced in adulthood, including the transcriptional regulator Hmga1. An additional set of genes displayed maturation associated laminar enrichment. Among these, we identified Bbox1 as a novel gene that specifically labels all neurons in the ganglion cell layer but is largely excluded from otherwise molecularly similar neurons in the inner retina. Finally, we established Dbn1 as a new marker enriched in amacrines and Fmnl3 as a marker for subsets of αRGCs. Together, these data provide a spatiotemporal map for laminae-specific molecules and suggest that diverse neuron types within a lamina share coordinating molecular features that may inform their fate or function.

Enhancer Dysfunction in 3D Genome and Disease.

Spatiotemporal patterns of gene expression depend on enhancer elements and other factors during individual development and disease progression. The rapid progress of high-throughput techniques has led to well-defined enhancer chromatin properties. Various genome-wide methods have revealed a large number of enhancers and the discovery of three-dimensional (3D) genome architecture showing the distant interacting mechanisms of enhancers that loop to target gene promoters. Whole genome sequencing projects directed at cancer have led to the discovery of substantial enhancer dysfunction in misregulating gene expression and in tumor initiation and progression. Results from genome-wide association studies (GWAS) combined with functional genomics analyses have elucidated the functional impacts of many cancer risk-associated variants that are enriched within the enhancer regions of chromatin. Risk variants dysregulate the expression of enhancer variant-associated genes via 3D genomic interactions. Moreover, these enhancer variants often alter the chromatin binding affinity for cancer-relevant transcription factors, which in turn leads to aberrant expression of the genes associated with cancer susceptibility. In this review, we investigate the extent to which these genetic regulatory circuits affect cancer predisposition and how the recent development of genome-editing methods have enabled the determination of the impacts of genomic variation and alteration on cancer phenotype, which will eventually lead to better management plans and treatment responses to human cancer in the clinic.

Transcriptional characterisation of human lung cells identifies novel mesenchymal lineage markers.

The lung mesenchyme gives rise to multiple distinct lineages of cells in the mature respiratory system, including smooth muscle cells of the airway and vasculature. However, a thorough understanding of the specification and mesenchymal cell diversity in the human lung is lacking. We completed single-cell RNA sequencing analysis of fetal human lung tissues. Canonical correlation analysis, clustering, cluster marker gene identification and t-distributed stochastic neighbour embedding representation was performed in Seurat. Cell populations were annotated using ToppFun. Immunohistochemistry and in situ hybridisation were used to validate spatiotemporal gene expression patterns for key marker genes. We identified molecularly distinct populations representing "committed" fetal human lung endothelial cells, pericytes and smooth muscle cells. Early endothelial lineages expressed "classic" endothelial cell markers (platelet endothelial cell adhesion molecule/CD31 and claudin 5), while pericytes expressed platelet-derived growth factor receptor-β, Thy-1 membrane glycoprotein and basement membrane molecules (collagen IV, laminin and proteoglycans). We observed a large population of "nonspecific" human lung mesenchymal progenitor cells characterised by expression of collagen I and multiple elastin fibre genes (ELN, MFAP2 and FBN1). We closely characterised the diversity of mesenchymal lineages defined by α2-smooth muscle actin (ACTA2) expression. Two cell populations, with the highest levels of ACTA2 transcriptional activity, expressed unique sets of markers associated with airway or vascular smooth muscle cells. Spatiotemporal analysis of these marker genes confirmed early and persistent spatial specification of airway (HHIP, MYLK and IGF1) and vascular (NTRK3 and MEF2C) smooth muscle cells in the developing human lung. Our data suggest that specification of distinct airway and vascular smooth muscle cell phenotypes is established early in development and can be identified using the markers we provide.

Alternative splicing, spatiotemporal expression of TEP family genes in Yesso scallop (Patinopecten yessoensis) and their disparity in responses to ocean acidification

The complement system constitutes a highly sophisticated and powerful body defense machinery acting in the innate immunity of both vertebrates and invertebrates. As central components of the complement system, significant effects of thioester-containing protein (TEP) family members on immunity have been reported in most vertebrates and in some invertebrates, but the spatiotemporal expression and regulatory patterns of TEP family genes under environmental stress have been less widely investigated in scallops. In this study, expression profiling of TEP family members in the Yesso scallop Patinopecten yessoensis (designated PyTEPs) was performed at all developmental stages, in different healthy adult tissues, and in mantles during exposure to different levels of acidification (pH = 6.5 and 7.5) for different time points (3, 6, 12 and 24 h); this profiling was accomplished through in silico analysis of transcriptome and genome databases. Spatiotemporal expression patterns revealed that PyTEPs had specific functional differentiation in all stages of growth and development of the scallop. Expression analysis confirmed the inducible expression patterns of PyTEPs during exposure to acidification. Gene duplication and alternative splicing events simultaneously occurred in PyTEP1. Seven different cDNA variants of PyTEP1 (designated PyTEP1-A–PyTEP1-G) were identified in the scallop mantle transcriptome during acidic stress. These variants were produced by the alternative splicing of seven differentially transcribed exons (exons 18–24), which encode the highly variable central region. The responses to immune stress may have arisen through the gene duplication and alternative splicing of PyTEP1. The sequence diversity of PyTEP1 isoforms and their different expression profiles in response to ocean acidification (OA) suggested a mechanism used by scallops to differentiate and regulate PyTEP1 gene expression. Collectively, these results demonstrate the gene duplication and alternative splicing of TEP family genes and provide valuable resources for elucidating their versatile roles in bivalve innate immune responses to OA challenge.

Immunostaining of the Embryonic and Larval Drosophila Brain.

Immunostaining is used to visualize the spatiotemporal expression pattern of developmental control genes that regulate the genesis and specification of the embryonic and larval brain of Drosophila. It is also used to visualize the effects of targeted misexpression or inactivation of disease-related genes. Immunostaining uses specific antibodies to mark expressed proteins and allows their localization to be traced. This method reveals insights into gene regulation, cell type specification, neuron and glial differentiation, axonal and synaptic scaffolding and posttranslational protein modifications underlying the patterning and specification of the maturing brain. Depending on the targeted protein, it is possible to visualize a multitude of regions of the Drosophila brain, such as small groups of neurons or glia, defined subcomponents of the brain's axon scaffold, or pre- and postsynaptic structures of neurons. Thus, antibody probes that recognize defined tissues, cells, or subcellular structures like axons or synaptic terminals can be used as markers to identify and analyze phenotypes in embryos and larvae. Several antibodies, combined with different labels can be used concurrently to examine protein colocalization. This protocol spans over 3-4days.

Spectral recording of gene expression history by fluorescent timer protein.

Monitoring spatio-temporal patterns of gene expression by fluorescent proteins requires longitudinal observation, which is often difficult to implement. Here, we fuse a fluorescent timer (FT) protein with an immediate early gene (IEG) promoter to track live gene expression in single cells. This results in a stimulus- and time-dependent spectral shift from blue to red for subsequent monitoring with fluorescence activated cell sorting (FACS) and live cell imaging. This spectral shift enables imputing the time point of activity post-hoc to dissociate early and late responders from a single snapshot in time. Thus, we provide a tool for tracking stimulus-driven IEG expression and demonstrate proof of concept exploiting promoter::FT fusions, adding new dimensions to experiments that require reconstructing spatio-temporal patterns of gene expression in cells, tissues or living organisms.

Spatiotemporal expression pattern of alpha, beta, and gamma genes during BoHV-5 infection

Bovine herpesvirus 5 is an alphaherpesvirus that causes nonsuppurative meningoencephalitis in cattle. This disease occurs naturally in either outbreaks or isolated cases, and exhibits low morbidity and high lethality. Although previous studies elucidated crucial aspects involved in the pathogenesis of the disease, there is a paucity of information regarding the molecular events contributing to infection and replication of BoHV-5. The objective of the present study was to determine the in vitro gene expression pattern of BoHV-5 (e.g., alpha, beta, and gamma genes) and host cells genes (GAPDH and 18S) over time utilizing different quantities of inoculated virus. Three BoHV-5 genes (bICP0, UL9, US4) and one structural bovine cell gene had their expression accessed by real-time PCR. While the expression of BoHV-5 genes increased during the course of infection, GAPDH gene expression decreased in the host cells, evidencing the effect of viral infection on the expression of bovine cell genes. The 18S ribosomal RNA (rRNA) gene was constitutively expressed throughout BoHV-5 infection. Our data clearly demonstrates that GAPDH gene should not be used as a reference gene in studies of BoHV-5 infection because it was influenced by viral infection. However, 18S rRNA was constitutively expressed and, therefore, is recommended for normalization of BoHV-5 infection studies in bovine cells. The expression of viral genes transcripts was not altered by increasing number of viral particles added to the culture. All viral genes included here demonstrated the same expression pattern over time and there was no difference in the expression of viral genes among the various time points. Our data show important differences comparing to classical studies regarding herpesvirus alpha, beta, and gamma genes expression. More research is necessary to improve our understanding about the BoHV-5 biology during infection. Studies employing next-generation sequencing (i.e., RNA-seq), using both in vitro and in vivo models, would be the next logical step to grasp the virus and host cell’s transcriptome changes over the course of infection.

Quantitative Comparison of the Anterior-Posterior Patterning System in the Embryos of Five Drosophila Species.

Complex spatiotemporal gene expression patterns direct the development of the fertilized egg into an adult animal. Comparisons across species show that, in spite of changes in the underlying regulatory DNA sequence, developmental programs can be maintained across millions of years of evolution. Reciprocally, changes in gene expression can be used to generate morphological novelty. Distinguishing between changes in regulatory DNA that lead to changes in gene expression and those that do not is therefore a central goal of evolutionary developmental biology. Quantitative, spatially-resolved measurements of developmental gene expression patterns play a crucial role in this goal, enabling the detection of subtle phenotypic differences between species and the development of computations models that link the sequence of regulatory DNA to expression patterns. Here we report the generation of two atlases of cellular resolution gene expression measurements for the primary anterior-posterior patterning genes in Drosophila simulans and Drosophila virilis. By combining these data sets with existing atlases for three other Drosophila species, we detect subtle differences in the gene expression patterns and dynamics driving the highly conserved axis patterning system and delineate inter-species differences in the embryonic morphology. These data sets will be a resource for future modeling studies of the evolution of developmental gene regulatory networks.

Genome-wide analysis of spatiotemporal gene expression patterns during floral organ development in Brassica rapa.

Flowering is a key agronomic trait that directly influences crop yield and quality and serves as a model system for elucidating the molecular basis that controls successful reproduction, adaptation, and diversification of flowering plants. Adequate knowledge of continuous series of expression data from the floral transition to maturation is lacking in Brassica rapa. To unravel the genome expression associated with the development of early small floral buds (< 2mm; FB2), early large floral buds (2-4mm; FB4), stamens (STs) and carpels (CPs), transcriptome profiling was carried out with a Br300K oligo microarray. The results showed that at least 6848 known nonredundant genes (30% of the genes of the Br300K) were differentially expressed during the floral transition from vegetative tissues to maturation. Functional annotation of the differentially expressed genes (DEGs) (fold change ≥ 5) by comparison with a close relative, Arabidopsis thaliana, revealed 6552 unigenes (4579 upregulated; 1973 downregulated), including 131 Brassica-specific and 116 functionally known floral Arabidopsis homologs. Additionally, 1723, 236 and 232 DEGs were preferentially expressed in the tissues of STs, FB2, and CPs. These DEGs also included 43 transcription factors, mainly AP2/ERF-ERF, NAC, MADS-MIKC, C2H2, bHLH, and WRKY members. The differential gene expression during flower development induced dramatic changes in activities related to metabolic processes (23.7%), cellular (22.7%) processes, responses to the stimuli (7.5%) and reproduction (1%). A relatively large number of DEGs were observed in STs and were overrepresented by photosynthesis-related activities. Subsequent analysis via semiquantitative RT-PCR, histological analysis performed with in situ hybridization of BrLTP1 and transgenic reporter lines (BrLTP promoter::GUS) of B. rapa ssp. pekinensis supported the spatiotemporal expression patterns. Together, these results suggest that a temporally and spatially regulated process of the selective expression of distinct fractions of the same genome leads to the development of floral organs. Interestingly, most of the differentially expressed floral transcripts were located on chromosomes 3 and 9. This study generated a genome expression atlas of the early floral transition to maturation that represented the flowering regulatory elements of Brassica rapa.

Fucosyltransferase gene expression in goat endometrium during the estrous cycle and early pregnancy

Regulation of the expression of the alpha(1,2)fucosyltransferase (FUT) genes and their enzymatic products, including the H-type 1 antigen (HT1), on the luminal surface of the uterus is believed to be critical for establishment of pregnancy in mammals. The FUT1 gene is a marker for conception rates in dairy cows and HT1 is a marker for uterine receptivity in rodents. To determine the spatiotemporal expression patterns of FUT1 and FUT2 genes in goats, endometrial tissues were obtained on six days spanning the estrous cycle (Days 5, 11, 13, 15, 17 and 19) and seven days spanning early pregnancy (Days 5, 11, 13, 15, 17, 19 and 25). In all data, we found no effect of status (cyclic or pregnant; P > 0.1) and pooled data where appropriate. We cloned FUT1 cDNA from goat endometrium and made probes from it for Northern and slot blot analyses. The analyses indicated that FUT1 gene expression was high until Day 13, and then declined. In situ hybridization revealed a change in the cell-specificity of FUT1 gene expression over the estrous cycle and early pregnancy. In situ hybridization signal intensity scores indicated that FUT1 expression by uterine epithelium was high on Day 5, moderate on Day 11, and minimal on subsequent days. In situ hybridization signals in uterine glandular epithelial cells remained high from Day 5 to Day 13, with weaker signals thereafter. Quantitative reverse transcription-PCR (RT-qPCR) assays were used for quantitation of FUT1 and FUT2 mRNAs. Quantitative RT-qPCR data were generated from endometrium collected from cyclic and pregnant animals on Days 5, 11 and 17. Relative levels of FUT1 mRNA were high on Days 5 and 11, but then fell 5-fold by Day 17 (P < 0.01). FUT2 mRNA concentrations were below the accurate detectable limit of the assay. High levels of HT1 were observed on the apical surface of uterine luminal epithelia on Days 5, 15, 17 and 19, with much lower levels on Days 11 and 13. Thus, data suggests that FUT1 is the primary enzyme responsible for the high levels of HT1 antigen present on the uterine luminal epithelium between Days 5 and 11 of the estrous cycle and early pregnancy. But changes in the expression of the FUT1 gene does not directly correlate to HT1 staining, which increased from Day 13–15. Future studies are required to understand the regulation of the HT1 antigen on the luminal surface of endometrium.

Expression analysis of four pseudo-response regulator (PRR) genes in Chrysanthemum morifolium under different photoperiods.

Genes encoding pseudo-response regulator (PRR) proteins play significant roles in plant circadian clocks. In this study, four genes related to flowering time were isolated from Chrysanthemum morifolium. Phylogenetic analysis showed that they are highly homologous to the counterparts of PRRs of Helianthus annuus and named as CmPRR2, CmPRR7, CmPRR37, and CmPRR73. Conserved motifs prediction indicated that most of the closely related members in the phylogenetic tree share common protein sequence motifs, suggesting functional similarities among the PRR proteins within the same subtree. In order to explore functions of the genes, we selected two Chrysanthemum varieties for comparison; that is, a short-day sensitive Zijiao and a short-day insensitive Aoyunbaixue. Compared to Aoyunbaixue, Zijiao needs 13 more days to complete the flower bud differentiation. Evidence from spatio-temporal gene expression patterns demonstrated that the CmPRRs are highly expressed in flower and stem tissues, with a growing trend across the Chrysanthemum developmental process. In addition, we also characterized the CmPRRs expression patterns and found that CmPRRs can maintain their circadian oscillation features to some extent under different photoperiod treatment conditions. These lines of evidence indicated that the four CmPRRs undergo circadian oscillation and possibly play roles in regulating the flowering time of C. morifolium.

Sequence based prediction of enhancer regions from DNA random walk

Regulatory elements play a critical role in development process of eukaryotic organisms by controlling the spatio-temporal pattern of gene expression. Enhancer is one of these elements which contributes to the regulation of gene expression through chromatin loop or eRNA expression. Experimental identification of a novel enhancer is a costly exercise, due to which there is an interest in computational approaches to predict enhancer regions in a genome. Existing computational approaches to achieve this goal have primarily been based on training of high-throughput data such as transcription factor binding sites (TFBS), DNA methylation, and histone modification marks etc. On the other hand, purely sequence based approaches to predict enhancer regions are promising as they are not biased by the complexity or context specificity of such datasets. In sequence based approaches, machine learning models are either directly trained on sequences or sequence features, to classify sequences as enhancers or non-enhancers. In this paper, we derived statistical and nonlinear dynamic features along with k-mer features from experimentally validated sequences taken from Vista Enhancer Browser through random walk model and applied different machine learning based methods to predict whether an input test sequence is enhancer or not. Experimental results demonstrate the success of proposed model based on Ensemble method with area under curve (AUC) 0.86, 0.89, and 0.87 in B cells, T cells, and Natural killer cells for histone marks dataset.

Local Auxin Biosynthesis Is a Key Regulator of Plant Development

Auxin is a major phytohormone that controls numerous aspects of plant development and coordinates plant responses to the environment. Morphogenic gradients of auxin govern cell fate decisions and underlie plant phenotypic plasticity. Polar auxin transport plays a central role in auxin maxima generation. The discovery of the exquisite spatiotemporal expression patterns of auxin biosynthesis genes of the WEI8/TAR and YUC families suggested that local auxin production may contribute to the formation of auxin maxima. Herein, we systematically addressed the role of local auxin biosynthesis in plant development and responses to the stress phytohormone ethylene by manipulating spatiotemporal patterns of WEI8. Our study revealed that local auxin biosynthesis and transport act synergistically and are individually dispensable for root meristem maintenance. In contrast, flower fertility and root responses to ethylene require local auxin production that cannot be fully compensated for by transport in the generation of morphogenic auxin maxima.