Homeotic genes and the evolution of arthropods and chordates

Evolutionarily new genes in humans with disease phenotypes reveal functional enrichment patterns shaped by adaptive innovation and sexual selection

New genes (or young genes) are genetic novelties pivotal in mammalian evolution. However, their phenotypic impacts and evolutionary patterns over time remain elusive in humans due to the technical and ethical complexities of functional studies. Integrating gene age dating with Mendelian disease phenotyping, our research shows a gradual rise in disease gene proportion as gene age increases. Logistic regression modeling indicates that this increase in older genes may be related to their longer sequence lengths and higher burdens of deleterious de novo germline variants (DNVs). We also find a steady integration of new genes with biomedical phenotypes into the human genome over macroevolutionary timescales (~0.07% per million years). Despite this stable pace, we observe distinct patterns in phenotypic enrichment, pleiotropy, and selective pressures across gene ages. Notably, young genes show significant enrichment in diseases related to the male reproductive system, indicating strong sexual selection. Young genes also exhibit disease-related functions in tissues and systems potentially linked to human phenotypic innovations, such as increased brain size, musculoskeletal phenotypes, and color vision. We further reveal a logistic growth pattern of pleiotropy over evolutionary time, indicating a diminishing marginal growth of new functions for older genes due to intensifying selective constraints over time. We propose a pleiotropy-barrier model that delineates higher potentials for phenotypic innovation in young genes compared to older genes, a process that is subject to natural selection. Our study demonstrates that evolutionarily new genes are critical in influencing human reproductive evolution and adaptive phenotypic innovations driven by sexual and natural selection, with low pleiotropy as a selective advantage.

The role of heterochronic gene expression and regulatory architecture in early developmental divergence

New developmental programs can evolve through adaptive changes to gene expression. The annelid Streblospio benedicti has a developmental dimorphism, which provides a unique intraspecific framework for understanding the earliest genetic changes that take place during developmental divergence. Using comparative RNAseq through ontogeny, we find that only a small proportion of genes are differentially expressed at any time, despite major differences in larval development and life history. These genes shift expression profiles across morphs by either turning off any expression in one morph or changing the timing or amount of gene expression. We directly connect the contributions of these mechanisms to differences in developmental processes. We examine F1 offspring - using reciprocal crosses - to determine maternal mRNA inheritance and the regulatory architecture of gene expression. These results highlight the importance of both novel gene expression and heterochronic shifts in developmental evolution, as well as the trans-acting regulatory factors in initiating divergence.

Body part-specific development in termite caste differentiation: crosstalk between hormonal actions and developmental toolkit genes

In social insects, interactions among colony members trigger caste differentiation with morphological modifications. During caste differentiation in termites, body parts and caste-specific morphologies are modified during postembryonic development under endocrine controls such as juvenile hormone (JH) and ecdysone. In addition to endocrine factors, developmental toolkit genes such as Hox- and appendage-patterning genes also contribute to the caste-specific body part modifications. These toolkits are thought to provide spatial information for specific morphogenesis. During social evolution, the complex crosstalks between physiological and developmental mechanisms should be established, leading to the sophisticated caste systems. This article reviews recent studies on these mechanisms underlying the termite caste differentiation and addresses implications for the evolution of caste systems in termites.

Adaptive tail-length evolution in deer mice is associated with differential Hoxd13 expression in early development

Variation in the size and number of axial segments underlies much of the diversity in animal body plans. Here we investigate the evolutionary, genetic and developmental mechanisms driving tail-length differences between forest and prairie ecotypes of deer mice (Peromyscus maniculatus). We first show that long-tailed forest mice perform better in an arboreal locomotion assay, consistent with tails being important for balance during climbing. We then identify six genomic regions that contribute to differences in tail length, three of which associate with caudal vertebra length and the other three with vertebra number. For all six loci, the forest allele increases tail length, indicative of the cumulative effect of natural selection. Two of the genomic regions associated with variation in vertebra number contain Hox gene clusters. Of those, we find an allele-specific decrease in Hoxd13 expression in the embryonic tail bud of long-tailed forest mice, consistent with its role in axial elongation. Additionally, we find that forest embryos have more presomitic mesoderm than prairie embryos and that this correlates with an increase in the number of neuromesodermal progenitors, which are modulated by Hox13 paralogues. Together, these results suggest a role for Hoxd13 in the development of natural variation in adaptive morphology on a microevolutionary timescale.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

Homeodomain complex formation and biomolecular condensates in Hox gene regulation

Hox genes are a family of homeodomain transcription factors that regulate specialized morphological structures along the anterior-posterior axis of metazoans. Over the past few decades, researchers have focused on defining how Hox factors with similar in vitro DNA binding activities achieve sufficient target specificity to regulate distinct cell fates in vivo. In this review, we highlight how protein interactions with other transcription factors, many of which are also homeodomain proteins, result in the formation of transcription factor complexes with enhanced DNA binding specificity. These findings suggest that Hox-regulated enhancers utilize distinct combinations of homeodomain binding sites, many of which are low-affinity, to recruit specific Hox complexes. However, low-affinity sites can only yield reproducible responses with high transcription factor concentrations. To overcome this limitation, recent studies revealed how transcription factors, including Hox factors, use intrinsically disordered domains (IDRs) to form biomolecular condensates that increase protein concentrations. Moreover, Hox factors with altered IDRs have been associated with altered transcriptional activity and human disease states, demonstrating the importance of IDRs in mediating essential Hox output. Collectively, these studies highlight how Hox factors use their DNA binding domains, protein-protein interaction domains, and IDRs to form specific transcription factor complexes that yield accurate gene expression.

CRISPR/Cas9-Mediated Mutagenesis of Antennapedia in Spodoptera frugiperda

