Evolution of transcription factor function

Genome-Wide Identification of the Whirly Gene Family and Its Potential Function in Low Phosphate Stress in Soybean (Glycine max)

The Whirly (WHY) gene family, functioning as transcription factors, plays an essential role in the regulation of plant metabolic responses, which has been demonstrated across multiple species. However, the WHY gene family and its functions in soybean remains unclear. In this paper, we conducted genome-wide screening and identification to characterize the WHY gene family. Seven WHY members were identified and randomly distributed across six chromosomes. The phylogenetic evolutionary tree of WHY genes in soybean and other species was divided into five clades. An in-depth analysis revealed that segmental duplications significantly contributed to the expansion of GmWHYs, and the GmWHY gene members may have experienced evolutionary pressure for purifying selection in soybeans. The analysis of promoter Cis-elements in GmWHYs suggested their potential significance in addressing diverse stress conditions. The expression patterns of GmWHYs exhibited tissue-specific variations throughout the different stages of soybean development. Additionally, six GmWHY genes exhibited different responses to low phosphate stress. These findings will provide a theoretical basis and valuable reference for the future exploration of WHY gene function.

Divergence of Grainy head affects chromatin accessibility, gene expression, and embryonic viability in Drosophila melanogaster

Pioneer factors are critical for gene regulation and development because they bind chromatin and make DNA more accessible for binding by other transcription factors. The pioneer factor Grainy head (Grh) is present across metazoans and has been shown to retain a role in epithelium development in fruit flies, nematodes, and mice despite extensive divergence in both amino acid sequence and length. Here, we investigate the evolution of Grh function by comparing the effects of the fly (Drosophila melanogaster) and worm (Caenorhabditis elegans) Grh orthologs on chromatin accessibility, gene expression, embryonic development, and viability in transgenic D. melanogaster. We found that the Caenorhabditis elegans ortholog rescued cuticle development but not full embryonic viability in Drosophila melanogaster grh null mutants. At the molecular level, the C. elegans ortholog only partially rescued chromatin accessibility and gene expression. Divergence in the disordered N-terminus of the Grh protein contributes to these differences in embryonic viability and molecular phenotypes. These data show how pioneer factors can diverge in sequence and function at the molecular level while retaining conserved developmental functions at the organismal level.

Identification, characterization, and prognosis investigation of pivotal genes shared in different stages of breast cancer

One of the leading causes of death (20.1 per 100,000 women per year), breast cancer is the most prevalent cancer in females. Statistically, 95% of breast cancer are categorized as adenocarcinomas, and 55% of all patients may go into invasive phases; however, it can be successfully treated in approximately 70-80% of cases if diagnosed in the nascent stages. The emergence of breast tumor cells which are intensely resistant to conventional therapies, along with the high rate of metastasis occurrence, has highlighted the importance of finding novel strategies and treatments. One of the most advantageous schemes to alleviate this complication is to identify the common differentially expressed genes (DEGs) among primary and metastatic cancerous cells to use resultants for designing new therapeutic agents which are able to target both primary and metastatic breast tumor cells. In this study, the gene expression dataset with accession number GSE55715 was analyzed containing two primary tumor samples, three bone-metastatic samples, and three normal samples to distinguish the up- and down regulated genes in each stage compared to normal cells as control. In the next step, the common upregulated genes between the two experimental groups were detected by Venny online tool. Moreover, gene ontology, functions and pathways, gene-targeting microRNA, and influential metabolites were determined using EnrichR 2021 GO, KEGG pathways miRTarbase 2017, and HMDB 2021, respectively. Furthermore, elicited from STRING protein-protein interaction networks were imported to Cytoscape software to identify the hub genes. Then, identified hub genes were checked to validate the study using oncological databases. The results of the present article disclosed 1263 critical common DEGs (573 upregulated + 690 downregulated), including 35 hub genes that can be broadly used as new targets for cancer treatment and as biomarkers for cancer detection by evaluation of expression level. Besides, this study opens a new horizon to reveal unknown aspects of cancer signaling pathways by providing raw data evoked from in silico experiments. This study's outcomes can also be widely utilized in further lab research since it contains diverse information on common DEGs of varied stages and metastases of breast cancer, their functions, structures, interactions, and associations.

Role of bZIP Transcription Factors in Plant Salt Stress

Soil salinity has become an increasingly serious problem worldwide, greatly limiting crop development and yield, and posing a major challenge to plant breeding. Basic leucine zipper (bZIP) transcription factors are the most widely distributed and conserved transcription factors and are the main regulators controlling various plant response processes against external stimuli. The bZIP protein contains two domains: a highly conserved, DNA-binding alkaline region, and a diverse leucine zipper, which is one of the largest transcription factor families in plants. Plant bZIP is involved in many biological processes, such as flower development, seed maturation, dormancy, and senescence, and plays an important role in abiotic stresses such as salt damage, drought, cold damage, osmotic stress, mechanical damage, and ABA signal response. In addition, bZIP is involved in the regulation of plant response to biological stresses such as insect pests and pathogen infection through salicylic acid, jasmonic acid, and ABA signal transduction pathways. This review summarizes and discusses the structural characteristics and functional characterization of the bZIP transcription factor group, the bZIP transcription factor complex and its molecular regulation mechanisms related to salt stress resistance, and the regulation of transcription factors in plant salt stress resistance. This review provides a theoretical basis and research ideas for further exploration of the salt stress-related functions of bZIP transcription factors. It also provides a theoretical basis for crop genetic improvement and green production in agriculture.

EDomics: a comprehensive and comparative multi-omics database for animal evo-devo

