Reversal of Hox1 gene subfunctionalization in the mouse

Hox genes and patterning the vertebrate body

The diversity of vertebrate body plans is dizzying, yet stunning for the many things they have in common. Vertebrates have inhabited virtually every part of the earth from its coldest to warmest climates. They locomote by swimming, flying, walking, slithering, or climbing, or combinations of these behaviors. And they exist in many different sizes, from the smallest of frogs, fish and lizards to giraffes, elephants, and blue whales. Despite these differences, vertebrates follow a remarkably similar blueprint for the establishment of their body plan. Within the relatively small amount of time required to complete gastrulation, the process through which the three germ layers, ectoderm, mesoderm, and endoderm are created, the embryo also generates its body axis and is simultaneously patterned. For the length of this axis, the genes that distinguish the neck from the rib cage or the trunk from the sacrum are the Hox genes. In vertebrates, there was evolutionary pressure to maintain this set of genes in the organism. Over the past decades, much has been learned regarding the regulatory mechanisms that ensure the appropriate expression of these genes along the main body axes. Genetic functions continue to be explored though much has been learned. Much less has been discerned on the identity of co-factors used by Hox proteins for the specificity of transcriptional regulation or what downstream targets and pathways are critical for patterning events, though there are notable exceptions. Current work in the field is demonstrating that Hox genes continue to function in many organs long after directing early patterning events. It is hopeful continued research will shed light on remaining questions regarding mechanisms used by this important and conserved set of transcriptional regulators.

Hox genes in development and beyond

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Essay the (unusual) heuristic value of Hox gene clusters; a matter of time?

Ever since their first report in 1984, Antennapedia-type homeobox (Hox) genes have been involved in such a series of interesting observations, in particular due to their conserved clustered organization between vertebrates and arthropods, that one may legitimately wonder about the origin of this heuristic value. In this essay, I first consider different examples where Hox gene clusters have been instrumental in providing conceptual advances, taken from various fields of research and mostly involving vertebrate embryos. These examples touch upon our understanding of genomic evolution, the revisiting of 19th century views on the relationships between development and evolution and the building of a new framework to understand long-range and pleiotropic gene regulation during development. I then discuss whether the high value of the Hox gene family, when considered as an epistemic object, is related to its clustered structure (and the absence thereof in some animal species) and, if so, what is it in such particular genetic oddities that made them so generous in providing the scientific community with interesting information.

Segmentation and patterning of the vertebrate hindbrain

During early development, the hindbrain is sub-divided into rhombomeres that underlie the organisation of neurons and adjacent craniofacial tissues. A gene regulatory network of signals and transcription factors establish and pattern segments with a distinct anteroposterior identity. Initially, the borders of segmental gene expression are imprecise, but then become sharply defined, and specialised boundary cells form. In this Review, we summarise key aspects of the conserved regulatory cascade that underlies the formation of hindbrain segments. We describe how the pattern is sharpened and stabilised through the dynamic regulation of cell identity, acting in parallel with cell segregation. Finally, we discuss evidence that boundary cells have roles in local patterning, and act as a site of neurogenesis within the hindbrain.

Dynamic Molecular Evolution of Mammalian Homeobox Genes: Duplication, Loss, Divergence and Gene Conversion Sculpt PRD Class Repertoires

The majority of homeobox genes are highly conserved across animals, but the eutherian-specific ETCHbox genes, embryonically expressed and highly divergent duplicates of CRX, are a notable exception. Here we compare the ETCHbox genes of 34 mammalian species, uncovering dynamic patterns of gene loss and tandem duplication, including the presence of a large tandem array of LEUTX loci in the genome of the European rabbit (Oryctolagus cuniculus). Despite extensive gene gain and loss, all sampled species possess at least two ETCHbox genes, suggesting their collective role is indispensable. We find evidence for positive selection and show that TPRX1 and TPRX2 have been the subject of repeated gene conversion across the Boreoeutheria, homogenising their sequences and preventing divergence, especially in the homeobox region. Together, these results are consistent with a model where mammalian ETCHbox genes are dynamic in evolution due to functional overlap, yet have collective indispensable roles.

The Hox protein conundrum: The "specifics" of DNA binding for Hox proteins and their partners

Homeotic genes (Hox genes) are homeodomain-transcription factors involved in conferring segmental identity along the anterior-posterior body axis. Molecular characterization of HOX protein function raises some interesting questions regarding the source of the binding specificity of the HOX proteins. How do HOX proteins regulate common and unique target specificity across space and time? This review attempts to summarize and interpret findings in this area, largely focused on results from in vitro and in vivo studies in Drosophila and mouse systems. Recent studies related to HOX protein binding specificity compel us to reconsider some of our current models for transcription factor-DNA interactions. It is crucial to study transcription factor binding by incorporating components of more complex, multi-protein interactions in concert with small changes in binding motifs that can significantly impact DNA binding specificity and subsequent alterations in gene expression. To incorporate the multiple elements that can determine HOX protein binding specificity, we propose a more integrative Cooperative Binding model.

Diverse Molecular Mechanisms Contribute to Differential Expression of Human Duplicated Genes

Emerging evidence links genes within human-specific segmental duplications (HSDs) to traits and diseases unique to our species. Strikingly, despite being nearly identical by sequence (>98.5%), paralogous HSD genes are differentially expressed across human cell and tissue types, though the underlying mechanisms have not been examined. We compared cross-tissue mRNA levels of 75 HSD genes from 30 families between humans and chimpanzees and found expression patterns consistent with relaxed selection on or neofunctionalization of derived paralogs. In general, ancestral paralogs exhibited greatest expression conservation with chimpanzee orthologs, though exceptions suggest certain derived paralogs may retain or supplant ancestral functions. Concordantly, analysis of long-read isoform sequencing data sets from diverse human tissues and cell lines found that about half of derived paralogs exhibited globally lower expression. To understand mechanisms underlying these differences, we leveraged data from human lymphoblastoid cell lines (LCLs) and found no relationship between paralogous expression divergence and post-transcriptional regulation, sequence divergence, or copy-number variation. Considering cis-regulation, we reanalyzed ENCODE data and recovered hundreds of previously unidentified candidate CREs in HSDs. We also generated large-insert ChIP-sequencing data for active chromatin features in an LCL to better distinguish paralogous regions. Some duplicated CREs were sufficient to drive differential reporter activity, suggesting they may contribute to divergent cis-regulation of paralogous genes. This work provides evidence that cis-regulatory divergence contributes to novel expression patterns of recent gene duplicates in humans.