The homeotic gene Antennapedia (Antp) has been identified as playing a pivotal role in the morphogenesis of the thorax and wings across various insect species. Leveraging insights from previous studies, the functional characterization of Antp in S. frugiperda was undertaken using RT-qPCR and the CRISPR/Cas9 genome-editing system. Phylogenetic analyses indicate that Antp shares a high degree of sequence homology among Lepidoptera species. The expression profile of SfAntp was detected by RT-qPCR. The results showed that SfAntp was expressed in the whole growth cycle of S. frugiperda, the expression level was the highest in the egg stage, and the expression level was higher from 12 h to 48 h. Tissue-specific expression profiling demonstrated that SfAntp was most abundantly expressed in the thoracic segments and legs. To functionally disrupt SfAntp, two sgRNA sites were designed at the first exon of SfAntp and the gene was knocked out by CRISPR/Cas9 via microinjection. The results showed that the deletion of SfAntp produced a mutant phenotype of thoracic fusion, thoracic leg defect, leg-like protrusions between the head and thoracic segments and pupation deformity. In addition, deletion of SfAntp resulted in high embryo mortality. Through DNA sequencing, it was found that the target site of the SfAntp mutant had different degrees of frameshift mutations, indicating that the mutant phenotype was indeed caused by the knockout of SfAntp.

Circular and unified analysis in network neuroscience

Genuinely new discovery transcends existing knowledge. Despite this, many analyses in systems neuroscience neglect to test new speculative hypotheses against benchmark empirical facts. Some of these analyses inadvertently use circular reasoning to present existing knowledge as new discovery. Here, I discuss that this problem can confound key results and estimate that it has affected more than three thousand studies in network neuroscience over the last decade. I suggest that future studies can reduce this problem by limiting the use of speculative evidence, integrating existing knowledge into benchmark models, and rigorously testing proposed discoveries against these models. I conclude with a summary of practical challenges and recommendations.

Integrating Complex Life Cycles in Comparative Developmental Biology

The goal of comparative developmental biology is identifying mechanistic differences in embryonic development between different taxa and how these evolutionary changes have led to morphological and organizational differences in adult body plans. Much of this work has focused on direct-developing species in which the adult forms straight from the embryo and embryonic modifications have direct effects on the adult. However, most animal lineages are defined by indirect development, in which the embryo gives rise to a larval body plan and the adult forms by transformation of the larva. Historically, much of our understanding of complex life cycles is viewed through the lenses of ecology and zoology. In this review, we discuss the importance of establishing developmental rather than morphological or ecological criteria for defining developmental mode and explicitly considering the evolutionary implications of incorporating complex life cycles into broad developmental comparisons of embryos across metazoans.

Evolutionarily new genes in humans with disease phenotypes reveal functional enrichment patterns shaped by adaptive innovation and sexual selection

New genes (or young genes) are structural novelties pivotal in mammalian evolution. Their phenotypic impact on humans, however, remains elusive due to the technical and ethical complexities in functional studies. Through combining gene age dating with Mendelian disease phenotyping, our research reveals that new genes associated with disease phenotypes steadily integrate into the human genome at a rate of ~ 0.07% every million years over macroevolutionary timescales. Despite this stable pace, we observe distinct patterns in phenotypic enrichment, pleiotropy, and selective pressures between young and old genes. Notably, young genes show significant enrichment in the male reproductive system, indicating strong sexual selection. Young genes also exhibit functions in tissues and systems potentially linked to human phenotypic innovations, such as increased brain size, bipedal locomotion, and color vision. Our findings further reveal increasing levels of pleiotropy over evolutionary time, which accompanies stronger selective constraints. We propose a "pleiotropy-barrier" model that delineates different potentials for phenotypic innovation between young and older genes subject to natural selection. Our study demonstrates that evolutionary new genes are critical in influencing human reproductive evolution and adaptive phenotypic innovations driven by sexual and natural selection, with low pleiotropy as a selective advantage.

Species differences in hormonally mediated gene expression underlie the evolutionary loss of sexually dimorphic coloration in Sceloporus lizards

Phenotypic sexual dimorphism often involves the hormonal regulation of sex-biased expression for underlying genes. However, it is generally unknown whether the evolution of hormonally mediated sexual dimorphism occurs through upstream changes in tissue sensitivity to hormone signals, downstream changes in responsiveness of target genes, or both. Here, we use comparative transcriptomics to explore these possibilities in 2 species of Sceloporus lizards exhibiting different patterns of sexual dichromatism. Sexually dimorphic S. undulatus develops blue and black ventral coloration in response to testosterone, while sexually monomorphic S. virgatus does not, despite exhibiting similar sex differences in circulating testosterone levels. We administered testosterone implants to juveniles of each species and used RNAseq to quantify gene expression in ventral skin. Transcriptome-wide responses to testosterone were stronger in S. undulatus than in S. virgatus, suggesting species differences in tissue sensitivity to this hormone signal. Species differences in the expression of genes for androgen metabolism and sex hormone-binding globulin were consistent with this idea, but expression of the androgen receptor gene was higher in S. virgatus, complicating this interpretation. Downstream of androgen signaling, we found clear species differences in hormonal responsiveness of genes related to melanin synthesis, which were upregulated by testosterone in S. undulatus, but not in S. virgatus. Collectively, our results indicate that hormonal regulation of melanin synthesis pathways contributes to the development of sexual dimorphism in S. undulatus, and that changes in the hormonal responsiveness of these genes in S. virgatus contribute to the evolutionary loss of ventral coloration.

A surrogate forelimb: Evolution, function and development of the avian cervical spine

The neck is a critical portion of the avian spine, one that works in tandem with the beak to act as a surrogate forelimb and allows birds to manipulate their surroundings despite the lack of a grasping capable hand. Birds display an incredible amount of diversity in neck morphology across multiple anatomical scales-from varying cervical counts down to intricate adaptations of individual vertebrae. Despite this morphofunctional disparity, little is known about the drivers of this enormous variation, nor how neck evolution has shaped avian macroevolution. To promote interest in this system, I review the development, function and evolution of the avian cervical spine. The musculoskeletal anatomy, basic kinematics and development of the avian neck are all documented, but focus primarily upon commercially available taxa. In addition, recent work has quantified the drivers of extant morphological variation across the avian neck, as well as patterns of integration between the neck and other skeletal elements. However, the evolutionary history of the avian cervical spine, and its contribution to the diversification and success of modern birds is currently unknown. Future work should aim to broaden our understanding of the cervical anatomy, development and kinematics to include a more diverse selection of extant birds, while also considering the macroevolutionary drivers and consequences of this important section of the avian spine.

A Micro-evolutionary Change in Target Binding Sites as a Key Determinant of Ultrabithorax Function in Drosophila