Evolutionary developmental biology (evo-devo) has been among the most fascinating interdisciplinary fields for decades, which aims to elucidate the origin and evolution of diverse developmental processes. The rapid accumulation of omics data provides unprecedented opportunities to answer many interesting but unresolved evo-devo questions. However, the access and utilization of these resources are hindered by challenges particularly in non-model animals. Here, we establish a comparative multi-omics database for animal evo-devo (EDomics, http://edomics.qnlm.ac) containing comprehensive genomes, bulk transcriptomes, and single-cell data across 40 representative species, many of which are generally used as model organisms for animal evo-devo study. EDomics provides a systematic view of genomic/transcriptomic information from various aspects, including genome assembly statistics, gene features and families, transcription factors, transposable elements, and gene expressional profiles/networks. It also exhibits spatiotemporal gene expression profiles at a single-cell level, such as cell atlas, cell markers, and spatial-map information. Moreover, EDomics provides highly valuable, customized datasets/resources for evo-devo research, including gene family expansion/contraction, inferred core gene repertoires, macrosynteny analysis for karyotype evolution, and cell type evolution analysis. EDomics presents a comprehensive and comparative multi-omics platform for animal evo-devo community to decipher the whole history of developmental evolution across the tree of life.

Genome-wide identification, phylogenetic, and expression analysis under abiotic stress conditions of Whirly (WHY) gene family in Medicago sativa L

The WHY family is a group of plant-specific transcription factors, that can bind to single-stranded DNA molecules and play a variety of functions in plant nuclei and organelles, participating in the regulation of plant leaf senescence. It has been identified and analyzed in many species, however, the systematic identification and analysis of the WHY genes family have not yet been reported in alfalfa (Medicago sativa L.). Therefore, to explore the function of alfalfa the WHY genes, and 10 MsWHY genes were identified and further characterized their evolutionary relationship and expression patterns by analyzing the recently published genome of alfalfa. Comprehensive analysis of the chromosome location, physicochemical properties of the protein, evolutionary relationship, conserved motifs, and responses to abiotic stresses of the WHY gene family in alfalfa using bioinformatics methods. The results showed that 10 MsWHY genes were distributed on 10 chromosomes, and collinearity analysis showed that many MsWHYs might be derived from segmental duplications, and these genes are under purifying selection. Based on phylogenetic analyses, the WHY gene family of alfalfa can be divided into four subfamilies: I-IV subfamily, and approximately all the WHY genes within the same subfamily share similar gene structures. The 10 MsWHY gene family members contained 10 motifs, of which motif 2 and motif 4 are the conserved motifs shared by these genes. Furthermore, the analysis of cis-regulatory elements indicated that regulatory elements related to transcription, cell cycle, development, hormone, and stress response are abundant in the promoter sequence of the MsWHY genes. Real-time quantitative PCR demonstrated that MsWHYs gene expression is induced by drought, salt, and methyl jasmonate. The present study serves as a basic foundation for future functional studies on the alfalfa WHY family.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

A combined RAD-Seq and WGS approach reveals the genomic basis of yellow color variation in bumble bee Bombus terrestris

Bumble bees exhibit exceptional diversity in their segmental body coloration largely as a result of mimicry. In this study we sought to discover genes involved in this variation through studying a lab-generated mutant in bumble bee Bombus terrestris, in which the typical black coloration of the pleuron, scutellum, and first metasomal tergite is replaced by yellow, a color variant also found in sister lineages to B. terrestris. Utilizing a combination of RAD-Seq and whole-genome re-sequencing, we localized the color-generating variant to a single SNP in the protein-coding sequence of transcription factor cut. This mutation generates an amino acid change that modifies the conformation of a coiled-coil structure outside DNA-binding domains. We found that all sequenced Hymenoptera, including sister lineages, possess the non-mutant allele, indicating different mechanisms are involved in the same color transition in nature. Cut is important for multiple facets of development, yet this mutation generated no noticeable external phenotypic effects outside of setal characteristics. Reproductive capacity was reduced, however, as queens were less likely to mate and produce female offspring, exhibiting behavior similar to that of workers. Our research implicates a novel developmental player in pigmentation, and potentially caste, thus contributing to a better understanding of the evolution of diversity in both of these processes.

MolluscDB: an integrated functional and evolutionary genomics database for the hyper-diverse animal phylum Mollusca

Mollusca represents the second largest animal phylum but remains poorly explored from a genomic perspective. While the recent increase in genomic resources holds great promise for a deep understanding of molluscan biology and evolution, access and utilization of these resources still pose a challenge. Here, we present the first comprehensive molluscan genomics database, MolluscDB (http://mgbase.qnlm.ac), which compiles and integrates current molluscan genomic/transcriptomic resources and provides convenient tools for multi-level integrative and comparative genomic analyses. MolluscDB enables a systematic view of genomic information from various aspects, such as genome assembly statistics, genome phylogenies, fossil records, gene information, expression profiles, gene families, transcription factors, transposable elements and mitogenome organization information. Moreover, MolluscDB offers valuable customized datasets or resources, such as gene coexpression networks across various developmental stages and adult tissues/organs, core gene repertoires inferred for major molluscan lineages, and macrosynteny analysis for chromosomal evolution. MolluscDB presents an integrative and comprehensive genomics platform that will allow the molluscan community to cope with ever-growing genomic resources and will expedite new scientific discoveries for understanding molluscan biology and evolution.

Gene Coexpression Network Reveals Insights into the Origin and Evolution of a Theanine-Associated Regulatory Module in Non-Camellia and Camellia Species

Theanine (thea) is one of the most important plant-derived characteristic secondary metabolites and a major healthcare product because of its beneficial biological activities, such as being an antianxiety agent, promoting memory, and lowering blood pressure. Thea mostly accumulates in Camellia plants and is especially rich in Camellia sinensis (tea plant). Although some functional genes (e.g., TS, GOGAT, and GS) attributed to thea accumulation have been separately well explored in tea plants, the evolution of a regulatory module (highly interacting gene group) related to thea metabolism remains to be elaborated. Herein, a thea-associated regulatory module (TARM) was mined by using a comprehensive analysis of a weighted gene coexpression network in Camellia and non-Camellia species. Comparative genomic analysis of 84 green plant species revealed that TARM originated from the ancestor of green plants (algae) and that TARM genes were recruited from different evolutionary nodes with the most gene duplication events at the early stage. Among the TARM genes, two core transcription factors of NAC080 and LBD38 were deduced, which may play a crucial role in regulating the biosynthesis of thea. Our findings provide the first insights into the origin and evolution of TARM and indicate a promising paradigm for identifying vital regulatory genes involved in thea metabolism.

Speed-Specificity Trade-Offs in the Transcription Factors Search for Their Genomic Binding Sites

Transcription factors (TFs) regulate gene expression by binding DNA sequences recognized by their DNA-binding domains (DBDs). DBD-recognized motifs are short and highly abundant in genomes. The ability of TFs to bind a specific subset of motif-containing sites, and to do so rapidly upon activation, is fundamental for gene expression in all eukaryotes. Despite extensive interest, our understanding of the TF-target search process is fragmented; although binding specificity and detection speed are two facets of this same process, trade-offs between them are rarely addressed. In this opinion article, we discuss potential speed-specificity trade-offs in the context of existing models. We further discuss the recently described 'distributed specificity' paradigm, suggesting that intrinsically disordered regions (IDRs) promote specificity while reducing the TF-target search time.

Capturing the transcription factor interactome in response to sub-lethal insecticide exposure

The increasing levels of pesticide resistance in agricultural pests and disease vectors represents a threat to both food security and global health. As insecticide resistance intensity strengthens and spreads, the likelihood of a pest encountering a sub-lethal dose of pesticide dramatically increases. Here, we apply dynamic Bayesian networks to a transcriptome time-course generated using sub-lethal pyrethroid exposure on a highly resistant Anopheles coluzzii population. The model accounts for circadian rhythm and ageing effects allowing high confidence identification of transcription factors with key roles in pesticide response. The associations generated by this model show high concordance with lab-based validation and identifies 44 transcription factors putatively regulating insecticide-responsive transcripts. We identify six key regulators, with each displaying differing enrichment terms, demonstrating the complexity of pesticide response. The considerable overlap of resistance mechanisms in agricultural pests and disease vectors strongly suggests that these findings are relevant in a wide variety of pest species.

Mechanistic philosophies of development: Theodor Boveri and Eric H. Davidson

Theodor Boveri's (1862-1915) and Eric Davidson's (1937-2015) achievements represent thoroughly two quite distant time frames in the history of the mechanistic approaches to development, that Jane Maienschein (2014) has characterized respectively as the era of the "experimental embryo" and of the "computed embryo". Nonetheless, Davidson's special bond to Boveri is meant to emphasize the genealogical continuity of an embryological tradition of mechanistic philosophy that, differently from molecular biology, is committed to an explanation of the hereditary transmission of organization. Davidson's genealogical claim is reconsidered through a contextualized analysis of the function of machine-like models and of the role of experiment in the making of their respective mechanistic philosophies. This analysis may help to shed light on resilience and change in the understanding of a mechanistic approach to development.

Divergence in DNA Specificity among Paralogous Transcription Factors Contributes to Their Differential In Vivo Binding

Paralogous transcription factors (TFs) are oftentimes reported to have identical DNA-binding motifs, despite the fact that they perform distinct regulatory functions. Differential genomic targeting by paralogous TFs is generally assumed to be due to interactions with protein co-factors or the chromatin environment. Using a computational-experimental framework called iMADS (integrative modeling and analysis of differential specificity), we show that, contrary to previous assumptions, paralogous TFs bind differently to genomic target sites even in vitro. We used iMADS to quantify, model, and analyze specificity differences between 11 TFs from 4 protein families. We found that paralogous TFs have diverged mainly at medium- and low-affinity sites, which are poorly captured by current motif models. We identify sequence and shape features differentially preferred by paralogous TFs, and we show that the intrinsic differences in specificity among paralogous TFs contribute to their differential in vivo binding. Thus, our study represents a step forward in deciphering the molecular mechanisms of differential specificity in TF families.

mTFkb: a knowledgebase for fundamental annotation of mouse transcription factors

Transcription factors (TFs) are well-known important regulators in cell biology and tissue development. However, in mouse, one of the most widely-used model species, currently the vast majority of the known TFs have not been functionally studied due to the lack of sufficient annotations. To this end, we collected and analyzed the whole transcriptome sequencing data from more than 30 major mouse tissues and used the expression profiles to annotate the TFs. We found that the expression patterns of the TFs are highly correlated with the histology of the tissue types thus can be used to infer the potential functions of the TFs. Furthermore, we found that as many as 30% TFs display tissue-specific expression pattern, and these tissue-specific TFs are among the key TFs in their corresponding tissues. We also observed signals of divergent transcription associated with many TFs with unique expression pattern. Lastly, we have integrated all the data, our analysis results as well as various annotation resources to build a web-based database named mTFkb freely accessible at http://www.myogenesisdb.org/mTFkb/. We believe that mTFkb could serve as a useful and valuable resource for TF studies in mouse.

Composition and genomic organization of arthropod Hox clusters

Background

The ancestral arthropod is believed to have had a clustered arrangement of ten Hox genes. Within arthropods, Hox gene mutations result in transformation of segment identities. Despite the fact that variation in segment number/character was common in the diversification of arthropods, few examples of Hox gene gains/losses have been correlated with morphological evolution. Furthermore, a full appreciation of the variation in the genomic arrangement of Hox genes in extant arthropods has not been recognized, as genome sequences from each major arthropod clade have not been reported until recently. Initial genomic analysis of the chelicerate Tetranychusurticae suggested that loss of Hox genes and Hox gene clustering might be more common than previously assumed. To further characterize the genomic evolution of arthropod Hox genes, we compared the genomic arrangement and general characteristics of Hox genes from representative taxa from each arthropod subphylum.

Results

In agreement with others, we find arthropods generally contain ten Hox genes arranged in a common orientation in the genome, with an increasing number of sampled species missing either Hox3 or abdominal-A orthologs. The genomic clustering of Hox genes in species we surveyed varies significantly, ranging from 0.3 to 13.6 Mb. In all species sampled, arthropod Hox genes are dispersed in the genome relative to the vertebrate Mus musculus. Differences in Hox cluster size arise from variation in the number of intervening genes, intergenic spacing, and the size of introns and UTRs. In the arthropods surveyed, Hox gene duplications are rare and four microRNAs are, in general, conserved in similar genomic positions relative to the Hox genes.

Conclusions

The tightly clustered Hox complexes found in the vertebrates are not evident within arthropods, and differential patterns of Hox gene dispersion are found throughout the arthropods. The comparative genomic data continue to support an ancestral arthropod Hox cluster of ten genes with a shared orientation, with four Hox gene-associated miRNAs, although the degree of dispersion between genes in an ancestral cluster remains uncertain. Hox3 and abdominal-A orthologs have been lost in multiple, independent lineages, and current data support a model in which inversions of the Abdominal-B locus that result in the loss of abdominal-A correlate with reduced trunk segmentation.

Electronic supplementary material

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

Comparative analyses of developmental transcription factor repertoires in sponges reveal unexpected complexity of the earliest animals

Developmental transcription factors (DTFs) control development of animals by affecting expression of target genes, some of which are transcription factors themselves. In bilaterians and cnidarians, conserved DTFs are involved in homologous processes such as gastrulation or specification of neurons. The genome of Amphimedon queenslandica, the first sponge to be sequenced, revealed that only a fraction of these conserved DTF families are present in demosponges. This finding was in line with the view that morphological complexity in the animal lineage correlates with developmental toolkit complexity. However, as the phylum Porifera is very diverse, Amphimedon's genome may not be representative of all sponges. The recently sequenced genomes of calcareous sponges Sycon ciliatum and Leucosolenia complicata allowed investigations of DTFs in a sponge lineage evolutionarily distant from demosponges. Surprisingly, the phylogenetic analyses of identified DTFs revealed striking differences between the calcareous sponges and Amphimedon. As these differences appear to be a result of independent gene loss events in the two sponge lineages, the last common ancestor of sponges had to possess a much more diverse repertoire of DTFs than extant sponges. Developmental expression of sponge homologs of genes involved in specification of the Bilaterian endomesoderm and the neurosensory cells suggests that roles of many DTFs date back to the last common ancestor of all animals. Strikingly, even DTFs displaying apparent pan-metazoan conservation of sequence and function are not immune to being lost from individual species genomes. The quest for a comprehensive picture of the developmental toolkit in the last common metazoan ancestor is thus greatly benefitting from the increasing accessibility of sequencing, allowing comparisons of multiple genomes within each phylum.

How do regulatory networks evolve and expand throughout evolution?

Throughout evolution, regulatory networks need to expand and adapt to accommodate novel genes and gene functions. However, the molecular details explaining how gene networks evolve remain largely unknown. Recent studies demonstrate that changes in transcription factors contribute to the evolution of regulatory networks. In particular, duplication of transcription factors followed by specific mutations in their DNA-binding or interaction domains propels the divergence and emergence of new networks. The innate promiscuity and modularity of regulatory networks contributes to their evolvability: duplicated promiscuous regulators and their target promoters can acquire mutations that lead to gradual increases in specificity, allowing neofunctionalization or subfunctionalization.

Hybrid incompatibility despite pleiotropic constraint in a sequence-based bioenergetic model of transcription factor binding

Hybrid incompatibility can result from gene misregulation produced by divergence in trans-acting regulatory factors and their cis-regulatory targets. However, change in trans-acting factors may be constrained by pleiotropy, which would in turn limit the evolution of incompatibility. We employed a mechanistically explicit bioenergetic model of gene expression wherein parameter combinations (number of transcription factor molecules, energetic properties of binding to the regulatory site, and genomic background size) determine the shape of the genotype-phenotype (G-P) map, and interacting allelic variants of mutable cis and trans sites determine the phenotype along that map. Misregulation occurs when the phenotype differs from its optimal value. We simulated a pleiotropic regulatory pathway involving a positively selected and a conserved trait regulated by a shared transcription factor (TF), with two populations evolving in parallel. Pleiotropic constraints shifted evolution in the positively selected trait to its cis-regulatory locus. We nevertheless found that the TF genotypes often evolved, accompanied by compensatory evolution in the conserved trait, and both traits contributed to hybrid misregulation. Compensatory evolution resulted in "developmental system drift," whereby the regulatory basis of the conserved phenotype changed although the phenotype itself did not. Pleiotropic constraints became stronger and in some cases prohibitive when the bioenergetic properties of the molecular interaction produced a G-P map that was too steep. Likewise, compensatory evolution slowed and hybrid misregulation was not evident when the G-P map was too shallow. A broad pleiotropic "sweet spot" nevertheless existed where evolutionary constraints were moderate to weak, permitting substantial hybrid misregulation in both traits. None of these pleiotropic constraints manifested when the TF contained nonrecombining domains independently regulating the respective traits.

Pervasive divergence of transcriptional gene regulation in Caenorhabditis nematodes

Because there is considerable variation in gene expression even between closely related species, it is clear that gene regulatory mechanisms evolve relatively rapidly. Because primary sequence conservation is an unreliable proxy for functional conservation of cis-regulatory elements, their assessment must be carried out in vivo. We conducted a survey of cis-regulatory conservation between C. elegans and closely related species C. briggsae, C. remanei, C. brenneri, and C. japonica. We tested enhancers of eight genes from these species by introducing them into C. elegans and analyzing the expression patterns they drove. Our results support several notable conclusions. Most exogenous cis elements direct expression in the same cells as their C. elegans orthologs, confirming gross conservation of regulatory mechanisms. However, the majority of exogenous elements, when placed in C. elegans, also directed expression in cells outside endogenous patterns, suggesting functional divergence. Recurrent ectopic expression of different promoters in the same C. elegans cells may reflect biases in the directions in which expression patterns can evolve due to shared regulatory logic of coexpressed genes. The fact that, despite differences between individual genes, several patterns repeatedly emerged from our survey, encourages us to think that general rules governing regulatory evolution may exist and be discoverable.

Transcriptome-wide profiling and expression analysis of transcription factor families in a liverwort, Marchantia polymorpha

BACKGROUND

Transcription factors (TFs) are vital elements that regulate transcription and the spatio-temporal expression of genes, thereby ensuring the accurate development and functioning of an organism. The identification of TF-encoding genes in a liverwort, Marchantia polymorpha, offers insights into TF organization in the members of the most basal lineages of land plants (embryophytes). Therefore, a comparison of Marchantia TF genes with other land plants (monocots, dicots, bryophytes) and algae (chlorophytes, rhodophytes) provides the most comprehensive view of the rates of expansion or contraction of TF genes in plant evolution.

RESULTS

In this study, we report the identification of TF-encoding transcripts in M. polymorpha for the first time, as evidenced by deep RNA sequencing data. In total, 3,471 putative TF encoding transcripts, distributed in 80 families, were identified, representing 7.4% of the generated Marchantia gametophytic transcriptome dataset. Overall, TF basic functions and distribution across families appear to be conserved when compared to other plant species. However, it is of interest to observe the genesis of novel sequences in 24 TF families and the apparent termination of 2 TF families with the emergence of Marchantia. Out of 24 TF families, 6 are known to be associated with plant reproductive development processes. We also examined the expression pattern of these TF-encoding transcripts in six male and female developmental stages in vegetative and reproductive gametophytic tissues of Marchantia.

CONCLUSIONS

The analysis highlighted the importance of Marchantia, a model plant system, in an evolutionary context. The dataset generated here provides a scientific resource for TF gene discovery and other comparative evolutionary studies of land plants.

The logic of gene regulatory networks in early vertebrate forebrain patterning

The vertebrate forebrain or prosencephalon is patterned at the beginning of neurulation into four major domains: the telencephalic, hypothalamic, retinal and diencephalic anlagen. These domains will then give rise to the majority of the brain structures involved in sensory integration and the control of higher intellectual and homeostatic functions. Understanding how forebrain pattering arises has thus attracted the interest of developmental neurobiologists for decades. As a result, most of its regulators have been identified and their hierarchical relationship is now the object of active investigation. Here, we summarize the main morphogenetic pathways and transcription factors involved in forebrain specification and propose the backbone of a possible gene regulatory network (GRN) governing its specification, taking advantage of the GRN principles elaborated by pioneer studies in simpler organisms. We will also discuss this GRN and its operational logic in the context of the remarkable morphological and functional diversification that the forebrain has undergone during evolution.

Sf-PHB2, a new transcription factor, drives WSSV Ie1 gene expression via a 12-bp DNA element

BACKGROUND

The WSSV immediate early gene ie1 is highly expressed throughout viral infection cycle and may play a central role in initiating viral replication during infection.

RESULTS

Here, a detailed characterization of the ie1 promoter was performed using deletion and mutation analyses to elucidate the role of the individual promoter motifs. Three results were obtained: 1) the ie1 promoter is a classical eukaryotic promoter that contains the initiator element (Inr) and TATA box responsible for the basal promoter activity; 2) mutation or truncation of a predicted Sp1 site decreased the level of promoter activity by about 3-fold, indicating that the Sp1 site is an important cis-element of the promoter; and 3) truncation of a 12-bp sequence that resides at -78/-67 of the ie1 promoter decreased the level of promoter activity by about 14-fold, indicating that the 12-bp motif is a critical upstream element of the ie1 promoter for binding of a strong transcription factor to drive the ie1 gene expression in the cells. Further, the 12-bp DNA binding protein was purified from the nuclear proteins of Sf9 cells using DNA affinity chromatography, and was identified as a homologue of the prohibitin2 protein (named as Sf-PHB2) using mass spectrometry. Furthermore, the DNA binding activity of Sf-PHB2 was verified using a super shift analysis.