Autosomal recessive diseases among the Athabaskans of the southwestern United States: anthropological, medical, and scientific aspects

The peopling of the Americas by Native Americans occurred in 4 waves of which the last was Nadene language speakers of whom Athabaskans are the largest group. As the Europeans were entering the Southwestern states of the USA, Athabaskan hunting-gathering tribes were migrating South from Canada along the Rocky Mountains and undergoing potential bottlenecks reflected in autosomal recessive diseases shared by Apaches and Navajos. About 300 years ago, the Navajo developing a sedentary culture learned from Pueblo Indians while the Apache remained hunter-gathers. Although most of the tribe was rounded up and forced to relocate to Bosque Redondo, the adult breeding population was large enough to prevent a genetic bottleneck. However, some Navajo underwent further population bottlenecks while hiding from the brutal US Army action (under Kit Carson’s guidance). This led to an increased frequency of other autosomal recessive diseases. Recent advances in population genetics, pathophysiology of the diseases, and social/ethical issues concerning their study are reviewed.

A six-amino-acid motif is a major determinant in functional evolution of HOX1 proteins

Gene duplication and divergence is a major driver in the emergence of evolutionary novelties. How variations in amino acid sequences lead to loss of ancestral activity and functional diversification of proteins is poorly understood. We used cross-species functional analysis of Drosophila Labial and its mouse HOX1 orthologs (HOXA1, HOXB1, and HOXD1) as a paradigm to address this issue. Mouse HOX1 proteins display low (30%) sequence similarity with Drosophila Labial. However, substituting endogenous Labial with the mouse proteins revealed that HOXA1 has retained essential ancestral functions of Labial, while HOXB1 and HOXD1 have diverged. Genome-wide analysis demonstrated similar DNA-binding patterns of HOXA1 and Labial in mouse cells, while HOXB1 binds to distinct targets. Compared with HOXB1, HOXA1 shows an enrichment in co-occupancy with PBX proteins on target sites and exists in the same complex with PBX on chromatin. Functional analysis of HOXA1-HOXB1 chimeric proteins uncovered a novel six-amino-acid C-terminal motif (CTM) flanking the homeodomain that serves as a major determinant of ancestral activity. In vitro DNA-binding experiments and structural prediction show that CTM provides an important domain for interaction of HOXA1 proteins with PBX. Our findings show that small changes outside of highly conserved DNA-binding regions can lead to profound changes in protein function.

GATA6 mutations in hiPSCs inform mechanisms for maldevelopment of the heart, pancreas, and diaphragm

Damaging GATA6 variants cause cardiac outflow tract defects, sometimes with pancreatic and diaphragmic malformations. To define molecular mechanisms for these diverse developmental defects, we studied transcriptional and epigenetic responses to GATA6 loss of function (LoF) and missense variants during cardiomyocyte differentiation of isogenic human induced pluripotent stem cells. We show that GATA6 is a pioneer factor in cardiac development, regulating SMYD1 that activates HAND2, and KDR that with HAND2 orchestrates outflow tract formation. LoF variants perturbed cardiac genes and also endoderm lineage genes that direct PDX1 expression and pancreatic development. Remarkably, an exon 4 GATA6 missense variant, highly associated with extra-cardiac malformations, caused ectopic pioneer activities, profoundly diminishing GATA4, FOXA1/2, and PDX1 expression and increasing normal retinoic acid signaling that promotes diaphragm development. These aberrant epigenetic and transcriptional signatures illuminate the molecular mechanisms for cardiovascular malformations, pancreas and diaphragm dysgenesis that arise in patients with distinct GATA6 variants.

Unexpected mutual regulation underlies paralogue functional diversification and promotes epithelial tissue maturation in Tribolium

Insect Hox3/zen genes represent an evolutionary hotspot for changes in function and copy number. Single orthologues are required either for early specification or late morphogenesis of the extraembryonic tissues, which protect the embryo. The tandemly duplicated zen paralogues of the beetle Tribolium castaneum present a unique opportunity to investigate both functions in a single species. We dissect the paralogues' expression dynamics (transcript and protein) and transcriptional targets (RNA-seq after RNAi) throughout embryogenesis. We identify an unexpected role of Tc-Zen2 in repression of Tc-zen1, generating a negative feedback loop that promotes developmental progression. Tc-Zen2 regulation is dynamic, including within co-expressed multigene loci. We also show that extraembryonic development is the major event within the transcriptional landscape of late embryogenesis and provide a global molecular characterization of the extraembryonic serosal tissue. Altogether, we propose that paralogue mutual regulation arose through multiple instances of zen subfunctionalization, leading to their complementary extant roles.

Segmental arithmetic: summing up the Hox gene regulatory network for hindbrain development in chordates

Organization and development of the early vertebrate hindbrain are controlled by a cascade of regulatory interactions that govern the process of segmentation and patterning along the anterior-posterior axis via Hox genes. These interactions can be assembled into a gene regulatory network that provides a framework to interpret experimental data, generate hypotheses, and identify gaps in our understanding of the progressive process of hindbrain segmentation. The network can be broadly separated into a series of interconnected programs that govern early signaling, segmental subdivision, secondary signaling, segmentation, and ultimately specification of segmental identity. Hox genes play crucial roles in multiple programs within this network. Furthermore, the network reveals properties and principles that are likely to be general to other complex developmental systems. Data from vertebrate and invertebrate chordate models are shedding light on the origin and diversification of the network. Comprehensive cis-regulatory analyses of vertebrate Hox gene regulation have enabled powerful cross-species gene regulatory comparisons. Such an approach in the sea lamprey has revealed that the network mediating segmental Hox expression was present in ancestral vertebrates and has been maintained across diverse vertebrate lineages. Invertebrate chordates lack hindbrain segmentation but exhibit conservation of some aspects of the network, such as a role for retinoic acid in establishing nested Hox expression domains. These comparisons lead to a model in which early vertebrates underwent an elaboration of the network between anterior-posterior patterning and Hox gene expression, leading to the gene-regulatory programs for segmental subdivision and rhombomeric segmentation. WIREs Dev Biol 2017, 6:e286. doi: 10.1002/wdev.286 For further resources related to this article, please visit the WIREs website.