Hox genes encode Homeodomain-containing transcription factors, which specify segmental identities along the anterior-posterior axis. Functional changes in Hox genes have been directly implicated in the evolution of body plans across the metazoan lineage. The Hox protein Ultrabithorax (Ubx) is expressed and required in developing third thoracic (T3) segments in holometabolous insects studied so far, particularly, of the order Coleoptera, Lepidoptera and Diptera. Ubx function is key to specify differential development of the second (T2) and T3 thoracic segments in these insects. While Ubx is expressed in the third thoracic segment in developing larvae of Hymenopteran Apis mellifera, the morphological differences between T2 and T3 are subtle. To identify evolutionary changes that are behind the differential function of Ubx in Drosophila and Apis, which are diverged for more than 350 million years, we performed comparative analyses of genome wide Ubx-binding sites between these two insects. Our studies reveal that a motif with a TAAAT core is a preferred binding site for Ubx in Drosophila, but not in Apis. Biochemical and transgenic assays suggest that in Drosophila, the TAAAT core sequence in the Ubx binding sites is required for Ubx-mediated regulation of two of its target genes studied here; CG13222, a gene that is normally upregulated by Ubx and vestigial (vg), whose expression is repressed by Ubx in T3. Interestingly, changing the TAAT site to a TAAAT site was sufficient to bring an otherwise unresponsive enhancer of the vg gene from Apis under the control of Ubx in a Drosophila transgenic assay. Taken together, our results suggest an evolutionary mechanism by which critical wing patterning genes might have come under the regulation of Ubx in the Dipteran lineage.

Linking Vertebrate Gene Duplications to the New Head Hypothesis

Vertebrates have diverse morphologies and various anatomical novelties that set them apart from their closest invertebrate relatives. A conspicuous head housing a large brain, paired sense organs, and protected by a skeleton of cartilage and bone is unique to vertebrates and is a defining feature of this taxon. Gans and Northcutt (1980s) proposed that the evolution of this "new head" was dependent on two key developmental innovations: neural crest cells (NCCs) and ectodermal placodes. NCCs are migratory embryonic cells that form bone, cartilage, and neurons in the new head. Based on genome size, Ohno (1970s) proposed a separate hypothesis, stating that vertebrate genome content was quadrupled via two rounds (2R) of whole genome duplications (WGDs), and the surplus of genetic material potentiated vertebrate morphological diversification. While both hypotheses offer explanations for vertebrate success, it is unclear if, and how, the "new head" and "2R" hypotheses are linked. Here, we consider both hypotheses and evaluate the experimental evidence connecting the two. Overall, evidence suggests that while the origin of the NC GRN predates the vertebrate WGDs, these genomic events may have potentiated the evolution of distinct genetic subnetworks in different neural crest subpopulations. We describe the general composition of the NC GRN and posit that its increased developmental modularity facilitated the independent evolution of NC derivatives and the diversification of the vertebrate head skeleton. Lastly, we discuss experimental strategies needed to test whether gene duplications drove the diversification of neural crest derivatives and the "new head".

Homeobox protein A1-like and DNA methylation regulate embryo-specific Zinc finger protein 615 gene expression and embryonic development in the silkworm Bombyx mori

DNA methylation and transcription factors play roles in gene expression and animal development. In insects, DNA methylation modifies gene bodies, but how DNA methylation and transcription factors regulate gene expression is unclear. In this study, we investigated the mechanism that regulates the expression of Bombyx mori Zinc finger protein 615 (ZnF 615), which is a downstream gene of DNA methyltransferase 1 (Dnmt1), and its effects on the regulation of embryonic development. By progressively truncating the ZnF 615 promoter, it was found that the -223 and -190 nt region, which contains homeobox (Hox) protein cis-regulatory elements (CREs), had the greatest impact on the transcription of ZnF 615. RNA interference (RNAi)-mediated knockdown and overexpression of Hox family genes showed that Hox A1-like can enhance the messenger RNA level of ZnF 615. Further studies showed that Hox A1-like regulates ZnF 615 expression by directly binding to the -223 and -190 nt region of its promoter. Simultaneous RNAi-mediated knockdown or overexpression of Hox A1-like and Dnmt1 significantly inhibited or enhanced the regulatory effect of either gene alone on ZnF 615 expression, suggesting that both DNA methylation of gene bodies and binding of transcription factors to promoters are essential for gene expression. RNAi-mediated knockdown of Hox A1-like and Dnmt1 showed that the embryonic development was retarded and the hatching rate was decreased. Taken together, these data suggest that Hox A1-like and DNA methylation enhance the expression of ZnF 615, thereby affecting the development of B. mori embryos.

Upregulation of Hox genes leading to caste-specific morphogenesis in a termite

BACKGROUND

In social insects, interactions among colony members trigger caste differentiation with morphological modifications. In termite caste differentiation, caste-specific morphologies (such as mandibles in soldiers, genital organs in reproductives or wings in alates) are well developed during post-embryonic development under endocrine controls (e.g., juvenile hormone and ecdysone). Since body part-specific morphogenesis in caste differentiation is hormonally regulated by global factors circulated throughout the body, positional information should be required for the caste-specific and also body part-specific morphogenesis. To identify factors providing the positional information, expression and functional analyses of eight Hox genes were carried out during the three types of caste differentiation (i.e., soldier, neotenic and alate differentiation) in a termite, Hodotermopsis sjostedti.

RESULTS

Spatio-temporal patterns of Hox gene expression during caste differentiation were elucidated by real-time qPCR, showing the caste-specific upregulations of Hox genes during the differentiation processes. Among eight Hox genes, Deformed (Dfd) was upregulated specifically in mandibles in soldier differentiation, abdominal-A (abd-A) and Abdominal-B (Abd-B) were upregulated in the abdomen in neotenic differentiation, while Sex-comb reduced (Scr) and Antennapedia (Antp) were upregulated during alate differentiation. Furthermore, RNAi knockdown of Dfd in soldier differentiation and of abd-A and Abd-B in neotenic differentiation distorted the modifications of caste-specific morphologies.

CONCLUSIONS

Gene expression and functional analyses in this study revealed that, in the caste differentiation in termites, upregulation of Hox genes provide positional identities of body segments, resulting in the caste-specific morphogenesis. The acquisition of such developmental modifications would have enabled the evolution of sophisticated caste systems in termites.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