CONCLUSION

These results support that the Sf-PHB2 is a novel transcription factor that drives WSSV ie1 gene expression by binding to the 12-bp DNA element.

Quantitative trait loci involved in regulating seed oil composition in Arabidopsis thaliana and their evolutionary implications

Fatty acid composition is an important determinant of seed oil quality. Overall, 72 QTL for 12 fatty acid traits that control seed oil composition were identified in four recombinant inbred line (RIL) populations (Ler-0 × Sha, Ler-0 × Col-4, Ler-2 × Cvi, Ler-0 × No-0) of Arabidopsis thaliana. The identified QTL explained 3.2-79.8% of the phenotypic variance; 33 of the 59 QTL identified in the Ler-0 × Sha and the Ler-0 × Col RIL populations co-located with several a priori candidate genes for seed oil composition. QTL for fatty acids 18:1, 18:2, 22:1, and fatty acids synthesized in plastids was identified in both Ler-0 × Sha and Ler-0 × Col-4 RIL populations, and QTL for 16:0 was identified in the Ler-0 × Sha and Ler-0 × No-0 RIL populations providing strong support for the importance of these QTL in determining seed oil composition. We identified melting point QTL in three RIL populations, and fatty acid QTL collocated with two of them, suggesting that the loci could be under selection for altering the melting point of seed oils to enhance adaptation and could be useful for breeding purposes. Nuclear-cytoplasmic interactions and epistasis were rare. Analysis of the genetic correlations between these loci and other fatty acids indicated that these correlations would tend to strongly enhance selection for desirable fatty acids.