Mouse fitness measures reveal incomplete functional redundancy of Hox paralogous group 1 proteins

Here we assess the fitness consequences of the replacement of the Hoxa1 coding region with its paralog Hoxb1 in mice (Mus musculus) residing in semi-natural enclosures. Previously, this Hoxa1^B1 swap was reported as resulting in no discernible embryonic or physiological phenotype (i.e., functionally redundant), despite the 51% amino acid sequence differences between these two Hox proteins. Within heterozygous breeding cages no differences in litter size nor deviations from Mendelian genotypic expectations were observed in the outbred progeny; however, within semi-natural population enclosures mice homozygous for the Hoxa1^B1 swap were out-reproduced by controls resulting in the mutant allele being only 87.5% as frequent as the control in offspring born within enclosures. Specifically, Hoxa1^B1 founders produced only 77.9% as many offspring relative to controls, as measured by homozygous pups, and a 22.1% deficiency of heterozygous offspring was also observed. These data suggest that Hoxa1 and Hoxb1 have diverged in function through either sub- or neo-functionalization and that the HoxA1 and HoxB1 proteins are not mutually interchangeable when expressed from the Hoxa1 locus. The fitness assays conducted under naturalistic conditions in this study have provided an ultimate-level assessment of the postulated equivalence of competing alleles. Characterization of these differences has provided greater understanding of the forces shaping the maintenance and diversifications of Hox genes as well as other paralogous genes. This fitness assay approach can be applied to any genetic manipulation and often provides the most sensitive way to detect functional differences.

New Insights Into the Roles of Retinoic Acid Signaling in Nervous System Development and the Establishment of Neurotransmitter Systems

Secreted chiefly from the underlying mesoderm, the morphogen retinoic acid (RA) is well known to contribute to the specification, patterning, and differentiation of neural progenitors in the developing vertebrate nervous system. Furthermore, RA influences the subtype identity and neurotransmitter phenotype of subsets of maturing neurons, although relatively little is known about how these functions are mediated. This review provides a comprehensive overview of the roles played by RA signaling during the formation of the central and peripheral nervous systems of vertebrates and highlights its effects on the differentiation of several neurotransmitter systems. In addition, the evolutionary history of the RA signaling system is discussed, revealing both conserved properties and alternate modes of RA action. It is proposed that comparative approaches should be employed systematically to expand our knowledge of the context-dependent cellular mechanisms controlled by the multifunctional signaling molecule RA.

Fitness Assays Reveal Incomplete Functional Redundancy of the HoxA1 and HoxB1 Paralogs of Mice

Gene targeting techniques have led to the phenotypic characterization of numerous genes; however, many genes show minimal to no phenotypic consequences when disrupted, despite many having highly conserved sequences. The standard explanation for these findings is functional redundancy. A competing hypothesis is that these genes have important ecological functions in natural environments that are not needed under laboratory settings. Here we discriminate between these hypotheses by competing mice (Mus musculus) whose Hoxb1 gene has been replaced by Hoxa1, its highly conserved paralog, against matched wild-type controls in seminatural enclosures. This Hoxb1(A1) swap was reported as a genetic manipulation resulting in no discernible embryonic or physiological phenotype under standard laboratory tests. We observed a transient decline in first litter size for Hoxb1(A1) homozygous mice in breeding cages, but their fitness was consistently and more dramatically reduced when competing against controls within seminatural populations. Specifically, males homozygous for the Hoxb1(A1) swap acquired 10.6% fewer territories and the frequency of the Hoxb1(A1) allele decreased from 0.500 in population founders to 0.419 in their offspring. The decrease in Hoxb1(A1) frequency corresponded with a deficiency of both Hoxb1(A1) homozygous and heterozygous offspring. These data suggest that Hoxb1 and Hoxa1 are more phenotypically divergent than previously reported and support that sub- and/or neofunctionalization has occurred in these paralogous genes leading to a divergence of gene function and incomplete redundancy. Furthermore, this study highlights the importance of obtaining fitness measures of mutants in ecologically relevant conditions to better understand gene function and evolution.

A distal 594 bp ECR specifies Hmx1 expression in pinna and lateral facial morphogenesis and is regulated by the Hox-Pbx-Meis complex

Hmx1 encodes a homeodomain transcription factor expressed in the developing lateral craniofacial mesenchyme, retina and sensory ganglia. Mutation or mis-regulation of Hmx1 underlies malformations of the eye and external ear in multiple species. Deletion or insertional duplication of an evolutionarily conserved region (ECR) downstream of Hmx1 has recently been described in rat and cow, respectively. Here, we demonstrate that the impact of Hmx1 loss is greater than previously appreciated, with a variety of lateral cranioskeletal defects, auriculofacial nerve deficits, and duplication of the caudal region of the external ear. Using a transgenic approach, we demonstrate that a 594 bp sequence encompassing the ECR recapitulates specific aspects of the endogenous Hmx1 lateral facial expression pattern. Moreover, we show that Hoxa2, Meis and Pbx proteins act cooperatively on the ECR, via a core 32 bp sequence, to regulate Hmx1 expression. These studies highlight the conserved role for Hmx1 in BA2-derived tissues and provide an entry point for improved understanding of the causes of the frequent lateral facial birth defects in humans.

Hoxa5: A Key Player in Development and Disease

A critical position in the developmental hierarchy is occupied by the Hox genes, which encode transcription factors. Hox genes are crucial in specifying regional identity along the embryonic axes and in regulating morphogenesis. In mouse, targeted mutations of Hox genes cause skeletal transformations and organ defects that can impair viability. Here, we present the current knowledge about the Hoxa5 gene, a paradigm for the function and the regulation of Hox genes. The phenotypic survey of Hoxa5-/- mice has unveiled its critical role in the regional specification of the skeleton and in organogenesis. Most Hoxa5-/- mice die at birth from respiratory distress due to tracheal and lung dysmorphogenesis and impaired diaphragm innervation. The severity of the phenotype establishes that Hoxa5 plays a predominant role in lung organogenesis versus other Hox genes. Hoxa5 also governs digestive tract morphogenesis, thyroid and mammary glands development, and ovary homeostasis. Deregulated Hoxa5 expression is reported in cancers, indicating Hoxa5 involvement in tumor predisposition and progression. The dynamic Hoxa5 expression profile is under the transcriptional control of multiple cis-acting sequences and trans-acting regulators. It is also modulated by epigenetic mechanisms, implicating chromatin modifications and microRNAs. Finally, lncRNAs originating from alternative splicing and distal promoters encompass the Hoxa5 locus.