Carcino-Evo-Devo, A Theory of the Evolutionary Role of Hereditary Tumors

A theory of the evolutionary role of hereditary tumors, or the carcino-evo-devo theory, is being developed. The main hypothesis of the theory, the hypothesis of evolution by tumor neofunctionalization, posits that hereditary tumors provided additional cell masses during the evolution of multicellular organisms for the expression of evolutionarily novel genes. The carcino-evo-devo theory has formulated several nontrivial predictions that have been confirmed in the laboratory of the author. It also suggests several nontrivial explanations of biological phenomena previously unexplained by the existing theories or incompletely understood. By considering three major types of biological development-individual, evolutionary, and neoplastic development-within one theoretical framework, the carcino-evo-devo theory has the potential to become a unifying biological theory.

Gsx2, but not Gsx1, is necessary for early forebrain patterning and long-term survival in zebrafish

BACKGROUND

Homeobox transcription factor encoding genes, genomic screen homeobox 1 and 2 (gsx1 and gsx2), are expressed during neurodevelopment in multiple vertebrates. However, we have limited knowledge of the dynamic expression of these genes through developmental time and the gene networks that they regulate in zebrafish.

RESULTS

We confirmed that gsx1 is expressed initially in the hindbrain and diencephalon and later in the optic tectum, pretectum, and cerebellar plate. gsx2 is expressed in the early telencephalon and later in the pallium and olfactory bulb. gsx1 and gsx2 are co-expressed in the hypothalamus, preoptic area, and hindbrain, however, rarely co-localize in the same cells. gsx1 and gsx2 mutant zebrafish were made with TALENs. gsx1 mutants exhibit stunted growth, however, they survive to adulthood and are fertile. gsx2 mutants experience swim bladder inflation failure that prevents survival. We also observed significantly reduced expression of multiple forebrain patterning distal-less homeobox genes in mutants, and expression of foxp2 was not significantly affected.

CONCLUSIONS

This work provides novel tools with which other target genes and functions of Gsx1 and Gsx2 can be characterized across the central nervous system to better understand the unique and overlapping roles of these highly conserved transcription factors.

Single-cell atavism reveals an ancient mechanism of cell type diversification in a sea anemone

Cnidocytes are the explosive stinging cells unique to cnidarians (corals, jellyfish, etc). Specialized for prey capture and defense, cnidocytes comprise a group of over 30 morphologically and functionally distinct cell types. These unusual cells are iconic examples of biological novelty but the developmental mechanisms driving diversity of the stinging apparatus are poorly characterized, making it challenging to understand the evolutionary history of stinging cells. Using CRISPR/Cas9-mediated genome editing in the sea anemone Nematostella vectensis, we show that a single transcription factor (NvSox2) acts as a binary switch between two alternative stinging cell fates. Knockout of NvSox2 causes a transformation of piercing cells into ensnaring cells, which are common in other species of sea anemone but appear to have been silenced in N. vectensis. These results reveal an unusual case of single-cell atavism and expand our understanding of the diversification of cell type identity.

Early effects of gene duplication on the robustness and phenotypic variability of gene regulatory networks

Background

Research on gene duplication is abundant and comes from a wide range of approaches, from high-throughput analyses and experimental evolution to bioinformatics and theoretical models. Notwithstanding, a consensus is still lacking regarding evolutionary mechanisms involved in evolution through gene duplication as well as the conditions that affect them. We argue that a better understanding of evolution through gene duplication requires considering explicitly that genes do not act in isolation. It demands studying how the perturbation that gene duplication implies percolates through the web of gene interactions. Due to evolution’s contingent nature, the paths that lead to the final fate of duplicates must depend strongly on the early stages of gene duplication, before gene copies have accumulated distinctive changes.

Methods

Here we use a widely-known model of gene regulatory networks to study how gene duplication affects network behavior in early stages. Such networks comprise sets of genes that cross-regulate. They organize gene activity creating the gene expression patterns that give cells their phenotypic properties. We focus on how duplication affects two evolutionarily relevant properties of gene regulatory networks: mitigation of the effect of new mutations and access to new phenotypic variants through mutation.

Results

Among other observations, we find that those networks that are better at maintaining the original phenotype after duplication are usually also better at buffering the effect of single interaction mutations and that duplication tends to enhance further this ability. Moreover, the effect of mutations after duplication depends on both the kind of mutation and genes involved in it. We also found that those phenotypes that had easier access through mutation before duplication had higher chances of remaining accessible through new mutations after duplication.

Conclusion

Our results support that gene duplication often mitigates the impact of new mutations and that this effect is not merely due to changes in the number of genes. The work that we put forward helps to identify conditions under which gene duplication may enhance evolvability and robustness to mutations.

Singular value decomposition of protein sequences as a method to visualize sequence and residue space

Singular value decomposition (SVD) of multiple sequence alignments (MSAs) is an important and rigorous method to identify subgroups of sequences within the MSA, and to extract consensus and covariance sequence features that define the alignment and distinguish the subgroups. This information can be correlated to structure, function, stability, and taxonomy. However, the mathematics of SVD is unfamiliar to many in the field of protein science. Here, we attempt to present an intuitive yet comprehensive description of SVD analysis of MSAs. We begin by describing the underlying mathematics of SVD in a way that is both rigorous and accessible. Next, we use SVD to analyze sequences generated with a simplified model in which the extent of sequence conservation and covariance between different positions is controlled, to show how conservation and covariance produce features in the decomposed coordinate system. We then use SVD to analyze alignments of two protein families, the homeodomain and the Ras superfamilies. Both families show clear evidence of sequence clustering when projected into singular value space. We use k-means clustering to group MSA sequences into specific clusters, show how the residues that distinguish these clusters can be identified, and show how these clusters can be related to taxonomy and function. We end by providing a description a set of Python scripts that can be used for SVD analysis of MSAs, displaying results, and identifying and analyzing sequence clusters. These scripts are freely available on GitHub.

Evolution of stickleback spines through independent cis-regulatory changes at HOXDB