Regulatory evolution through divergence of a phosphoswitch in the transcription factor CEBPB

There is an emerging consensus that gene regulation evolves through changes in cis-regulatory elements and transcription factors. Although it is clear how nucleotide substitutions in cis-regulatory elements affect gene expression, it is not clear how amino-acid substitutions in transcription factors influence gene regulation. Here we show that amino-acid changes in the transcription factor CCAAT/enhancer binding protein-β (CEBPB, also known as C/EBP-β) in the stem-lineage of placental mammals changed the way it responds to cyclic AMP/protein kinase A (cAMP/PKA) signalling. By functionally analysing resurrected ancestral proteins, we identify three amino-acid substitutions in an internal regulatory domain of CEBPB that are responsible for the novel function. These amino-acid substitutions reorganize the location of key phosphorylation sites, introducing a new site and removing two ancestral sites, reversing the response of CEBPB to GSK-3β-mediated phosphorylation from repression to activation. We conclude that changing the response of transcription factors to signalling pathways can be an important mechanism of gene regulatory evolution.

PS integrins and laminins: key regulators of cell migration during Drosophila embryogenesis

During embryonic development, there are numerous cases where organ or tissue formation depends upon the migration of primordial cells. In the Drosophila embryo, the visceral mesoderm (vm) acts as a substrate for the migration of several cell populations of epithelial origin, including the endoderm, the trachea and the salivary glands. These migratory processes require both integrins and laminins. The current model is that αPS1βPS (PS1) and/or αPS3βPS (PS3) integrins are required in migrating cells, whereas αPS2βPS (PS2) integrin is required in the vm, where it performs an as yet unidentified function. Here, we show that PS1 integrins are also required for the migration over the vm of cells of mesodermal origin, the caudal visceral mesoderm (CVM). These results support a model in which PS1 might have evolved to acquire the migratory function of integrins, irrespective of the origin of the tissue. This integrin function is highly specific and its specificity resides mainly in the extracellular domain. In addition, we have identified the Laminin α1,2 trimer, as the key extracellular matrix (ECM) component regulating CVM migration. Furthermore, we show that, as it is the case in vertebrates, integrins, and specifically PS2, contributes to CVM movement by participating in the correct assembly of the ECM that serves as tracks for migration.