Cellular and molecular insights into Hox protein action

Hox genes encode homeodomain transcription factors that control morphogenesis and have established functions in development and evolution. Hox proteins have remained enigmatic with regard to the molecular mechanisms that endow them with specific and diverse functions, and to the cellular functions that they control. Here, we review recent examples of Hox-controlled cellular functions that highlight their versatile and highly context-dependent activity. This provides the setting to discuss how Hox proteins control morphogenesis and organogenesis. We then summarise the molecular modalities underlying Hox protein function, in particular in light of current models of transcription factor function. Finally, we discuss how functional divergence between Hox proteins might be achieved to give rise to the many facets of their action.

Shox2 is required for the proper development of the facial motor nucleus and the establishment of the facial nerves

Background

Axons from the visceral motor neurons (vMNs) project from nuclei in the hindbrain to innervate autonomic ganglia and branchial arch-derived muscles. Although much is known about the events that govern specification of somatic motor neurons, the genetic pathways responsible for the development of vMNs are less well characterized. We know that vMNs, like all motor neurons, depend on sonic hedgehog signaling for their generation. Similarly, the paired-like homeobox 2b (Phox2b) gene, which is expressed in both proliferating progenitors and post-mitotic motor neurons, is essential for the development of vMNs. Given that our previous study identified a novel role for the short stature homeobox 2 (Shox2) gene in the hindbrain, and since SHOX2 has been shown to regulate transcription of islet 1 (Isl1), an important regulator of vMN development, we sought to determine whether Shox2 is required for the proper development of the facial motor nucleus.

Results

Using a Nestin-Cre driver, we show that elimination of Shox2 throughout the brain results in elevated cell death in the facial motor nucleus at embryonic day 12.5 (E12.5) and E14.5, which correlates with impaired axonal projection properties of vMNs. We also observed changes in the spatial expression of the vMN cell fate factors Isl1 and Phox2b, and concomitant defects in Shh and Ptch1 expression in Shox2 mutants. Furthermore, we demonstrate that elimination of Shox2 results in the loss of dorsomedial and ventromedial subnuclei by postnatal day 0 (P0), which may explain the changes in physical activity and impaired feeding/nursing behavior in Shox2 mutants.

Conclusions

Combined, our data show that Shox2 is required for development of the facial motor nucleus and its associated facial (VII) nerves, and serves as a new molecular tool to probe the genetic programs of this understudied hindbrain region.

Electronic supplementary material

The online version of this article (doi:10.1186/s12868-015-0176-0) contains supplementary material, which is available to authorized users.

Imaging activity in neurons and glia with a Polr2a-based and cre-dependent GCaMP5G-IRES-tdTomato reporter mouse

New strategies for introducing genetically encoded activity indicators into animal models facilitate the investigation of nervous system function. We have developed the PC::G5-tdT mouse line that expresses the GCaMP5G calcium indicator in a Cre-dependent fashion. Instead of targeting the ROSA26 locus, we inserted the reporter cassette nearby the ubiquitously expressed Polr2a gene without disrupting locus integrity. The indicator was tagged with IRES-tdTomato to aid detection of positive cells. This reporter system is effective in a wide range of developmental and cellular contexts. We recorded spontaneous cortical calcium waves in intact awake newborns and evaluated concentration-dependent responses to odorants in the adult olfactory bulb. Moreover, PC::G5-tdT effectively reports intracellular calcium dynamics in somas and fine processes of astrocytes and microglial cells. Through electrophysiological and behavioral analyses, we determined that GCaMP5G expression had no major impact on nervous system performance. PC::G5-tdT will be instrumental for a variety of brain mapping experiments.

VIDEO ABSTRACT

Expression evolution facilitated the convergent neofunctionalization of a sodium channel gene

Ion channels have played a substantial role in the evolution of novel traits across all of the domains of life. A fascinating example of a novel adaptation is the convergent evolution of electric organs in the Mormyroid and Gymnotiform electric fishes. The regulated currents that flow through ion channels directly generate the electrical signals which have evolved in these fish. Here, we investigated how the expression evolution of two sodium channel paralogs (Scn4aa and Scn4ab) influenced their convergent molecular evolution following the teleost-specific whole-genome duplication. We developed a reliable assay to accurately measure the expression stoichiometry of these genes and used this technique to analyze relative expression of the duplicate genes in a phylogenetic context. We found that before a major shift in expression from skeletal muscle and neofunctionalization in the muscle-derived electric organ, Scn4aa was first downregulated in the ancestors of both electric lineages. This indicates that underlying the convergent evolution of this gene, there was a greater propensity toward neofunctionalization due to its decreased expression relative to its paralog Scn4ab. We investigated another derived muscle tissue, the sonic organ of Porichthys notatus, and show that, as in the electric fishes, Scn4aa again shows a radical shift in expression away from the ancestral muscle cells into the evolutionarily novel muscle-derived tissue. This study presents evidence that expression downregulation facilitates neofunctionalization after gene duplication, a pattern that may often set the stage for novel trait evolution after gene duplication.

Hox targets and cellular functions

Hox genes are a group of genes that specify structures along the anteroposterior axis in bilaterians. Although in many cases they do so by modifying a homologous structure with a different (or no) Hox input, there are also examples of Hox genes constructing new organs with no homology in other regions of the body. Hox genes determine structures though the regulation of targets implementing cellular functions and by coordinating cell behavior. The genetic organization to construct or modify a certain organ involves both a genetic cascade through intermediate transcription factors and a direct regulation of targets carrying out cellular functions. In this review I discuss new data from genome-wide techniques, as well as previous genetic and developmental information, to describe some examples of Hox regulation of different cell functions. I also discuss the organization of genetic cascades leading to the development of new organs, mainly using Drosophila melanogaster as the model to analyze Hox function.

Evolution of transcriptional enhancers and animal diversity