Understanding the mechanisms leading to new traits or additional features in organisms is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different fish genera to alter the length and number of the prominent dorsal spines used to classify stickleback species. In Gasterosteus aculeatus (typically 'three-spine sticklebacks'), a variant HOXDB allele is genetically linked to shortening an existing spine and adding an additional spine. In Apeltes quadracus (typically 'four-spine sticklebacks'), a variant HOXDB allele is associated with lengthening a spine and adding an additional spine in natural populations. The variant alleles alter the same non-coding enhancer region in the HOXDB locus but do so by diverse mechanisms, including single-nucleotide polymorphisms, deletions and transposable element insertions. The independent regulatory changes are linked to anterior expansion or contraction of HOXDB expression. We propose that associated changes in spine lengths and numbers are partial identity transformations in a repeating skeletal series that forms major defensive structures in fish. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes.

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

Characterizing Hox genes in mayflies (Ephemeroptera), with Hexagenia limbata as a new mayfly model

BACKGROUND

Hox genes are key regulators of appendage development in the insect body plan. The body plan of mayfly (Ephemeroptera) nymphs differs due to the presence of abdominal appendages called gills. Despite mayflies' phylogenetic position in Paleoptera and novel morphology amongst insects, little is known of their developmental genetics, such as the appendage-regulating Hox genes. To address this issue we present an annotated, early instar transcriptome and embryonic expression profiles for Antennapedia, Ultrabithorax, and Abdominal A proteins in the mayfly Hexagenia limbata, identify putative Hox protein sequences in the mayflies H. limbata, Cloeon dipterum, and Ephemera danica, and describe the genomic organization of the Hox gene cluster in E. danica.

RESULTS

Transcriptomic sequencing of early instar H. limbata nymphs yielded a high-quality assembly of 83,795 contigs, of which 22,975 were annotated against Folsomia candida, Nilaparvata lugens, Zootermopsis nevadensis and UniRef90 protein databases. Homeodomain protein phylogeny and peptide annotations identified coding sequences for eight of the ten canonical Hox genes (excluding zerknüllt/Hox3 and fushi tarazu) in H. limbata and C. dipterum, and all ten in E. danica. Mayfly Hox protein sequences and embryonic expression patterns of Antp, Ubx, and Abd-A appear highly conserved with those seen in other non-holometabolan insects. Similarly, the genomic organization of the Hox cluster in E. danica resembles that seen in most insects.

CONCLUSIONS

We present evidence that mayfly Hox peptide sequences and the embryonic expression patterns for Antp, Ubx, and Abd-A are extensively conserved with other insects, as is organization of the mayfly Hox gene cluster. The protein data suggest mayfly Antp, Ubx, and Abd-A play appendage promoting and repressing roles during embryogenesis in the thorax and abdomen, respectively, as in other insects. The identified expression of eight Hox genes, including Ubx and abd-A, in early instar nymphs further indicates a post-embryonic role, possibly in gill development. These data provide a basis for H. limbata as a complementary Ephemeridae model to the growing repertoire of mayfly model species and molecular techniques.

Hox genes collaborate with helix-loop-helix factor Grainyhead to promote neuroblast apoptosis along the anterior-posterior axis of the Drosophila larval central nervous system

Hox genes code for a family of a homeodomain (HD) containing transcription factors that use TALE-HD containing factors Pbx/Exd and Meis/Hth to specify the development of the anterior-posterior (AP) axis of an organism. However, the absence of TALE-HD containing factors from specific tissues emphasizes the need to identify and validate new Hox cofactors. In Drosophila central nervous system (CNS), Hox execute segment-specific apoptosis of neural stem cells (neuroblasts-NBs) and neurons. In abdominal segments of larval CNS, Hox gene Abdominal-A (AbdA) mediates NB apoptosis with the help of Exd and bHLH factor Grainyhead (Grh) using a 717 bp apoptotic enhancer. In this study, we show that this enhancer is critical for abdominal NB apoptosis and relies on two separable set of DNA binding motifs responsible for its initiation and maintenance. Our results also show that AbdA and Grh interact through their highly conserved DNA binding domains, and the DNA binding specificity of AbdA-HD is important for it to interact with Grh and essential for it to execute NB apoptosis in CNS. We also establish that Grh is required for Hox-dependent NB apoptosis in Labial and Sex Combs Reduced (Scr) expressing regions of the CNS, and it can physically interact with all the Hox proteins in vitro. Our biochemical and functional data collectively support the idea that Grh can function as a Hox cofactor and help them carry out their in vivo roles during development.

Identification of Hox genes and their expression profiles during embryonic development of the emerging model organism, Macrobrachium olfersii

Hox genes encode transcription factors that specify the body segment identity during development, including crustaceans, such as amphipods and decapods, that possess a remarkable diversity of segments and specialized appendages. In amphipods, alterations of specialized appendages have been obtained using knockout experiment of Hox genes, which suggests that these genes are involved in the evolution of morphology within crustaceans. However, studies of Hox genes in crustaceans have been limited to a few species. Here, we identified the homeodomain of nine Hox genes: labial (lab), proboscipedia (pb), Deformed (Dfd), Sex combs reduced (Scr), fushi tarazu (ftz), Antennapedia (Antp), Ultrabithorax (Ubx), abdominal-A (abdA), and Abdominal-B (AbdB), and evaluated their expression by RT-qPCR and RT-PCR in the ovary, during embryonic development, and at the first larval stage (Zoea I) of the decapod Macrobrachium olfersii. The transcript levels of lab, Dfd, and ftz decreased and transcripts of pb, Scr, Antp, Ubx, abdA, and AbdB increased during embryonic development. Hox genes were expressed in mature ovaries and Zoea I larval stages, except Scr and ftz, respectively. In addition, isoforms of Dfd, Scr, Ubx, and abdA, which have been scarcely reported in crustaceans, were described. New partial sequences of 87 Hox genes from other crustaceans were identified from the GenBank database. Our results are interesting for future studies to determine the specific function of Hox genes and their isoforms in the freshwater prawn M. olfersii and to contribute to the understanding of the diversity and evolution of body plans and appendages in Crustaceans.

Embryonic Programs in Cancer and Metastasis—Insights From the Mammary Gland