Proteome-wide evidence for enhanced positive Darwinian selection within intrinsically disordered regions in proteins

BACKGROUND

Understanding the adaptive changes that alter the function of proteins during evolution is an important question for biology and medicine. The increasing number of completely sequenced genomes from closely related organisms, as well as individuals within species, facilitates systematic detection of recent selection events by means of comparative genomics.

RESULTS

We have used genome-wide strain-specific single nucleotide polymorphism data from 64 strains of budding yeast (Saccharomyces cerevisiae or Saccharomyces paradoxus) to determine whether adaptive positive selection is correlated with protein regions showing propensity for different classes of structure conformation. Data from phylogenetic and population genetic analysis of 3,746 gene alignments consistently shows a significantly higher degree of positive Darwinian selection in intrinsically disordered regions of proteins compared to regions of alpha helix, beta sheet or tertiary structure. Evidence of positive selection is significantly enriched in classes of proteins whose functions and molecular mechanisms can be coupled to adaptive processes and these classes tend to have a higher average content of intrinsically unstructured protein regions.

CONCLUSIONS

We suggest that intrinsically disordered protein regions may be important for the production and maintenance of genetic variation with adaptive potential and that they may thus be of central significance for the evolvability of the organism or cell in which they occur.

Divergence in cis-regulatory sequences surrounding the opsin gene arrays of African cichlid fishes

Background

Divergence within cis-regulatory sequences may contribute to the adaptive evolution of gene expression, but functional alleles in these regions are difficult to identify without abundant genomic resources. Among African cichlid fishes, the differential expression of seven opsin genes has produced adaptive differences in visual sensitivity. Quantitative genetic analysis suggests that cis-regulatory alleles near the SWS2-LWS opsins may contribute to this variation. Here, we sequence BACs containing the opsin genes of two cichlids, Oreochromis niloticus and Metriaclima zebra. We use phylogenetic footprinting and shadowing to examine divergence in conserved non-coding elements, promoter sequences, and 3'-UTRs surrounding each opsin in search of candidate cis-regulatory sequences that influence cichlid opsin expression.