Deciphering the genetic bases that drive animal diversity is one of the major challenges of modern biology. Although four decades ago it was proposed that animal evolution was mainly driven by changes in cis-regulatory DNA elements controlling gene expression rather than in protein-coding sequences, only now are powerful bioinformatics and experimental approaches available to accelerate studies into how the evolution of transcriptional enhancers contributes to novel forms and functions. In the introduction to this Theme Issue, we start by defining the general properties of transcriptional enhancers, such as modularity and the coexistence of tight sequence conservation with transcription factor-binding site shuffling as different mechanisms that maintain the enhancer grammar over evolutionary time. We discuss past and current methods used to identify cell-type-specific enhancers and provide examples of how enhancers originate de novo, change and are lost in particular lineages. We then focus in the central part of this Theme Issue on analysing examples of how the molecular evolution of enhancers may change form and function. Throughout this introduction, we present the main findings of the articles, reviews and perspectives contributed to this Theme Issue that together illustrate some of the great advances and current frontiers in the field.

Hox genes: choreographers in neural development, architects of circuit organization

The neural circuits governing vital behaviors, such as respiration and locomotion, are comprised of discrete neuronal populations residing within the brainstem and spinal cord. Work over the past decade has provided a fairly comprehensive understanding of the developmental pathways that determine the identity of major neuronal classes within the neural tube. However, the steps through which neurons acquire the subtype diversities necessary for their incorporation into a particular circuit are still poorly defined. Studies on the specification of motor neurons indicate that the large family of Hox transcription factors has a key role in generating the subtypes required for selective muscle innervation. There is also emerging evidence that Hox genes function in multiple neuronal classes to shape synaptic specificity during development, suggesting a broader role in circuit assembly. This Review highlights the functions and mechanisms of Hox gene networks and their multifaceted roles during neuronal specification and connectivity.

Partial functional redundancy between Hoxa5 and Hoxb5 paralog genes during lung morphogenesis

Hox genes encode transcription factors governing complex developmental processes in several organs. A subset of Hox genes are expressed in the developing lung. Except for Hoxa5, the lack of overt lung phenotype in single mutants suggests that Hox genes may not play a predominant role in lung ontogeny or that functional redundancy may mask anomalies. In the Hox5 paralog group, both Hoxa5 and Hoxb5 genes are expressed in the lung mesenchyme whereas Hoxa5 is also expressed in the tracheal mesenchyme. Herein, we generated Hoxa5;Hoxb5 compound mutant mice to evaluate the relative contribution of each gene to lung development. Hoxa5;Hoxb5 mutants carrying the four mutated alleles displayed an aggravated lung phenotype, resulting in the death of the mutant pups at birth. Characterization of the phenotype highlighted the role of Hoxb5 in lung formation, the latter being involved in branching morphogenesis, goblet cell specification, and postnatal air space structure, revealing partial functional redundancy with Hoxa5. However, the Hoxb5 lung phenotypes were less severe than those seen in Hoxa5 mutants, likely because of Hoxa5 compensation. New specific roles for Hoxa5 were also unveiled, demonstrating the extensive contribution of Hoxa5 to the developing respiratory system. The exclusive expression of Hoxa5 in the trachea and the phrenic motor column likely underlies the Hoxa5-specific trachea and diaphragm phenotypes. Altogether, our observations establish that the Hoxa5 and Hoxb5 paralog genes shared some functions during lung morphogenesis, Hoxa5 playing a predominant role.

The genetic and epigenetic journey of embryonic stem cells into mature neural cells

Epigenetic changes occur throughout life from embryonic development into adulthood. This results in the timely expression of developmentally important genes, determining the morphology and identity of different cell types and tissues within the body. Epigenetics regulate gene expression and cellular morphology through multiple mechanisms without alteration in the underlying DNA sequences. Different epigenetic mechanisms include chromatin condensation, post-translational modification of histone proteins, DNA cytosine marks, and the activity of non-coding RNA molecules. Epigenetics play key roles in development, stem cell differentiation, and have high impact in human disease. In this review, we will discuss our current knowledge about these epigenetic mechanisms, with a focus on histone and DNA marks. We will then talk about the genetics and epigenetics of embryonic stem cell self-renewal and differentiation into neural stem cells, and further into specific neuronal cell types.

A genetic approach to the transcriptional regulation of Hox gene clusters

The evolution of vertebrate genomes was accompanied by an astounding increase in the complexity of their regulatory modalities. Genetic redundancy resulting from large-scale genome duplications at the base of the chordate tree was repeatedly exploited by the functional redeployment of paralogous genes via innovations in their regulatory circuits. As a paradigm of such regulatory evolution, we have extensively studied those control mechanisms at work in-cis over vertebrate Hox gene clusters. Here, we review the portfolio of genetic strategies that have been developed to tackle the intricate relationship between genomic topography and the transcriptional activities in this gene family, and we describe some of the mechanistic insights we gained by using the HoxD cluster as an example. We discuss the high heuristic value of this system in our general understanding of how changes in transcriptional regulation can diversify gene function and thereby fuel morphological evolution.

Ancient Pbx-Hox signatures define hundreds of vertebrate developmental enhancers

BACKGROUND

Gene regulation through cis-regulatory elements plays a crucial role in development and disease. A major aim of the post-genomic era is to be able to read the function of cis-regulatory elements through scrutiny of their DNA sequence. Whilst comparative genomics approaches have identified thousands of putative regulatory elements, our knowledge of their mechanism of action is poor and very little progress has been made in systematically de-coding them.

RESULTS

Here, we identify ancient functional signatures within vertebrate conserved non-coding elements (CNEs) through a combination of phylogenetic footprinting and functional assay, using genomic sequence from the sea lamprey as a reference. We uncover a striking enrichment within vertebrate CNEs for conserved binding-site motifs of the Pbx-Hox hetero-dimer. We further show that these predict reporter gene expression in a segment specific manner in the hindbrain and pharyngeal arches during zebrafish development.

CONCLUSIONS