Cancer is characterized as a reversion of a differentiated cell to a primitive cell state that recapitulates, in many aspects, features of embryonic cells. This review explores the current knowledge of developmental mechanisms that are essential for embryonic mouse mammary gland development, with a particular focus on genes and signaling pathway components that are essential for the induction, morphogenesis, and lineage specification of the mammary gland. The roles of these same genes and signaling pathways in mammary gland or breast tumorigenesis and metastasis are then summarized. Strikingly, key embryonic developmental pathways are often reactivated or dysregulated during tumorigenesis and metastasis in processes such as aberrant proliferation, epithelial-to-mesenchymal transition (EMT), and stem cell potency which affects cellular lineage hierarchy. These observations are in line with findings from recent studies using lineage tracing as well as bulk- and single-cell transcriptomics that have uncovered features of embryonic cells in cancer and metastasis through the identification of cell types, cell states and characterisation of their dynamic changes. Given the many overlapping features and similarities of the molecular signatures of normal development and cancer, embryonic molecular signatures could be useful prognostic markers for cancer. In this way, the study of embryonic development will continue to complement the understanding of the mechanisms of cancer and aid in the discovery of novel therapeutic targets and strategies.

Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

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.

Recent Advances in Syndactyly: Basis, Current Status and Future Perspectives

A comprehensive summary of recent knowledge in syndactyly (SD) is important for understanding the genetic etiology of SD and disease management. Thus, this review article provides background information on SD, as well as insights into phenotypic and genetic heterogeneity, newly identified gene mutations in various SD types, the role of HOXD13 in limb deformities, and recently introduced modern surgical techniques for SD. This article also proposes a procedure for genetic analysis to obtain a clearer genotype-phenotype correlation for SD in the future. We briefly describe the classification of non-syndromic SD based on variable phenotypes to explain different phenotypic features and mutations in the various genes responsible for the pathogenesis of different types of SD. We describe how different types of mutation in HOXD13 cause various types of SD, and how a mutation in HOXD13 could affect its interaction with other genes, which may be one of the reasons behind the differential phenotypes and incomplete penetrance. Furthermore, we also discuss some recently introduced modern surgical techniques, such as free skin grafting, improved flap techniques, and dermal fat grafting in combination with the Z-method incision, which have been successfully practiced clinically with no post-operative complications.

Aberrant overexpression of HOTAIR inhibits abdominal adipogenesis through remodelling of genome-wide DNA methylation and transcription

Objective

Abdominal adiposity is strongly associated with diabetic and cardiovascular comorbidities. The long noncoding RNA HOTAIR (HOX Transcript Antisense Intergenic RNA) is an important epigenetic regulator with fat depot-specific expression. Its functional roles and epigenetic regulation in abdominal adipogenesis remain uncertain.

Methods

We collected different fat depots from healthy, severely obese, and uraemic subjects to measure fat-depot specific gene expression and quantify regional adiposity via dual-energy X-ray absorptiometry (DXA). HOTAIR was overexpressed to evaluate its functional roles. Reduced representation bisulfite sequencing (RRBS), RNA-sequencing, real-time qPCR and RNA/chromatin immunoprecipitation were performed to analyse HOTAIR-mediated epigenetic regulation.

Results

A negative correlation between adipose tissue HOTAIR expression (arm or abdominal subcutaneous fat depots) and regional adiposity under the status of severe obesity or uraemia was observed. HOTAIR overexpression using human immortalized abdominal preadipocytes further revealed its anti-adipogenic effects. Integrative analysis of genome-wide DNA methylation by reduced representation bisulfite sequencing (RRBS) and gene expression was performed. Overall, the differentially methylated genes were functionally enriched for nervous system development, suggesting that HOTAIR may be epigenetically associated with cell lineage commitment. We specifically found that HOTAIR-mediated genes showed strong changes in both DNA methylation and gene expression during abdominal adipogenesis. We observed that two HOTAIR-repressed genes, SLITRK4 and PITPNC1, present an obesity-driven fat-depot specific expression pattern that is positively correlated with the central body fat distribution.

Conclusions

Our study indicated that HOTAIR is a key regulator of abdominal adipogenesis via intricate DNA methylation and is likely to be associated with the transcriptional regulation of genes involved in nervous system development and lipid metabolism, such as SLITRK4 and PITPNC1.

•

HOTAIR was lowly expressed in abdominal and arm fats compared to the gluteal fat.

•

Fat-depot-specific HOTAIR expression could be altered in the obese or uraemic status.

•

HOTAIR overexpression suppressed abdominal adipogenesis and modulated methylome.

•

HOTAIR-suppressed genes were associated with neural development and lipid metabolism.

The Evolution of Biomineralization through the Co-Option of Organic Scaffold Forming Networks

Biomineralization is the process in which organisms use minerals to generate hard structures like teeth, skeletons and shells. Biomineralization is proposed to have evolved independently in different phyla through the co-option of pre-existing developmental programs. Comparing the gene regulatory networks (GRNs) that drive biomineralization in different species could illuminate the molecular evolution of biomineralization. Skeletogenesis in the sea urchin embryo was extensively studied and the underlying GRN shows high conservation within echinoderms, larval and adult skeletogenesis. The organic scaffold in which the calcite skeletal elements form in echinoderms is a tubular compartment generated by the syncytial skeletogenic cells. This is strictly different than the organic cartilaginous scaffold that vertebrates mineralize with hydroxyapatite to make their bones. Here I compare the GRNs that drive biomineralization and tubulogenesis in echinoderms and in vertebrates. The GRN that drives skeletogenesis in the sea urchin embryo shows little similarity to the GRN that drives bone formation and high resemblance to the GRN that drives vertebrates’ vascular tubulogenesis. On the other hand, vertebrates’ bone-GRNs show high similarity to the GRNs that operate in the cells that generate the cartilage-like tissues of basal chordate and invertebrates that do not produce mineralized tissue. These comparisons suggest that biomineralization in deuterostomes evolved through the phylum specific co-option of GRNs that control distinct organic scaffolds to mineralization.

Roles of Drosophila Hox Genes in the Assembly of Neuromuscular Networks and Behavior