Results

We identified 20 conserved non-coding elements surrounding the opsins of cichlids and other teleosts, including one known enhancer and a retinal microRNA. Most conserved elements contained computationally-predicted binding sites that correspond to transcription factors that function in vertebrate opsin expression; O. niloticus and M. zebra were significantly divergent in two of these. Similarly, we found a large number of relevant transcription factor binding sites within each opsin's proximal promoter, and identified five opsins that were considerably divergent in both expression and the number of transcription factor binding sites shared between O. niloticus and M. zebra. We also found several microRNA target sites within the 3'-UTR of each opsin, including two 3'-UTRs that differ significantly between O. niloticus and M. zebra. Finally, we examined interspecific divergence among 18 phenotypically diverse cichlids from Lake Malawi for one conserved non-coding element, two 3'-UTRs, and five opsin proximal promoters. We found that all regions were highly conserved with some evidence of CRX transcription factor binding site turnover. We also found three SNPs within two opsin promoters and one non-coding element that had weak association with cichlid opsin expression.

Conclusions

This study is the first to systematically search the opsins of cichlids for putative cis-regulatory sequences. Although many putative regulatory regions are highly conserved across a large number of phenotypically diverse cichlids, we found at least nine divergent sequences that could contribute to opsin expression differences in cis and stand out as candidates for future functional analyses.

Evolutionarily conserved function of Gbx2 in anterior hindbrain development

The amino acid sequence across the DNA-binding homeodomain of Gbx2 is highly conserved across multiple species. In mice, Gbx2 is essential for establishment of the midbrain-hindbrain boundary (MHB), and in development of anterior hindbrain structures, rhombomeres (r) 1-r3, and the r2/r3-derived cranial nerve V. In contrast, studies in zebrafish have implicated gbx1 in establishment of the MHB. Therefore, we tested potential roles for gbx2 in anterior hindbrain development in zebrafish. gbx2 knockdown with antisense morpholino results in increased cell death in r2, r3, and r5 and a truncation of the anterior hindbrain, similar to the defect in Gbx2(-/-) mice. Moreover, there is abnormal clustering of cranial nerve V cell bodies in r2 and r3 indicative of defects in aspects of anterior hindbrain patterning. These phenotypes can be rescued by expression of the mouse GBX2 protein. These results suggest that gbx2/Gbx2 has an evolutionarily conserved role in anterior hindbrain development.

A catalogue of eukaryotic transcription factor types, their evolutionary origin, and species distribution

Transcription factors (TFs) play key roles in the regulation of gene expression by binding in a sequence-specific manner to genomic DNA. In eukaryotes, DNA binding is achieved by a wide range of structural forms and motifs. TFs are typically classified by their DNA-binding domain (DBD) type. In this chapter, we catalogue and survey 91 different TF DBD types in metazoa, plants, fungi, and protists. We briefly discuss well-characterized TF families representing the major DBD superclasses. We also examine the species distributions and inferred evolutionary histories of the various families, and the potential roles played by TF family expansion and dimerization.

Surprising flexibility in a conserved Hox transcription factor over 550 million years of evolution

Although metazoan body plans are remarkably diverse, the structure and function of many embryonic regulatory genes are conserved because large changes would be detrimental to development. However, the fushi tarazu (ftz) gene has changed dramatically during arthropod evolution from Hox-like to a pair-rule segmentation gene in Drosophila. Changes in both expression and protein sequence contributed to this new function: ftz expression switched from Hox-like to stripes and changes in Ftz cofactor interaction motifs led to loss of homeotic and gain of segmentation potential. Here, we reconstructed ftz changes in a rigorous phylogenetic context. We found that ftz did not simply switch from Hox-like to segmentation function; rather, ftz is remarkably labile, having undergone multiple changes in sequence and expression. The segmentation LXXLL motif was stably acquired in holometabolous insects after the appearance of striped expression in early insect lineages. The homeotic YPWM motif independently degenerated multiple times. These "degen-YPWMs" showed varying degrees of homeotic potential when expressed in Drosophila, suggesting variable loss of Hox function in different arthropods. Finally, the intensity of ftz Hox-like expression decreased to marginal levels in some crustaceans. We propose that decreased expression levels permitted ftz variants to arise and persist in populations without disadvantaging organismal development. This process, in turn, allowed evolutionary transitions in protein function, as weakly expressed "hopeful gene variants" were coopted into alternative developmental pathways. Our findings show that variation of a pleiotropic transcription factor is more extensive than previously imagined, suggesting that evolutionary plasticity may be widespread among regulatory genes.

Silencing of an abdominal Hox gene during early development is correlated with limb development in a crustacean trunk

We tested whether Artemia abd-A could repress limbs in Drosophila embryos, and found that although abd-A transcripts were produced, ABD-A protein was not. Similarly, developing Artemia epidermal cells showed expression of abd-A transcripts without accumulation of ABD-A protein. This finding in Artemia reveals a new variation in Hox gene function that is associated with morphological evolution. In this case, a HOX protein expression pattern is completely absent during early development, although the HOX protein is expressed at later stages in the central nervous system in a "homeotic-like" pattern. The combination of an absence of ABD-A protein expression in the Artemia limb primordia and the weak repressive function of Artemia UBX protein on the limb-promoting gene Dll are likely to be two reasons why homonomous limbs develop throughout the entire Artemia trunk.

Genome-wide phylogenetic comparative analysis of plant transcriptional regulation: a timeline of loss, gain, expansion, and correlation with complexity

Evolutionary retention of duplicated genes encoding transcription-associated proteins (TAPs, comprising transcription factors and other transcriptional regulators) has been hypothesized to be positively correlated with increasing morphological complexity and paleopolyploidizations, especially within the plant kingdom. Here, we present the most comprehensive set of classification rules for TAPs and its application for genome-wide analyses of plants and algae. Using a dated species tree and phylogenetic comparative (PC) analyses, we define the timeline of TAP loss, gain, and expansion among Viridiplantae and find that two major bursts of gain/expansion occurred, coinciding with the water-to-land transition and the radiation of flowering plants. For the first time, we provide PC proof for the long-standing hypothesis that TAPs are major driving forces behind the evolution of morphological complexity, the latter in Plantae being shaped significantly by polyploidization and subsequent biased paleolog retention. Principal component analysis incorporating the number of TAPs per genome provides an alternate and significant proxy for complexity, ideally suited for PC genomics. Our work lays the ground for further interrogation of the shaping of gene regulatory networks underlying the evolution of organism complexity.

Expression patterns of AaDL, a CRABS CLAW ortholog in Asparagus asparagoides (Asparagaceae), demonstrate a stepwise evolution of CRC/DL subfamily of YABBY genes