These findings evoke an evolutionary scenario in which many CNEs evolved early in the vertebrate lineage to co-ordinate Hox-dependent gene-regulatory interactions that pattern the vertebrate head. In a broader context, our evolutionary analyses reveal that CNEs are composed of tightly linked transcription-factor binding-sites (TFBSs), which can be systematically identified through phylogenetic footprinting approaches. By placing a large number of ancient vertebrate CNEs into a developmental context, our findings promise to have a significant impact on efforts toward de-coding gene-regulatory elements that underlie vertebrate development, and will facilitate building general models of regulatory element evolution.

Identification of novel Hoxa1 downstream targets regulating hindbrain, neural crest and inner ear development

Hox genes play a crucial role during embryonic patterning and organogenesis. Of the 39 Hox genes, Hoxa1 is the first to be expressed during embryogenesis and the only anterior Hox gene linked to a human syndrome. Hoxa1 is necessary for the proper development of the brainstem, inner ear and heart in humans and mice; however, almost nothing is known about the molecular downstream targets through which it exerts its function. To gain insight into the transcriptional network regulated by this protein, we performed microarray analysis on tissue microdissected from the prospective rhombomere 3-5 region of Hoxa1 null and wild type embryos. Due to the very early and transient expression of this gene, dissections were performed on early somite stage embryos during an eight-hour time window of development. Our array yielded a list of around 300 genes differentially expressed between the two samples. Many of the identified genes play a role in a specific developmental or cellular process. Some of the validated targets regulate early neural crest induction and specification. Interestingly, three of these genes, Zic1, Hnf1b and Foxd3, were down-regulated in the posterior hindbrain, where cardiac neural crest cells arise, which pattern the outflow tract of the heart. Other targets are necessary for early inner ear development, e.g. Pax8 and Fgfr3 or are expressed in specific hindbrain neurons regulating respiration, e.g. Lhx5. These findings allow us to propose a model where Hoxa1 acts in a genetic cascade upstream of genes controlling specific aspects of embryonic development, thereby providing insight into possible mechanisms underlying the human HoxA1-syndrome.

Reshuffling genomic landscapes to study the regulatory evolution of Hox gene clusters

The emergence of Vertebrata was accompanied by two rounds of whole-genome duplications. This enabled paralogous genes to acquire novel functions with high evolutionary potential, a process suggested to occur mostly by changes in gene regulation, rather than in protein sequences. In the case of Hox gene clusters, such duplications favored the appearance of distinct global regulations. To assess the impact of such "regulatory evolution" upon neo-functionalization, we developed PANTHERE (PAN-genomic Translocation for Heterologous Enhancer RE-shuffling) to bring the entire megabase-scale HoxD regulatory landscape in front of the HoxC gene cluster via a targeted translocation in vivo. At this chimeric locus, Hoxc genes could both interpret this foreign regulation and functionally substitute for their Hoxd counterparts. Our results emphasize the importance of evolving regulatory modules rather than their target genes in the process of neo-functionalization and offer a genetic tool to study the complexity of the vertebrate regulatory genome.

Directed neural differentiation of mouse embryonic stem cells is a sensitive system for the identification of novel Hox gene effectors

The evolutionarily conserved Hox family of homeodomain transcription factors plays fundamental roles in regulating cell specification along the anterior posterior axis during development of all bilaterian animals by controlling cell fate choices in a highly localized, extracellular signal and cell context dependent manner. Some studies have established downstream target genes in specific systems but their identification is insufficient to explain either the ability of Hox genes to direct homeotic transformations or the breadth of their patterning potential. To begin delineating Hox gene function in neural development we used a mouse ES cell based system that combines efficient neural differentiation with inducible Hoxb1 expression. Gene expression profiling suggested that Hoxb1 acted as both activator and repressor in the short term but predominantly as a repressor in the long run. Activated and repressed genes segregated in distinct processes suggesting that, in the context examined, Hoxb1 blocked differentiation while activating genes related to early developmental processes, wnt and cell surface receptor linked signal transduction and cell-to-cell communication. To further elucidate aspects of Hoxb1 function we used loss and gain of function approaches in the mouse and chick embryos. We show that Hoxb1 acts as an activator to establish the full expression domain of CRABPI and II in rhombomere 4 and as a repressor to restrict expression of Lhx5 and Lhx9. Thus the Hoxb1 patterning activity includes the regulation of the cellular response to retinoic acid and the delay of the expression of genes that commit cells to neural differentiation. The results of this study show that ES neural differentiation and inducible Hox gene expression can be used as a sensitive model system to systematically identify Hox novel target genes, delineate their interactions with signaling pathways in dictating cell fate and define the extent of functional overlap among different Hox genes.

Epigenomic reorganization of the clustered Hox genes in embryonic stem cells induced by retinoic acid

Retinoic acid (RA) regulates clustered Hox gene expression during embryogenesis and is required to establish the anterior-posterior body plan. Using mutant embryonic stem cell lines deficient in the RA receptor γ (RARγ) or Hoxa1 3'-RA-responsive element, we studied the kinetics of transcriptional and epigenomic patterning responses to RA. RARγ is essential for RA-induced Hox transcriptional activation, and deletion of its binding site in the Hoxa1 enhancer attenuates transcriptional and epigenomic activation of both Hoxa and Hoxb gene clusters. The kinetics of epigenomic reorganization demonstrate that complete erasure of the polycomb repressive mark H3K27me3 is not necessary to initiate Hox transcription. RARγ is not required to establish the bivalent character of Hox clusters, but RA/RARγ signaling is necessary to erase H3K27me3 from activated Hox genes during embryonic stem cell differentiation. Highly coordinated, long range epigenetic Hox cluster reorganization is closely linked to transcriptional activation and is triggered by RARγ located at the Hoxa1 3'-RA-responsive element.

An evolving NGF-Hoxd1 signaling pathway mediates development of divergent neural circuits in vertebrates

Species are endowed with unique sensory capabilities encoded by divergent neural circuits. One potential explanation for how divergent circuits have evolved is that conserved extrinsic signals are differentially interpreted by developing neurons of different species to yield unique patterns of axonal connections. Although NGF controls survival, maturation and axonal projections of nociceptors of different vertebrates, whether the NGF signal is differentially transduced in different species to yield unique features of nociceptor circuits is unclear. We identified a species-specific signaling module induced by NGF and mediated by a rapidly evolving Hox transcription factor, Hoxd1. Mice lacking Hoxd1 display altered nociceptor circuitry which resembles that normally found in chicks. Conversely, ectopic expression of Hoxd1 in developing chick nociceptors promotes a pattern of axonal projections reminiscent of the mouse. We propose that conserved growth factors control divergent neuronal transcriptional events which mediate interspecies differences in neural circuits and the behaviors they control.