Hox genes have been known for specifying the anterior-posterior axis (AP) in bilaterian body plans. Studies in vertebrates have shown their importance in developing region-specific neural circuitry and diversifying motor neuron pools. In Drosophila, they are instrumental for segment-specific neurogenesis and myogenesis early in development. Their robust expression in differentiated neurons implied their role in assembling region-specific neuromuscular networks. In the last decade, studies in Drosophila have unequivocally established that Hox genes go beyond their conventional functions of generating cellular diversity along the AP axis of the developing central nervous system. These roles range from establishing and maintaining the neuromuscular networks to controlling their function by regulating the motor neuron morphology and neurophysiology, thereby directly impacting the behavior. Here we summarize the limited knowledge on the role of Drosophila Hox genes in the assembly of region-specific neuromuscular networks and their effect on associated behavior.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

A chromosome-level genome of Portunus trituberculatus provides insights into its evolution, salinity adaptation and sex determination

Portunus trituberculatus (Crustacea: Decapoda: Brachyura), commonly known as the swimming crab, is of major ecological importance, as well as being important to the fisheries industry. P. trituberculatus is also an important farmed species in China due to its rapid growth rate and high economic value. Here, we report the genome sequence of the swimming crab, which was assembled at the chromosome scale, covering ~1.2 Gb, with 79.99% of the scaffold sequences assembled into 53 chromosomes. The contig and scaffold N50 values were 108.7 kb and 15.6 Mb, respectively, with 19,981 protein-coding genes. Based on comparative genomic analyses of crabs and shrimps, the C2H2 zinc finger protein family was found to be the only gene family expanded in crab genomes, suggesting it was closely related to the evolution of crabs. The combination of transcriptome and bulked segregant analysis provided insights into the genetic basis of salinity adaptation and rapid growth in P. trituberculatus. In addition, the specific region of the Y chromosome was located for the first time in the genome of P. trituberculatus, and three genes were preliminarily identified as candidate genes for sex determination in this region. Decoding the swimming crab genome not only provides a valuable genomic resource for further biological and evolutionary studies, but is also useful for molecular breeding of swimming crabs.

Dissecting the evolutionary role of the Hox gene proboscipedia in Drosophila mouthpart diversification by full locus replacement

Hox genes determine positional codes along the head-to-tail axis. Here, we replaced the entire Drosophila melanogaster proboscipedia (pb) Hox locus, which controls the development of the proboscis and maxillary palps, with that from Drosophila mimica, a related species with highly modified mouthparts. The D. mimica replacement rescues most aspects of adult proboscis morphology; however, the shape and orientation of maxillary palps were modified, resembling D. mimica and closely related species. Expressing the D. mimica Pb protein in the D. melanogaster pattern fully rescued D. melanogaster morphology. However, the expression pattern directed by D. mimica pb cis-regulatory sequences differed from that of D. melanogaster pb in cells that produce altered maxillary structures, indicating that pb regulatory sequences can evolve in related species to alter morphology.

Gene regulation: Context is everything

How transcription factors regulate a distinct set of target genes in different cell types is a fundamental question. A new study demonstrates how Ultrabithorax, a Hox transcription factor, acts as both a repressor and an activator in a cell type-specific manner to alter chromatin accessibility and gene regulation.

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.

An adverse outcome pathway on the disruption of retinoic acid metabolism leading to developmental craniofacial defects

Adverse outcome pathway (AOP) is a conceptual framework that links a molecular initiating event (MIE) via intermediate key events (KEs) with adverse effects (adverse outcomes, AO) relevant for risk assessment, through defined KE relationships (KERs). The aim of the present work is to describe a linear AOP, supported by experimental data, for skeletal craniofacial defects as the AO. This AO was selected in view of its relative high incidence in humans and the suspected relation to chemical exposure. We focused on inhibition of CYP26, a retinoic acid (RA) metabolizing enzyme, as MIE, based on robust previously published data. Conazoles were selected as representative stressors. Intermediate KEs are RA disbalance, aberrant HOX gene expression, disrupted specification, migration, and differentiation of neural crest cells, and branchial arch dysmorphology. We described the biological basis of the postulated events and conducted weight of evidence (WoE) assessments. The biological plausibility and the overall empirical evidence were assessed as high and moderate, respectively, the latter taking into consideration the moderate evidence for concordance of dose-response and temporal relationships. Finally, the essentiality assessment of the KEs, considered as high, supported the robustness of the presented AOP. This AOP, which appears of relevance to humans, thus contributes to mechanistic underpinning of selected test methods, thereby supporting their application in integrated new approach test methodologies and strategies and application in a regulatory context.

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.

Panarthropod tiptop/teashirt and spalt orthologs and their potential role as "trunk"-selector genes

BACKGROUND

In the vinegar fly Drosophila melanogaster, the homeodomain containing transcription factor Teashirt (Tsh) appears to specify trunk identity in concert with the function of the Hox genes. While in Drosophila there is a second gene closely related to tsh, called tiptop (tio), in other arthropods species only one copy exists (called tio/tsh). The expression of tsh and tio/tsh, respectively, is surprisingly similar among arthropods suggesting that its function as trunk selector gene may be conserved. Other research, for example on the beetle Tribolium castaneum, questions even conservation of Tsh function among insects. The zinc-finger transcription factor Spalt (Sal) is involved in the regulation of Drosophila tsh, but this regulatory interaction does not appear to be conserved in Tribolium either. Whether the function and interaction of tsh and sal as potential trunk-specifiers, however, is conserved is still unclear because comparative studies on sal expression (except for Tribolium) are lacking, and functional data are (if at all existing) restricted to Insecta.

RESULTS

Here, we provide additional data on arthropod tsh expression, show the first data on onychophoran tio/tsh expression, and provide a comprehensive investigation on sal expression patterns in arthropods and an onychophoran.

CONCLUSIONS

Our data support the idea that tio/tsh genes are involved in the development of "trunk" segments by regulating limb development. Our data suggest further that the function of Sal is indeed unlikely to be conserved in trunk vs head development like in Drosophila, but early expression of sal is in line with a potential homeotic function, at least in Arthropoda.

Genome of the fatal tapeworm Sparganum proliferum uncovers mechanisms for cryptic life cycle and aberrant larval proliferation