Understanding the evolutionary path of transcription factors is essential for the elucidation of plant evolution. The CRC/DL subfamily of the YABBY gene family are functionally diverse, plant-specific, putative transcription factors. In Arabidopsis thaliana, CRABS CLAW (CRC) is expressed in the abaxial region of carpel primordia and in floral nectaries, where it regulates carpel morphology and nectary development. By contrast, in Oryza sativa, DROOPING LEAF (DL) is expressed in the entire carpel primordium and in the central undifferentiated cells of leaves, where it regulates carpel identity and midrib development. Recent studies suggest that abaxial expression and functional roles in the carpel are ancestral characters, although when and how neo-functionalizations occurred remains unclear. To elucidate the evolutionary processes of the CRC/DL subfamily, we examined in situ expression patterns of a CRC ortholog (AaDL) in Asparagus asparagoides (Asparagales). Like CRC in Arabidopsis thaliana, AaDL was clearly expressed in the abaxial region of the ovary wall. Expression was also detected in the phloem of leaves, but not in the septal nectary. Thus, expression in the entire carpel primordium might have been acquired after the divergence of Asparagus, with expression competence in the leaves acquired before the divergence of Asparagus in monocots. Our data indicate that the evolution of CRC/DL subfamily genes occurred in a stepwise manner.

Rapid sequence evolution of transcription factors controlling neuron differentiation in Caenorhabditis

Whether phenotypic evolution proceeds predominantly through changes in regulatory sequences is a controversial issue in evolutionary genetics. Ample evidence indicates that the evolution of gene regulatory networks via changes in cis-regulatory sequences is an important determinant of phenotypic diversity. However, recent experimental work suggests that the role of transcription factor (TF) divergence in developmental evolution may be underestimated. In order to help understand what levels of constraints are acting on the coding sequence of developmental regulatory genes, evolutionary rates were investigated among 48 TFs required for neuronal development in Caenorhabditis elegans. Allelic variation was then sampled for 28 of these genes within a population of the related species Caenorhabditis remanei. Neuronal TFs are more divergent, both within and between species, than structural genes. TFs affecting different neuronal classes are under different levels of selective constraints. The regulatory genes controlling the differentiation of chemosensory neurons evolve particularly fast and exhibit higher levels of within- and between-species nucleotide variation than TFs required for the development of several neuronal classes and TFs required for motorneuron differentiation. The TFs affecting chemosensory neuron development are also more divergent than chemosensory genes expressed in the neurons they differentiate. These results illustrate that TFs are not as highly constrained as commonly thought and suggest that the role of divergence in developmental regulatory genes during the evolution of gene regulatory networks requires further attention.

Transcription factor function and promoter architecture govern the evolution of bacterial regulons

Evolutionary changes in ancestral regulatory circuits can bring about phenotypic differences between related organisms. Studies of regulatory circuits in eukaryotes suggest that these modifications result primarily from changes in cis-regulatory elements (as opposed to alterations in the transcription factors that act upon these sequences). It is presently unclear how the evolution of gene regulatory circuits has proceeded in bacteria, given the rampant effects of horizontal gene transfer, which has significantly altered the composition of bacterial regulons. We now demonstrate that the evolution of the regulons governed by the regulatory protein PhoP in the related human pathogens Salmonella enterica and Yersinia pestis has entailed functional changes in the PhoP protein as well as in the architecture of PhoP-dependent promoters. These changes have resulted in orthologous PhoP proteins that differ both in their ability to promote transcription and in their role as virulence regulators. We posit that these changes allow bacterial transcription factors to incorporate newly acquired genes into ancestral regulatory circuits and yet retain control of the core members of a regulon.

Evolution of an insect-specific GROUCHO-interaction motif in the ENGRAILED selector protein

Animal morphology evolves through alterations in the genetic regulatory networks that control development. Regulatory connections are commonly added, subtracted, or modified via mutations in cis-regulatory elements, but several cases are also known where transcription factors have gained or lost activity-modulating peptide motifs. In order to better assess the role of novel transcription factor peptide motifs in evolution, we searched for synapomorphic motifs in the homeotic selectors of Drosophila melanogaster and related insects. Here, we describe an evolutionarily novel GROUCHO (GRO)-interaction motif in the ENGRAILED (EN) selector protein. This "ehIFRPF" motif is not homologous to the previously characterized "engrailed homology 1" (eh1) GRO-interaction motif of EN. This second motif is an insect-specific "WRPW"-type motif that has been maintained by purifying selection in at least the dipteran/lepidopteran lineage. We demonstrate that this motif contributes to in vivo repression of the wingless (wg) target gene and to interaction with GRO in vitro. The acquisition and conservation of this auxiliary peptide motif shows how the number and activity of short peptide motifs can evolve in transcription factors while existing regulatory functions are maintained.

The role of bZIP transcription factors in green plant evolution: adaptive features emerging from four founder genes

Background

Transcription factors of the basic leucine zipper (bZIP) family control important processes in all eukaryotes. In plants, bZIPs are regulators of many central developmental and physiological processes including photomorphogenesis, leaf and seed formation, energy homeostasis, and abiotic and biotic stress responses. Here we performed a comprehensive phylogenetic analysis of bZIP genes from algae, mosses, ferns, gymnosperms and angiosperms.

Methodology/Principal Findings

We identified 13 groups of bZIP homologues in angiosperms, three more than known before, that represent 34 Possible Groups of Orthologues (PoGOs). The 34 PoGOs may correspond to the complete set of ancestral angiosperm bZIP genes that participated in the diversification of flowering plants. Homologous genes dedicated to seed-related processes and ABA-mediated stress responses originated in the common ancestor of seed plants, and three groups of homologues emerged in the angiosperm lineage, of which one group plays a role in optimizing the use of energy.

Conclusions/Significance

Our data suggest that the ancestor of green plants possessed four bZIP genes functionally involved in oxidative stress and unfolded protein responses that are bZIP-mediated processes in all eukaryotes, but also in light-dependent regulations. The four founder genes amplified and diverged significantly, generating traits that benefited the colonization of new environments.

An empirical test for branch-specific positive selection

The use of phylogenetic analysis to predict positive selection specific to human genes is complicated by the very close evolutionary relationship with our nearest extant primate relatives, chimpanzees. To assess the power and limitations inherent in use of maximum-likelihood (ML) analysis of codon substitution patterns in such recently diverged species, a series of simulations was performed to assess the impact of several parameters of the evolutionary model on prediction of human-specific positive selection, including branch length and d(N)/d(S) ratio. Parameters were varied across a range of values observed in alignments of 175 transcription factor (TF) genes that were sequenced in 12 primate species. The ML method largely lacks the power to detect positive selection that has occurred since the most recent common ancestor between humans and chimpanzees. An alternative null model was developed on the basis of gene-specific evaluation of the empirical distribution of ML results, using simulated neutrally evolving sequences. This empirical test provides greater sensitivity to detect lineage-specific positive selection in the context of recent evolutionary divergence.