Genome digging: insight into the mitochondrial genome of Homo

BACKGROUND

A fraction of the Neanderthal mitochondrial genome sequence has a similarity with a 5,839-bp nuclear DNA sequence of mitochondrial origin (numt) on the human chromosome 1. This fact has never been interpreted. Although this phenomenon may be attributed to contamination and mosaic assembly of Neanderthal mtDNA from short sequencing reads, we explain the mysterious similarity by integration of this numt (mtAncestor-1) into the nuclear genome of the common ancestor of Neanderthals and modern humans not long before their reproductive split.

PRINCIPAL FINDINGS

Exploiting bioinformatics, we uncovered an additional numt (mtAncestor-2) with a high similarity to the Neanderthal mtDNA and indicated that both numts represent almost identical replicas of the mtDNA sequences ancestral to the mitochondrial genomes of Neanderthals and modern humans. In the proteins, encoded by mtDNA, the majority of amino acids distinguishing chimpanzees from humans and Neanderthals were acquired by the ancestral hominins. The overall rate of nonsynonymous evolution in Neanderthal mitochondrial protein-coding genes is not higher than in other lineages. The model incorporating the ancestral hominin mtDNA sequences estimates the average divergence age of the mtDNAs of Neanderthals and modern humans to be 450,000-485,000 years. The mtAncestor-1 and mtAncestor-2 sequences were incorporated into the nuclear genome approximately 620,000 years and 2,885,000 years ago, respectively.

CONCLUSIONS

This study provides the first insight into the evolution of the mitochondrial DNA in hominins ancestral to Neanderthals and humans. We hypothesize that mtAncestor-1 and mtAncestor-2 are likely to be molecular fossils of the mtDNAs of Homo heidelbergensis and a stem Homo lineage. The d(N)/d(S) dynamics suggests that the effective population size of extinct hominins was low. However, the hominin lineage ancestral to humans, Neanderthals and H. heidelbergensis, had a larger effective population size and possessed genetic diversity comparable with those of chimpanzee and gorilla.

Integration of RNA processing and expression level control modulates the function of the Drosophila Hox gene Ultrabithorax during adult development

Although most metazoan genes undergo alternative splicing, the functional relevance of the majority of alternative splicing products is still unknown. Here we explore this problem in the Drosophila Hox gene Ultrabithorax (Ubx). Ubx produces a family of six protein isoforms through alternative splicing. To investigate the functional specificity of the Ubx isoforms, we studied their role during the formation of the Drosophila halteres, small dorsal appendages that are essential for normal flight. Our work shows that isoform Ia, which is encoded by all Ubx exons, is more efficient than isoform IVa, which lacks the amino acids coded by two small exons, in controlling haltere development and regulating Ubx downstream targets. However, our experiments also demonstrate that the functional differences among the Ubx isoforms can be compensated for by increasing the expression levels of the less efficient form. The analysis of the DNA-binding profiles of Ubx isoforms to a natural Ubx target, spalt, shows no major differences in isoform DNA-binding activities, suggesting that alternative splicing might primarily affect the regulatory capacity of the isoforms rather than their DNA-binding patterns. Our results suggest that to obtain distinct functional outputs during normal development genes must integrate the generation of qualitative differences by alternative splicing to quantitative processes affecting isoform protein expression levels.

The C. elegans peroxidasin PXN-2 is essential for embryonic morphogenesis and inhibits adult axon regeneration

Peroxidasins form a highly conserved family of extracellular peroxidases of unknown cellular function. We identified the C. elegans peroxidasin PXN-2 in screens for mutants defective in embryonic morphogenesis. We find that PXN-2 is essential for specific stages of embryonic morphogenesis and muscle-epidermal attachment, and is also required postembryonically for basement membrane integrity. The peroxidase catalytic activity of PXN-2 is necessary for these developmental roles. pxn-2 mutants display aberrant ultrastructure of the extracellular matrix, suggesting a role in basement membrane consolidation. PXN-2 affects specific axon guidance choice points in the developing nervous system but is dispensable for maintenance of process positions. In adults, loss of pxn-2 function promotes regrowth of axons after injury, providing the first evidence that C. elegans extracellular matrix can play an inhibitory role in axon regeneration. Loss of function in the closely related C. elegans peroxidasin pxn-1 does not cause overt developmental defects. Unexpectedly, pxn-2 mutant phenotypes are suppressed by loss of function in pxn-1 and exacerbated by overexpression of wild-type pxn-1, indicating that PXN-1 and PXN-2 have antagonistic functions. These results demonstrate that peroxidasins play crucial roles in development and reveal a new role for peroxidasins as extracellular inhibitors of axonal regeneration.

Duplicate dmbx1 genes regulate progenitor cell cycle and differentiation during zebrafish midbrain and retinal development

Background

The Dmbx1 gene is important for the development of the midbrain and hindbrain, and mouse gene targeting experiments reveal that this gene is required for mediating postnatal and adult feeding behaviours. A single Dmbx1 gene exists in terrestrial vertebrate genomes, while teleost genomes have at least two paralogs. We compared the loss of function of the zebrafish dmbx1a and dmbx1b genes in order to gain insight into the molecular mechanism by which dmbx1 regulates neurogenesis, and to begin to understand why these duplicate genes have been retained in the zebrafish genome.

Results

Using gene knockdown experiments we examined the function of the dmbx1 gene paralogs in zebrafish, dmbx1a and dmbx1b in regulating neurogenesis in the developing retina and midbrain. Dose-dependent loss of dmbx1a and dmbx1b function causes a significant reduction in growth of the midbrain and retina that is evident between 48-72 hpf. We show that this phenotype is not due to patterning defects or persistent cell death, but rather a deficit in progenitor cell cycle exit and differentiation. Analyses of the morphant retina or anterior hindbrain indicate that paralogous function is partially diverged since loss of dmbx1a is more severe than loss of dmbx1b. Molecular evolutionary analyses of the Dmbx1 genes suggest that while this gene family is conservative in its evolution, there was a dramatic change in selective constraint after the duplication event that gave rise to the dmbx1a and dmbx1b gene families in teleost fish, suggestive of positive selection. Interestingly, in contrast to zebrafish dmbx1a, over expression of the mouse Dmbx1 gene does not functionally compensate for the zebrafish dmbx1a knockdown phenotype, while over expression of the dmbx1b gene only partially compensates for the dmbx1a knockdown phenotype.