The cryptic parasite Sparganum proliferum proliferates in humans and invades tissues and organs. Only scattered cases have been reported, but S. proliferum infection is always fatal. However, S. proliferum’s phylogeny and life cycle remain enigmatic. To investigate the phylogenetic relationships between S. proliferum and other cestode species, and to examine the mechanisms underlying pathogenicity, we sequenced the entire genomes of S. proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei. Additionally, we performed larvae transcriptome analyses of S. proliferum plerocercoid to identify genes involved in asexual reproduction in the host. The genome sequences confirmed that the S. proliferum has experienced a clearly distinct evolutionary history from S. erinaceieuropaei. Moreover, we found that nonordinal extracellular matrix coordination allows asexual reproduction in the host, and loss of sexual maturity in S. proliferum are responsible for its fatal pathogenicity to humans. Our high-quality reference genome sequences should be valuable for future studies of pseudophyllidean tapeworm biology and parasitism.

Kikuchi et al. sequence the genome of the fatal tapeworm Sparganum proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei, and describe its genomic features suggesting the natural history and molecular mechanisms underlying pathogenicity. Their findings indicate that nonordinal extracellular matrix coordination is important for its asexual reproduction, and suggest that loss of sexual maturity contributes to the fatal pathogenicity of S. proliferum to humans.

Evolutionary divergence of a Hoxa2b hindbrain enhancer in syngnathids mimics results of functional assays

Hoxa2 genes provide critical patterning signals during development, and their regulation and function have been extensively studied. We report a previously uncharacterized significant sequence divergence of a highly conserved hindbrain hoxa2b enhancer element in the family syngnathidae (pipefishes, seahorses, pipehorses, seadragons). We compared the hox cis-regulatory element variation in the Gulf pipefish and two species of seahorse against eight other species of fish, as well as human and mouse. We annotated the hoxa2b enhancer element binding sites across three species of seahorse, four species of pipefish, and one species of ghost pipefish. Finally, we performed in situ hybridization analysis of hoxa2b expression in Gulf pipefish embryos. We found that all syngnathid fish examined share a modified rhombomere 4 hoxa2b enhancer element, despite the fact that this element has been found to be highly conserved across all vertebrates examined previously. Binding element sequence motifs and spacing between binding elements have been modified for the hoxa2b enhancer in several species of pipefish and seahorse, and that the loss of the Prep/Meis binding site and further space shortening happened after ghost pipefish split from the rest of the syngnathid clade. We showed that expression of this gene in rhombomere 4 is lower relative to the surrounding rhombomeres in developing Gulf pipefish embryos, reflecting previously published functional tests for this enhancer. Our findings highlight the benefits of studying highly derived, diverse taxa for understanding of gene regulatory evolution and support the hypothesis that natural mutations can occur in deeply conserved pathways in ways potentially related to phenotypic diversity.

The Chinese mitten crab genome provides insights into adaptive plasticity and developmental regulation

The infraorder Brachyura (true or short-tailed crabs) represents a successful group of marine invertebrates yet with limited genomic resources. Here we report a chromosome-anchored reference genome and transcriptomes of the Chinese mitten crab Eriocheir sinensis, a catadromous crab and invasive species with wide environmental tolerance, strong osmoregulatory capacity and high fertility. We show the expansion of specific gene families in the crab, including F-ATPase, which enhances our knowledge on the adaptive plasticity of this successful invasive species. Our analysis of spatio-temporal transcriptomes and the genome of E. sinensis and other decapods shows that brachyurization development is associated with down-regulation of Hox genes at the megalopa stage when tail shortening occurs. A better understanding of the molecular mechanism regulating sexual development is achieved by integrated analysis of multiple omics. These genomic resources significantly expand the gene repertoire of Brachyura, and provide insights into the biology of this group, and Crustacea in general.

Blueprint of human thymopoiesis reveals molecular mechanisms of stage-specific TCR enhancer activation

Cell differentiation is accompanied by epigenetic changes leading to precise lineage definition and cell identity. Here we present a comprehensive resource of epigenomic data of human T cell precursors along with an integrative analysis of other hematopoietic populations. Although T cell commitment is accompanied by large scale epigenetic changes, we observed that the majority of distal regulatory elements are constitutively unmethylated throughout T cell differentiation, irrespective of their activation status. Among these, the TCRA gene enhancer (Eα) is in an open and unmethylated chromatin structure well before activation. Integrative analyses revealed that the HOXA5-9 transcription factors repress the Eα enhancer at early stages of T cell differentiation, while their decommission is required for TCRA locus activation and enforced αβ T lineage differentiation. Remarkably, the HOXA-mediated repression of Eα is paralleled by the ectopic expression of homeodomain-related oncogenes in T cell acute lymphoblastic leukemia. These results highlight an analogous enhancer repression mechanism at play in normal and cancer conditions, but imposing distinct developmental constraints.

A developmental perspective of homology and evolutionary novelty

The development and evolution of multicellular body plans is complex. Many distinct organs and body parts must be reproduced at each generation, and those that are traceable over long time scales are considered homologous. Among the most pressing and least understood phenomena in evolutionary biology is the mode by which new homologs, or "novelties" are introduced to the body plan and whether the developmental changes associated with such evolution deserve special treatment. In this chapter, we address the concepts of homology and evolutionary novelty through the lens of development. We present a series of case studies, within insects and vertebrates, from which we propose a developmental model of multicellular organ identity. With this model in hand, we make predictions regarding the developmental evolution of body plans and highlight the need for more integrative analysis of developing systems.

Accumulation of cis- and trans-regulatory variations is associated with phenotypic divergence of a complex trait between yeast species

Gene regulatory variations accumulate during evolution and alter gene expression. While the importance of expression variation in phenotypic evolution is well established, the molecular basis remains largely unknown. Here, we examine two closely related yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus, which show phenotypical differences in morphology and cell cycle progression when grown in the same environment. By profiling the cell cycle transcriptome and binding of key transcription factors (TFs) in the two species and their hybrid, we show that changes in expression levels and dynamics of oscillating genes are dominated by upstream trans-variations. We find that multiple cell cycle regulators show both cis- and trans-regulatory variations, which alters their expression in favor of the different cell cycle phenotypes. Moreover, we show that variations in the cell cycle TFs, Fkh1, and Fkh2 affect both the expression of target genes, and the binding specificity of an interacting TF, Ace2. Our study reveals how multiple variations accumulate and propagate through the gene regulatory network, alter TFs binding, contributing to phenotypic changes in cell cycle progression.