Resurrecting the role of transcription factor change in developmental evolution

A long-standing question in evolutionary and developmental biology concerns the relative contribution of cis-regulatory and protein changes to developmental evolution. Central to this argument is which mutations generate evolutionarily relevant phenotypic variation? A review of the growing body of evolutionary and developmental literature supports the notion that many developmentally relevant differences occur in the cis-regulatory regions of protein-coding genes, generally to the exclusion of changes in the protein-coding region of genes. However, accumulating experimental evidence demonstrates that many of the arguments against a role for proteins in the evolution of gene regulation, and the developmental evolution in general, are no longer supported and there is an increasing number of cases in which transcription factor protein changes have been demonstrated in evolution. Here, we review the evidence that cis-regulatory evolution is an important driver of phenotypic evolution and provide examples of protein-mediated developmental evolution. Finally, we present an argument that the evolution of proteins may play a more substantial, but thus far underestimated, role in developmental evolution.

cDNA sequences for transcription factors and signaling proteins of the hemichordate Saccoglossus kowalevskii: efficacy of the expressed sequence tag (EST) approach for evolutionary and developmental studies of a new organism

We describe a collection of expressed sequence tags (ESTs) for Saccoglossus kowalevskii, a direct-developing hemichordate valuable for evolutionary comparisons with chordates. The 202,175 ESTs represent 163,633 arrayed clones carrying cDNAs prepared from embryonic libraries, and they assemble into 13,677 continuous sequences (contigs), leaving 10,896 singletons (excluding mitochondrial sequences). Of the contigs, 53% had significant matches when BLAST was used to query the NCBI databases (< or = 10(-10)), as did 51% of the singletons. Contigs most frequently matched sequences from amphioxus (29%), chordates (67%), and deuterostomes (87%). From the clone array, we isolated 400 full-length sequences for transcription factors and signaling proteins of use for evolutionary and developmental studies. The set includes sequences for fox, pax, tbx, hox, and other homeobox-containing factors, and for ligands and receptors of the TGFbeta, Wnt, Hh, Delta/Notch, and RTK pathways. At least 80% of key sequences have been obtained, when judged against gene lists of model organisms. The median length of these cDNAs is 2.3 kb, including 1.05 kb of 3' untranslated region (UTR). Only 30% are entirely matched by single contigs assembled from ESTs. We conclude that an EST collection based on 150,000 clones is a rich source of sequences for molecular developmental work, and that the EST approach is an efficient way to initiate comparative studies of a new organism.

Context-dependent regulation of Hox protein functions by CK2 phosphorylation sites

Variations in Hox protein sequences and functions have been proposed to contribute to evolutionary changes in appendage shape and number in crustaceans and insects. One model is that insect Hox proteins of the Ultrabithorax (UBX) ortholog class evolved increased abilities to repress Distal-less (Dll) transcription and appendage development in part through the loss of serine and threonine residues in casein kinase 2 (CK2) phosphorylation sites. To explore this possibility, we constructed and tested the appendage repression function of chimeric proteins with insertions of different CK2 consensus sites or phosphomimetics of CK2 sites in C-terminal regions of Drosophila melanogaster UBX. Our results indicate that CK2 sites C-terminal to the homeodomain can inhibit the appendage repression functions of UBX proteins, but only in the context of specific amino acid sequences. Our results, combined with previous findings on evolutionary changes in Hox protein, suggest how intra-protein regulatory changes can diversify Hox protein function, and thus animal morphology.

Proximal-distal patterns of transcription factor gene expression during Arabidopsis root development

The expression pattern of genes can identify the cells in which the respective proteins are active during development. As a step towards defining the genetic network that controls the development of roots, a high-throughput method of whole-mount in situ hybridization has been developed that does not require expensive equipment and allows the definition of the expression patterns of 137 transcription factor genes in young developing roots. Of the 137 transcription factors, 81.8% were expressed in the root while 18.2% showed no detectable expression. In all three proximal distal zones (meristem, elongation, and differentiation) of the root, 52.6% were expressed whereas 21.2% were expressed in only two zones. Eight percent of the genes were expressed in a single proximal distal zone. Cell-specific gene expression patterns were also detected. This rapid approach identified potential key regulators of cell differentiation and provides important spatial information for the expression patterns of a large number of transcriptional regulators that function during root development.

Genome clashes in hybrids: insights from gene expression

In interspecific hybrids, novel phenotypes often emerge from the interaction of two divergent genomes. Interactions between the two transcriptional networks are assumed to contribute to these unpredicted new phenotypes by inducing novel patterns of gene expression. Here we provide a review of the recent literature on the accumulation of regulatory incompatibilities. We review specific examples of regulatory incompatibilities reported at particular loci as well as genome-scale surveys of gene expression in interspecific hybrids. Finally, we consider and preview novel technologies that could help decipher how divergent transcriptional networks interact in hybrids between species.

The locus of evolution: evo devo and the genetics of adaptation

An important tenet of evolutionary developmental biology ("evo devo") is that adaptive mutations affecting morphology are more likely to occur in the cis-regulatory regions than in the protein-coding regions of genes. This argument rests on two claims: (1) the modular nature of cis-regulatory elements largely frees them from deleterious pleiotropic effects, and (2) a growing body of empirical evidence appears to support the predominant role of gene regulatory change in adaptation, especially morphological adaptation. Here we discuss and critique these assertions. We first show that there is no theoretical or empirical basis for the evo devo contention that adaptations involving morphology evolve by genetic mechanisms different from those involving physiology and other traits. In addition, some forms of protein evolution can avoid the negative consequences of pleiotropy, most notably via gene duplication. In light of evo devo claims, we then examine the substantial data on the genetic basis of adaptation from both genome-wide surveys and single-locus studies. Genomic studies lend little support to the cis-regulatory theory: many of these have detected adaptation in protein-coding regions, including transcription factors, whereas few have examined regulatory regions. Turning to single-locus studies, we note that the most widely cited examples of adaptive cis-regulatory mutations focus on trait loss rather than gain, and none have yet pinpointed an evolved regulatory site. In contrast, there are many studies that have both identified structural mutations and functionally verified their contribution to adaptation and speciation. Neither the theoretical arguments nor the data from nature, then, support the claim for a predominance of cis-regulatory mutations in evolution. Although this claim may be true, it is at best premature. Adaptation and speciation probably proceed through a combination of cis-regulatory and structural mutations, with a substantial contribution of the latter.