Conclusion

Our data suggest that both zebrafish dmbx1a and dmbx1b genes are retained in the fish genome due to their requirement during midbrain and retinal neurogenesis, although their function is partially diverged. At the cellular level, Dmbx1 regulates cell cycle exit and differentiation of progenitor cells. The unexpected observation of putative post-duplication positive selection of teleost Dmbx1 genes, especially dmbx1a, and the differences in functionality between the mouse and zebrafish genes suggests that the teleost Dmbx1 genes may have evolved a diverged function in the regulation of neurogenesis.

Uner tan syndrome: history, clinical evaluations, genetics, and the dynamics of human quadrupedalism

This review includes for the first time a dynamical systems analysis of human quadrupedalism in Uner Tan syndrome, which is characterized by habitual quadrupedalism, impaired intelligence, and rudimentary speech. The first family was discovered in a small village near Iskenderun, and families were later found in Adana and two other small villages near Gaziantep and Canakkale. In all the affected individuals dynamic balance was impaired during upright walking, and they habitually preferred walking on all four extremities. MRI scans showed inferior cerebellovermian hypoplasia with slightly simplified cerebral gyri in three of the families, but appeared normal in the fourth. PET scans showed a decreased glucose metabolic activity in the cerebellum, vermis and, to a lesser extent the cerebral cortex, except for one patient, whose MRI scan also appeared to be normal. All four families had consanguineous marriages in their pedigrees, suggesting autosomal recessive transmission. The syndrome was genetically heterogeneous. Since the initial discoveries more cases have been found, and these exhibit facultative quadrupedal locomotion, and in one case, late childhood onset. It has been suggested that the human quadrupedalism may, at least, be a phenotypic example of reverse evolution. From the viewpoint of dynamic systems theory, it was concluded there may not be a single factor that predetermines human quadrupedalism in Uner Tan syndrome, but that it may involve self-organization, brain plasticity, and rewiring, from the many decentralized and local interactions among neuronal, genetic, and environmental subsystems.

Improving Hox protein classification across the major model organisms

The family of Hox-proteins has been a major focus of research for over 30 years. Hox-proteins are crucial to the correct development of bilateral organisms, however, some uncertainty remains as to which Hox-proteins are functionally equivalent across different species. Initial classification of Hox-proteins was based on phylogenetic analysis of the 60 amino acid homeodomain. This approach was successful in classifying Hox-proteins with differing homeodomains, but the relationships of Hox-proteins with nearly identical homeodomains, yet distinct biological functions, could not be resolved. Correspondingly, these ‘problematic’ proteins were classified into one large unresolved group. Other classifications used the relative location of the Hox-protein coding genes on the chromosome (synteny) to further resolve this group. Although widely used, this synteny-based classification is inconsistent with experimental evidence from functional equivalence studies. These inconsistencies led us to re-examine and derive a new classification for the Hox-protein family using all Hox-protein sequences available in the GenBank non-redundant protein database (NCBI-nr). We compare the use of the homeodomain, the homeodomain with conserved flanking regions (the YPWM and linker region), and full length Hox-protein sequences as a basis for classification of Hox-proteins. In contrast to previous attempts, our approach is able to resolve the relationships for the ‘problematic’ as well as ABD-B-like Hox-proteins. We highlight differences to previous classifications and clarify the relationships of Hox-proteins across the five major model organisms, Caenorhabditis elegans, Drosophila melanogaster, Branchiostoma floridae, Mus musculus and Danio rerio. Comparative and functional analysis of Hox-proteins, two fields crucial to understanding the development of bilateral organisms, have been hampered by difficulties in predicting functionally equivalent Hox-proteins across species. Our classification scheme offers a higher-resolution classification that is in accordance with phylogenetic as well as experimental data and, thereby, provides a novel basis for experiments, such as comparative and functional analyses of Hox-proteins.

Mouse and zebrafish Hoxa3 orthologues have nonequivalent in vivo protein function

Hox genes play evolutionarily conserved roles in specifying axial position during embryogenesis. A prevailing paradigm is that changes in Hox gene expression drive evolution of metazoan body plans. Conservation of Hox function across species, and among paralogous Hox genes within a species, supports a model of functional equivalence. In this report, we demonstrate that zebrafish hoxa3a (zfhoxa3a) expressed from the mouse Hoxa3 locus can substitute for mouse Hoxa3 in some tissues, but has distinct or null phenotypes in others. We further show, by using an allele encoding a chimeric protein, that this difference maps primarily to the zfhoxa3a C-terminal domain. Our data imply that the mouse and zebrafish proteins have diverged considerably since their last common ancestor, and that the major difference between them resides in the C-terminal domain. Our data further show that Hox protein function can evolve independently in different cell types or for specific functions. The inability of zfhoxa3a to perform all of the normal roles of mouse Hoxa3 illustrates that Hox orthologues are not always functionally interchangeable.

Evolution of stress-regulated gene expression in duplicate genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >∼0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments.

Plants have developed a multitude of response mechanisms to survive stressful environments. Since the environment is highly variable, these stress response mechanisms are expected to undergo frequent innovation. Duplicate genes represent a potential source for such innovation. In this paper, we explored the evolutionary changes in stress responses at the transcriptional level among duplicated genes in the model plant Arabidopsis thaliana. We found that after gene duplication, ancestral stress responses tend to be retained by only one of the gene duplicates (partitioning). In addition, the pattern of partitioning of multiple stress responses is extremely asymmetric, where one duplicate tends to inherit most or all of the ancestral stress responses. We present evidence that the asymmetric loss of stress responses is correlated with the asymmetric loss of putative transcription factor binding sites. Interestingly, those duplicate genes inheriting few or no ancestral responses tend to have gained new stress responses, providing support for the model that gene duplicates are a source of innovation. Our findings provide important insight into the mechanisms of gene function evolution and lay the foundation for experimental studies to determine the significance of gain of stress responses in plant adaptation.