Hox clusters as models for vertebrate genome evolution

Evolution and expression of TLR5a and TLR5b in lamprey (Lampetra japonica)

The lamprey serves as a key model organism for studying the origin and evolution of species, embryonic development, and the immune system. The immune system primarily relies on pattern recognition receptors, including Toll-like receptors (TLRs), with Toll-like receptor 5 (TLR5) having a particularly complex evolutionary history. Currently, although TLR5 is being identified in an expanding array of taxonomic groups, a comprehensive study on its evolutionary aspects is yet to be conducted. In this study, we identified Lj-TLR5a and Lj-TLR5b in Lampetra japonica and examined their distribution across various tissues in this species. Furthermore, we conducted preliminary investigations into their immune functions and discovered that, as primitive genes, they are highly sensitive to various pathogens. Upon recognizing flagellar proteins, both Lj-TLR5a and Lj-TLR5b work together; however, these TLRs may function independently in response to other stimuli. Subsequently, we performed comprehensive structural and evolutionary analyses of the TLR5 family, incorporating TLR5 data from various species at different evolutionary stages. Our findings revealed that TLR5a and TLR5b in lamprey are in a relatively primitive evolutionary state. Meanwhile, TLR5L differentiated during the early stages of evolution and exhibits a pseudogenic trend throughout this evolutionary process; notably, this TLR is currently preserved only in certain amphibian and reptile species. In cartilaginous fish, only one type of TLR5 is usually retained, whereas bony fish typically possess both TLR5 and TLR5S. TLR5S has a relatively simple structure, likely arising from repetitive whole-genome events in bony fishes. Among bony fishes, flesh-finned fish were found to retain only one TLR5, which eventually evolved into the TLR5 found in quadrupeds. In summary, this study provides significant insights into the origin and evolution of the TLR5 family by analyzing the evolutionary status and immune functions of Lj-TLR5a and Lj-TLR5b in Japanese lampreys.

HOXC12 promotes the invasion and migration of gastric cancer cells by upregulating SALL4 and activating Wnt/β-catenin signaling pathway

BACKGROUND

Gastric cancer is one of the most common malignant tumors in the world, with a poor prognosis. HOXC is a family of transcription factors that are up-regulated in gastric cancer tissues. However, the relationship between Homeobox C12 (HOXC12) and gastric cancer is still unclear.

METHODS

TCGA-STAD and HPA data were analyzed to explore HOXC12 level. Kaplan-Meier Plotter was used to analyze the relationship between HOXC12 level and the prognosis of gastric cancer patients. The HOXC12 was knocked down or overexpressed by shRNA or overexpression vector to explore its functions. Cell migration/invasion assays and wound healing assay were used to assess the invasion/migration ability of gastric cancer cells. Western blot and qPCR were used to detect gene expression and the activation of Wnt/β-catenin signaling pathway. Dual-luciferase reporter assay was used to detect the active region bound by HOXC12 in the promoter of Spalt-like transcription factor 4 (SALL4).

RESULTS

HOXC12 was highly expressed in gastric cancer and was positively correlated with the poor prognosis of gastric cancer patients. HOXC12 promotes the invasion and migration of gastric cancer cells via Wnt/β-catenin signaling pathway. HOXC12 upregulated the transcription of SALL4 by binding to its promoter. HOXC12 was negatively correlated with both the levels of CD8+ T cells and T cell cytotoxicity-related genes.

CONCLUSION

HOXC12 promotes the invasion/migration of gastric cancer cells via SALL4/Wnt/β-catenin axis, and is negatively correlated with the infiltration of CD8+ T cells, suggesting HOXC12 as a diagnostic marker and a potential therapeutic target for gastric cancer.

Evolutionary origin and distribution of leucine-rich repeat-containing G protein-coupled receptors in insects

Leucine-rich repeat-containing G protein-coupled receptors (LGRs) are crucial for animal growth and development. They were categorized into four types (A, B, C1, and C2) based on their sequence and domain structures. Despite the widespread distribution of LGRs across bilaterians, a thorough investigation of their distribution and evolutionary history remains elusive. Recent studies insect LGRs, especially the emergence of type C2 LGRs in various hemimetabolous insects, had prompted our study to address these problems. Initially, we traced the origins of LGRs by exploiting data from 99 species spanning 11 metazoan phyla, and discovered that type A and B LGRs originated from sponges, while type C LGRs originated from cnidarians. Subsequently, through comprehensive genomic and transcriptomic analyses across 565 species across 25 orders of insects, we found that both type A and C1 LGRs divided into two gene clusters. These clusters can be traced back to basal Insecta and an early ancestor of the Arthropoda, respectively. Furthermore, the absence of type B LGRs in wingless insects suggests a role in wing development, while the absence of type C2 LGRs in holometabolous insects hints at novel functions unrelated to insect metamorphosis. According to the origin of LGRs and the investigation of LGRs in insects, we speculated that type A and B LGRs appeared first among four types of LGRs, type A evolved into type C LGRs later, and type A and C1 LGRs independently duplicated during the evolutionary process. This study provides a more comprehensive view of the evolution of LGR genes than previously available, and sheds light on the evolutionary history and significance of LGRs in insect biology.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

Continuity and change in neural plasticity through embryonic morphogenesis, fetal activity-dependent synaptogenesis, and infant memory consolidation

There is an apparent continuity in human neural development that can be traced to venerable themes of vertebrate morphogenesis that have shaped the evolution of the reptilian telencephalon (including both primitive three-layered cortex and basal ganglia) and then the subsequent evolution of the mammalian six-layered neocortex. In this theoretical analysis, we propose that an evolutionary-developmental analysis of these general morphogenetic themes can help to explain the embryonic development of the dual divisions of the limbic system that control the dorsal and ventral networks of the human neocortex. These include the archicortical (dorsal limbic) Papez circuits regulated by the hippocampus that organize spatial, contextual memory, as well as the paleocortical (ventral limbic) circuits that organize object memory. We review evidence that these dorsal and ventral limbic divisions are controlled by the differential actions of brainstem lemnothalamic and midbrain collothalamic arousal control systems, respectively, thereby traversing the vertebrate subcortical neuraxis. These dual control systems are first seen shaping the phyletic morphogenesis of the archicortical and paleocortical foundations of the forebrain in embryogenesis. They then provide dual modes of activity-dependent synaptic organization in the active (lemnothalamic) and quiet (collothalamic) stages of fetal sleep. Finally, these regulatory systems mature to form the major systems of memory consolidation of postnatal development, including the rapid eye movement (lemnothalamic) consolidation of implicit memory and social attachment in the first year, and then-in a subsequent stage-the non-REM (collothalamic) consolidation of explicit memory that is integral to the autonomy and individuation of the second year of life.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Comparative Studies on Duplicated foxl2 Paralogs in Spotted Knifejaw Oplegnathus punctatus Show Functional Diversification

As a member of the forkhead box L gene family, foxl2 plays a significant role in gonadal development and the regulation of reproduction. During the evolution of deuterostome, whole genome duplication (WGD)-enriched lineage diversifications and regulation mechanisms occurs. However, only limited research exists on foxl2 duplication in teleost or other vertebrate species. In this study, two foxl2 paralogs, foxl2 and foxl2l, were identified in the transcriptome of spotted knifejaw (Oplegnathus punctatus), which had varying expressions in the gonads. The foxl2 was expressed higher in the ovary, while foxl2l was expressed higher in the testis. Phylogenetic reconstruction, synteny analysis, and the molecular evolution test confirmed that foxl2 and foxl2l likely originated from the first two WGD. The expression patterns test using qRT-PCR and ISH as well as motif scan analysis revealed evidence of potentially functional divergence between the foxl2 and foxl2l paralogs in spotted knifejaw. Our results indicate that foxl2 and foxl2l may originate from the first two WGD, be active in transcription, and have undergone functional divergence. These results shed new light on the evolutionary trajectories of foxl2 and foxl2l and highlights the need for further detailed functional analysis of these two duplicated paralogs.

Chicken HOXC8 and HOXC10 genes may play a role in the altered skull morphology associated with the Crest phenotype

One of the most intriguing traits found in domestic chickens is the Crest phenotype. This trait, characterized by a tuft of elongated feathers sprouted from the head, is found in breeds such as Polish chickens and Silkie chickens. Moreover, some crested chicken breeds also exhibit a protuberance in their anterodorsal skull region. Previous studies have strived to identify the causative factors of this trait. This study aimed to elucidate the role of chicken HOXC8 and HOXC10 in the formation of the Crest phenotype. We explored the effect of ectopic expression of HOXC8 or HOXC10 on the chicken craniofacial morphology using the RCAS retrovirus transformation system. Microcomputed tomography scanning was conducted to measure the 3D structure of the cranial bone of transgenic embryos for geometric morphometric analysis. We found that the ectopic expression of HOXC8 or HOXC10 in chicken heads caused mild morphological changes in the skull compared with the GFP-transgenic control group. Geometric morphometric analysis showed that HOXC8 and HOXC10 transgenic groups expressed a mild upward shape change in the frontal region of the skull compared with the control group, which is similar to what is seen in the crested chicken breeds. In conclusion, this study supports findings in previous studies in which HOX genes play a role in the formation of the altered skull morphology related to the Crest phenotype. It also supports that mutations in HOX genes may contribute to intra- and inter-specific variation in morphological traits in vertebrates.

Zebrafish as a Model to Study Retinoic Acid Signaling in Development and Disease

Retinoic acid (RA) is a metabolite of vitamin A (retinol) that plays various roles in development to influence differentiation, patterning, and organogenesis. RA also serves as a crucial homeostatic regulator in adult tissues. The role of RA and its associated pathways are well conserved from zebrafish to humans in both development and disease. This makes the zebrafish a natural model for further interrogation into the functions of RA and RA-associated maladies for the sake of basic research, as well as human health. In this review, we explore both foundational and recent studies using zebrafish as a translational model for investigating RA from the molecular to the organismal scale.

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.

Genome-wide identification of endosialin family of C-type lectins in common carp (Cyprinus carpio) and their response following Aeromonas hydrophila infection

The endosialin family is the group XIV of C-type lectin, regulating several processes involved in innate immunity and inflammation. Endosialin family genes have been extensively studied in human and mammals, however, rarely reported in teleost. In the present study, a set of 8 endosialin family genes was identified across the entire common carp genome. Functional domain and motif prediction and phylogenetic analysis supported their annotation and orthologies. Through examining gene copy number across several vertebrates, endosialin family genes were found have undergone gene duplication. Most of the endosialin family genes were ubiquitously expressed during common carp early developmental stages, and presented tissue-specific expression patterns in various healthy tissues, with relatively high expression in intestine, liver, gill, spleen and kidney, indicating their likely essential roles in maintaining homeostasis and host immune response. After Aeromonas hydrophila infection, gene thbd-1, thbd-2 and cd93-2 were significantly up-regulated at one or more timepoints in spleen and kidney, while gene cd248a-1, cd248a-2, cd248b-1, cd248b-2, and cd93-1 were significantly down-regulated. Taken together, all these results suggested that endosialin family genes were involved in host immune response to A. hydrophila infection in common carp, and provided fundamental genomic resources for better understanding the critical roles of endosialin family on the primary innate immune processes in teleost.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

Parallel Genomic Changes Drive Repeated Evolution of Placentas in Live-Bearing Fish

The evolutionary origin of complex organs challenges empirical study because most organs evolved hundreds of millions of years ago. The placenta of live-bearing fish in the family Poeciliidae represents a unique opportunity to study the evolutionary origin of complex organs, because in this family a placenta evolved at least nine times independently. It is currently unknown whether this repeated evolution is accompanied by similar, repeated, genomic changes in placental species. Here, we compare whole genomes of 26 poeciliid species representing six out of nine independent origins of placentation. Evolutionary rate analysis revealed that the evolution of the placenta coincides with convergent shifts in the evolutionary rate of 78 protein-coding genes, mainly observed in transporter- and vesicle-located genes. Furthermore, differences in sequence conservation showed that placental evolution coincided with similar changes in 76 noncoding regulatory elements, occurring primarily around genes that regulate development. The unexpected high occurrence of GATA simple repeats in the regulatory elements suggests an important function for GATA repeats in developmental gene regulation. The distinction in molecular evolution observed, with protein-coding parallel changes more often found in metabolic and structural pathways, compared with regulatory change more frequently found in developmental pathways, offers a compelling model for complex trait evolution in general: changing the regulation of otherwise highly conserved developmental genes may allow for the evolution of complex traits.

Highly conserved and cis-acting lncRNAs produced from paralogous regions in the center of HOXA and HOXB clusters in the endoderm lineage

Long noncoding RNAs (lncRNAs) have been shown to play important roles in gene regulatory networks acting in early development. There has been rapid turnover of lncRNA loci during vertebrate evolution, with few human lncRNAs conserved beyond mammals. The sequences of these rare deeply conserved lncRNAs are typically not similar to each other. Here, we characterize HOXA-AS3 and HOXB-AS3, lncRNAs produced from the central regions of the HOXA and HOXB clusters. Sequence-similar orthologs of both lncRNAs are found in multiple vertebrate species and there is evident sequence similarity between their promoters, suggesting that the production of these lncRNAs predates the duplication of the HOX clusters at the root of the vertebrate lineage. This conservation extends to similar expression patterns of the two lncRNAs, in particular in cells transiently arising during early development or in the adult colon. Functionally, the RNA products of HOXA-AS3 and HOXB-AS3 regulate the expression of their overlapping HOX5-7 genes both in HT-29 cells and during differentiation of human embryonic stem cells. Beyond production of paralogous protein-coding and microRNA genes, the regulatory program in the HOX clusters therefore also relies on paralogous lncRNAs acting in restricted spatial and temporal windows of embryonic development and cell differentiation.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

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.

Grand Challenges in Comparative Tooth Biology

Teeth are a model system for integrating developmental genomics, functional morphology, and evolution. We are at the cusp of being able to address many open issues in comparative tooth biology and we outline several of these newly tractable and exciting research directions. Like never before, technological advances and methodological approaches are allowing us to investigate the developmental machinery of vertebrates and discover both conserved and excitingly novel mechanisms of diversification. Additionally, studies of the great diversity of soft tissues, replacement teeth, and non-trophic functions of teeth are providing new insights into dental diversity. Finally, we highlight several emerging model groups of organisms that are at the forefront of increasing our appreciation of the mechanisms underlying tooth diversification.

The Invariant Nature of a Morphological Character and Character State: Insights from Gene Regulatory Networks

What constitutes a discrete morphological character versus character state has been long discussed in the systematics literature but the consensus on this issue is still missing. Different methods of classifying organismal features into characters and character states (CCSs) can dramatically affect the results of phylogenetic analyses. Here, I show that, in the framework of Markov models, the modular structure of the gene regulatory network (GRN) underlying trait development, and the hierarchical nature of GRN evolution, essentially remove the distinction between morphological CCS, thus endowing the CCS with an invariant property with respect to each other. This property allows the states of one character to be represented as several individual characters and vice versa. In practice, this means that a phenotype can be encoded using a set of characters or just one complex character with numerous states. The representation of a phenotype using one complex character can be implemented in Markov models of trait evolution by properly structuring transition rate matrix.

Hox genes reveal variations in the genomic DNA of allotetraploid hybrids derived from Carassius auratus red var. (female) × Cyprinus carpio L. (male)

Background

Hox transcription factors are master regulators of animal development. Although highly conserved, they can contribute to the formation of novel biological characteristics when modified, such as during the generation of hybrid species, thus potentially serving as species-specific molecular markers. Here, we systematically studied the evolution of genomic sequences of Hox loci in an artificial allotetraploid lineage (4nAT, 4n = 200) derived from a red crucian carp (♀, RCC, 2n = 100) × common carp (♂, CC, 2n = 100) cross and its parents (RCC and CC).

Results

PCR amplification yielded 23 distinct Hox gene fragments from 160 clones in 4nAT, 22 fragments from 90 clones in RCC, and 19 fragments from 90 clones in CC. Sequence alignment of the HoxA3a and HoxC10a genes indicated both the inheritance and loss of paternal genomic DNA in 4nAT. The HoxA5a gene from 4nAT consisted of two subtypes from RCC and two subtypes from CC, indicating that homologous recombination occurred in the 4nAT hybrid genome. Moreover, 4nAT carried genomic pseudogenization in the HoxA10b and HoxC13a loci. Interestingly, a new type of HoxC9a gene was found in 4nAT as a hybrid sequence of CC and RCC by recombination in the intronic region.

Conclusion

The results revealed the influence of Hox genes during polyploidization in hybrid fish. The data provided insight into the evolution of vertebrate genomes and might be benefit for artificial breeding programs.

Homeobox Gene Duplication and Divergence in Arachnids

Homeobox genes are key toolkit genes that regulate the development of metazoans and changes in their regulation and copy number have contributed to the evolution of phenotypic diversity. We recently identified a whole genome duplication (WGD) event that occurred in an ancestor of spiders and scorpions (Arachnopulmonata), and that many homeobox genes, including two Hox clusters, appear to have been retained in arachnopulmonates. To better understand the consequences of this ancient WGD and the evolution of arachnid homeobox genes, we have characterized and compared the homeobox repertoires in a range of arachnids. We found that many families and clusters of these genes are duplicated in all studied arachnopulmonates (Parasteatoda tepidariorum, Pholcus phalangioides, Centruroides sculpturatus, and Mesobuthus martensii) compared with nonarachnopulmonate arachnids (Phalangium opilio, Neobisium carcinoides, Hesperochernes sp., and Ixodes scapularis). To assess divergence in the roles of homeobox ohnologs, we analyzed the expression of P. tepidariorum homeobox genes during embryogenesis and found pervasive changes in the level and timing of their expression. Furthermore, we compared the spatial expression of a subset of P. tepidariorum ohnologs with their single copy orthologs in P. opilio embryos. We found evidence for likely subfunctionlization and neofunctionalization of these genes in the spider. Overall our results show a high level of retention of homeobox genes in spiders and scorpions post-WGD, which is likely to have made a major contribution to their developmental evolution and diversification through pervasive subfunctionlization and neofunctionalization, and paralleling the outcomes of WGD in vertebrates.

Long non-coding RNA HOTAIR and HOTTIP as potential biomarkers for hepatitis C virus genotype 4-induced hepatocellular carcinoma

Background

Long non-coding RNAs (lncRNAs) homeobox (Hox) transcript antisense intergenic RNA (HOTAIR) and HOXA transcript at the distal tip (HOTTIP) have been suggested to be implicated in liver cancer tumorigenesis and progression; however, little is known about the role of the plasma HOTAIR and HOTTIP in liver cancer diagnosis and prognosis. The current study aimed at measuring the plasma levels of long non-coding RNAs (HOTAIR and HOTTIP) expression in chronic liver disease (CLD) due to HCV genotype 4 infection with/without cirrhosis and HCC patients in an attempt to evaluate the potential benefits of these new circulating as non-invasive diagnostic biomarkers and a novel therapeutic strategy for liver cirrhosis and carcinogenesis of Egyptian patients. Hundred subjects were included in this study, divided into two groups; group I (50 patients) were classified into subgroup Ia (CLD without cirrhosis, n = 25) and subgroup Ib (CLD with cirrhosis, n = 25), group II (CLD patients with HCC, n = 25), and control (healthy volunteer, n = 25). The expression of lncRNAs (HOTAIR and HOTTIP) genes was analyzed by real-time PCR.

Results

LncRNAs (HOTAIR and HOTTIP) showed upregulation in all diseased groups, which was in consistent with the progression of the disease toward the HCC stage. In addition, HOTAIR and HOTTIP showed a diagnostic ability to discriminate between cases of cirrhosis and HCC compared with healthy control (p < 0.001), while HOTAIR and HOTTIP did not show a discrimination significant differences between cirrhotic cases and non-cirrhotic cases. By using receiver operating characteristic curve (ROC) analysis, it was found that LncRNAs (HOTAIR and HOTTIP) could diagnose liver cancer with 64.0% sensitivity and 86.0% specificity and 48.0% sensitivity and 88.0% specificity. Furthermore, both genes can be considered as the predictor and prognostic parameters for cirrhosis (OR = 1.111, p = 0.05) and (OR = 1.07, p = 0.05) respectively, and HCC (OR = 1.047, p = 0.01) and (OR = 1.05, p = 0.003). The increased HOTAIR and HOTTIP expression were associated with advanced tumor stages and higher grades.

Conclusion

These results strongly prompt us that HOTAIR and HOTTIP genes can be used as non-invasive prognostic biomarkers and new therapeutic targets for HCV genotype 4-induced HCC.

Deciphering The Potential Role of Hox Genes in Pancreatic Cancer

The Hox gene family plays an important role in organogenesis and animal development. Currently, 39 Hox genes that are clustered in four chromosome regions have been identified in humans. Emerging evidence suggests that Hox genes are involved in the development of the pancreas. However, the expression of Hox genes in pancreatic tumor tissues has been investigated in only a few studies. In addition, whether specific Hox genes can promote or suppress cancer metastasis is not clear. In this article, we first review the recent progress in studies on the role of Hox genes in pancreatic cancer. By comparing the expression profiles of pancreatic cancer cells isolated from genetically engineered mice established in our laboratory with three different proliferative and metastatic abilities, we identified novel Hox genes that exhibited tumor-promoting activity in pancreatic cancer. Finally, a potential oncogenic mechanism of the Hox genes was hypothesized.

Unravelling paralogous gene expression dynamics during three-spined stickleback embryogenesis

Development requires the implementation of a plethora of molecular mechanisms, involving a large set of genes to ensure proper cell differentiation, morphogenesis of tissues and organs as well as the growth of the organism. Genome duplication and resulting paralogs are considered to provide the raw genetic materials important for new adaptation opportunities and boosting evolutionary innovation. The present study investigated paralogous genes, involved in three-spined stickleback (Gasterosteus aculeatus) development. Therefore, the transcriptomes of five early stages comprising developmental leaps were explored. Obtained expression profiles reflected the embryo's needs at different stages. Early stages, such as the morula stage comprised transcripts mainly involved in energy requirements while later stages were mostly associated with GO terms relevant to organ development and morphogenesis. The generated transcriptome profiles were further explored for differential expression of known and new paralogous genes. Special attention was given to hox genes, with hoxa13a being of particular interest and to pigmentation genes where itgb1, involved in the melanophore development, displayed a complementary expression pattern throughout studied stages. Knowledge obtained by untangling specific paralogous gene functions during development might not only significantly contribute to the understanding of teleost ontogenesis but might also shed light on paralogous gene evolution.

Clinical significance of long non-coding RNA HOTTIP in early-stage non-small-cell lung cancer

Background

HOTTIP, a long non-coding RNA located in the HOXA cluster, plays a role in the patterning of tissues with mesodermal components, including the lung. Overexpression of HOXA genes, including HOTTIP, has been associated with a more aggressive phenotype in several cancers. However, the prognostic impact of HOTTIP has not yet been explored in non-small-cell lung cancer (NSCLC). We have correlated HOTTIP expression with time to relapse (TTR) and overall survival (OS) in early-stage NSCLC patients.

Methods

Ninety-nine early-stage NSCLC patients who underwent surgical resection in our center from June 2007 to November 2013 were included in the study. Mean age was 66; 77.8% were males; 73.7% had stage I disease; and 55.5% had adenocarcinoma. A validation data set comprised stage I-II patients from The Cancer Genome Atlas (TCGA) Research Network.

Results

HOTTIP was expressed in all tumor samples and was overexpressed in squamous cell carcinoma (p = 0.007) and in smokers (p = 0.018). Patients with high levels of HOTTIP had shorter TTR (78.3 vs 58 months; p = 0.048) and shorter OS (81.2 vs 61 months; p = 0.023) than those with low levels. In the multivariate analysis, HOTTIP emerged as an independent prognostic marker for TTR (OR: 2.05, 95%CI: 1–4.2; p = 0.05), and for OS (OR: 2.31, 95%CI: 1.04–5.1; p = 0.04). HOTTIP was validated as a prognostic marker for OS in the TCGA adenocarcinoma cohort (p = 0.025). Moreover, we identified a 1203-mRNA and a 61-miRNA signature that correlated with HOTTIP expression.

Conclusions

The lncRNA HOTTIP can be considered a prognostic biomarker in early-stage NSCLC.

Electronic supplementary material

The online version of this article (10.1186/s12890-019-0816-8) contains supplementary material, which is available to authorized users.

De novo draft assembly of the Botrylloides leachii genome provides further insight into tunicate evolution

Tunicates are marine invertebrates that compose the closest phylogenetic group to the vertebrates. These chordates present a particularly diverse range of regenerative abilities and life-history strategies. Consequently, tunicates provide an extraordinary perspective into the emergence and diversity of these traits. Here we describe the genome sequencing, annotation and analysis of the Stolidobranchian Botrylloides leachii. We have produced a high-quality 159 Mb assembly, 82% of the predicted 194 Mb genome. Analysing genome size, gene number, repetitive elements, orthologs clustering and gene ontology terms show that B. leachii has a genomic architecture similar to that of most solitary tunicates, while other recently sequenced colonial ascidians have undergone genome expansion. In addition, ortholog clustering has identified groups of candidate genes for the study of colonialism and whole-body regeneration. By analysing the structure and composition of conserved gene linkages, we observed examples of cluster breaks and gene dispersions, suggesting that several lineage-specific genome rearrangements occurred during tunicate evolution. We also found lineage-specific gene gain and loss within conserved cell-signalling pathways. Such examples of genetic changes within conserved cell-signalling pathways commonly associated with regeneration and development that may underlie some of the diverse regenerative abilities observed in tunicates. Overall, these results provide a novel resource for the study of tunicates and of colonial ascidians.

Evolution of the α2-adrenoreceptors in vertebrates: ADRA2D is absent in mammals and crocodiles

Evolutionary studies of genes that have been functionally characterized and whose variation has been associated with pathological conditions represent an opportunity to understand the genetic basis of pathologies. α₂-Adrenoreceptors (ADRA2) are a class of G protein-coupled receptors that regulate several physiological processes including blood pressure, platelet aggregation, insulin secretion, lipolysis, and neurotransmitter release. This gene family has been extensively studied from a molecular/physiological perspective, yet much less is known about its evolutionary history. Accordingly, the goal of this study was to investigate the evolutionary history of α₂-adrenoreceptors (ADRA2) in vertebrates. Our results show that in addition to the three well-recognized α₂-adrenoreceptor genes (ADRA2A, ADRA2B and ADRA2C), we recovered a clade that corresponds to the fourth member of the α₂-adrenoreceptor gene family (ADRA2D). We also recovered a clade that possesses two ADRA2 sequences found in two lamprey species. Furthermore, our results show that mammals and crocodiles are characterized by possessing three α₂-adrenoreceptor genes, whereas all other vertebrate groups possess the full repertoire of α₂-adrenoreceptor genes. Among vertebrates ADRA2D seems to be a dispensable gene, as it was lost two independent times during the evolutionary history of the group. Additionally, we found that most examined species possess the most common alleles described for humans; however, there are cases in which non-human mammals possess the alternative variant. Finally, transcript abundance profiles revealed that during the early evolutionary history of gnathostomes, the expression of ADRA2D in different taxonomic groups became specialized to different tissues, but in the ancestor of sarcopterygians this specialization would have been lost.

Reprimo tissue-specific expression pattern is conserved between zebrafish and human

Reprimo (RPRM), a member of the RPRM gene family, is a tumor-suppressor gene involved in the regulation of the p53-mediated cell cycle arrest at G2/M. RPRM has been associated with malignant tumor progression and proposed as a potential biomarker for early cancer detection. However, the expression and role of RPRM, as well as its family, are poorly understood and their physiology is as yet unstudied. In this scenario, a model system like the zebrafish could serve to dissect the role of the RPRM family members in vivo. Phylogenetic analysis reveals that RPRM and RPRML have been differentially retained by most species throughout vertebrate evolution, yet RPRM3 has been retained only in a small group of distantly related species, including zebrafish. Herein, we characterized the spatiotemporal expression of RPRM (present in zebrafish as an infraclass duplication rprma/rprmb), RPRML and RPRM3 in the zebrafish. By whole-mount in situ hybridization (WISH) and fluorescent in situ hybridization (FISH), we demonstrate that rprm (rprma/rprmb) and rprml show a similar spatiotemporal expression profile during zebrafish development. At early developmental stages rprmb is expressed in somites. After one day post-fertilization, rprm (rprma/rprmb) and rprml are expressed in the notochord, brain, blood vessels and digestive tube. On the other hand, rprm3 shows the most unique expression profile, being expressed only in the central nervous system (CNS). We assessed the expression patterns of RPRM gene transcripts in adult zebrafish and human RPRM protein product in tissue samples by RT-qPCR and immunohistochemistry (IHC) staining, respectively. Strikingly, tissue-specific expression patterns of the RPRM transcripts and protein are conserved between zebrafish and humans. We propose the zebrafish as a powerful tool to elucidate the both physiological and pathological roles of the RPRM gene family.

How does chromatin package DNA within nucleus and regulate gene expression?

The human body is made up of 60 trillion cells, each cell containing 2 millions of genomic DNA in its nucleus. How is this genomic deoxyribonucleic acid [DNA] organised into nuclei? Around 1880, W. Flemming discovered a nuclear substance that was clearly visible on staining under primitive light microscopes and named it 'chromatin'; this is now thought to be the basic unit of genomic DNA organization. Since long before DNA was known to carry genetic information, chromatin has fascinated biologists. DNA has a negatively charged phosphate backbone that produces electrostatic repulsion between adjacent DNA regions, making it difficult for DNA to fold upon itself. In this article, we will try to shed light on how does chromatin package DNA within nucleus and regulate gene expression?

Evolution of the β-adrenoreceptors in vertebrates

The study of the evolutionary history of genes related to human disease lies at the interface of evolution and medicine. These studies provide the evolutionary context on which medical researchers should work, and are also useful in providing information to suggest further genetic experiments, especially in model species where genetic manipulations can be made. Here we studied the evolution of the β-adrenoreceptor gene family in vertebrates with the aim of adding an evolutionary framework to the already abundant physiological information. Our results show that in addition to the three already described vertebrate β-adrenoreceptor genes there is an additional group containing cyclostome sequences. We suggest that β-adrenoreceptors diversified as a product of the two whole genome duplications that occurred in the ancestor of vertebrates. Gene expression patterns are in general consistent across species, suggesting that expression dynamics were established early in the evolutionary history of vertebrates, and have been maintained since then. Finally, amino acid polymorphisms that are associated to pathological conditions in humans appear to be common in non-human mammals, suggesting that the phenotypic effects of these mutations depend on epistatic interaction with other positions. The evolutionary analysis of the β-adrenoreceptors delivers new insights about the diversity of these receptors in vertebrates, the evolution of the expression patterns and a comparative perspective regarding the polymorphisms that in humans are linked to pathological conditions.

Comparative Physiology of Energy Metabolism: Fishing for Endocrine Signals in the Early Vertebrate Pool

Energy is the common currency of life. To guarantee a homeostatic supply of energy, multiple neuro-endocrine systems have evolved in vertebrates; systems that regulate food intake, metabolism, and distribution of energy. Even subtle (lasting) dysregulation of the delicate balance of energy intake and expenditure may result in severe pathologies. Feeding-related pathologies have fueled research on mammals, including of course the human species. The mechanisms regulating food intake and body mass are well-characterized in these vertebrates. The majority of animal life is ectothermic, only birds and mammals are endotherms. What can we learn from a (comparative) study on energy homeostasis in teleostean fishes, ectotherms, with a very different energy budget and expenditure? We present several adaptation strategies in fish. In recent years, the components that regulate food intake in fishes have been identified. Although there is homology of the major genetic machinery with mammals (i.e., there is a vertebrate blueprint), in many cases this does not imply analogy. Although both mammals and fish must gain their energy from food, the expenditure of the energy obtained is different. Mammals need to spend vast amounts of energy to maintain body temperature; fishes seem to utilize a broader metabolic range to their advantage. In this review, we briefly discuss ecto- and endothermy and their consequences for energy balance. Next, we argue that the evolution of endothermy and its (dis-)advantages may explain very different strategies in endocrine regulation of energy homeostasis among vertebrates. We follow a comparative and evolutionary line of thought: we discuss similarities and differences between fish and mammals. Moreover, given the extraordinary radiation of teleostean fishes (with an estimated number of 33,400 contemporary species, or over 50% of vertebrate life forms), we also compare strategies in energy homeostasis between teleostean species. We present recent developments in the field of (neuro)endocrine regulation of energy balance in teleosts, with a focus on leptin.

HOXB4 knockdown reverses multidrug resistance of human myelogenous leukemia K562/ADM cells by downregulating P-gp, MRP1 and BCRP expression via PI3K/Akt signaling pathway

Multidrug resistance (MDR) plays a pivotal role in human chronic myelogenous leukemia (CML) chemotherapy failure. MDR is mainly associated with the overexpression of drug efflux transporters of the ATP-binding cassette (ABC) proteins. Phosphoinositide 3-kinase (PI3K)/Akt signaling cascade is involved in the MDR phenotype and is correlated with multidrug resistance 1 (MDR1)/P-glycoprotein (P-gp), multidrug resistance-associated protein 1 (MRP1) and breast cancer resistance protein (BCRP) expression in many human malignancies. Homeobox (HOX) B4, a member of the HOX gene family, has been reported to be correlated with occurrence, development, poor prognosis and drug resistance of human leukemia. In the present study, HOXB4 expression was analyzed in K562 cell line and its MDR subline K562/ADM. Compared with K562 cells, drug-resistant K562/ADM cells demonstrated evidently higher HOXB4 expression. In addition, we firstly investigated the reversal effect of HOXB4 deletion on K562/ADM cells and the underlying mechanism. The Cell Counting kit-8 (CCK-8) and flow cytometry assays showed that knockdown of HOXB4 enhanced chemosensitivity and decreased drug efflux in K562/ADM cells. Moreover, HOXB4 knockout led to downregulation of P-gp, MRP1 and BCRP expression and PI3K/Akt signaling activity, suggesting that repression of HOXB4 might be a key point to reverse MDR of K562/ADM cells.

Evolution history of duplicated smad3 genes in teleost: insights from Japanese flounder, Paralichthys olivaceus

Following the two rounds of whole-genome duplication (WGD) during deuterosome evolution, a third genome duplication occurred in the ray-fined fish lineage and is considered to be responsible for the teleost-specific lineage diversification and regulation mechanisms. As a receptor-regulated SMAD (R-SMAD), the function of SMAD3 was widely studied in mammals. However, limited information of its role or putative paralogs is available in ray-finned fishes. In this study, two SMAD3 paralogs were first identified in the transcriptome and genome of Japanese flounder (Paralichthys olivaceus). We also explored SMAD3 duplication in other selected species. Following identification, genomic structure, phylogenetic reconstruction, and synteny analyses performed by MrBayes and online bioinformatic tools confirmed that smad3a/3b most likely originated from the teleost-specific WGD. Additionally, selection pressure analysis and expression pattern of the two genes performed by PAML and quantitative real-time PCR (qRT-PCR) revealed evidence of subfunctionalization of the two SMAD3 paralogs in teleost. Our results indicate that two SMAD3 genes originate from teleost-specific WGD, remain transcriptionally active, and may have likely undergone subfunctionalization. This study provides novel insights to the evolution fates of smad3a/3b and draws attentions to future function analysis of SMAD3 gene family.

Evolutionary history of the reprimo tumor suppressor gene family in vertebrates with a description of a new reprimo gene lineage

Genes related to human diseases should be natural targets for evolutionary studies, since they could provide clues regarding the genetic bases of pathologies and potential treatments. Here we studied the evolution of the reprimo gene family, a group of tumor-suppressor genes that are implicated in p53-mediated cell cycle arrest. These genes, especially the reprimo duplicate located on human chromosome 2, have been associated with epigenetic modifications correlated with transcriptional silencing and cancer progression. We demonstrate the presence of a third reprimo lineage that, together with the reprimo and reprimo-like genes, appears to have been differentially retained during the evolutionary history of vertebrates. We present evidence that these reprimo lineages originated early in vertebrate evolution and expanded as a result of the two rounds of whole genome duplications that occurred in the last common ancestor of vertebrates. The reprimo gene has been lost in birds, and the third reprimo gene lineage has been retained in only a few distantly related species, such as coelacanth and gar. Expression analyses revealed that the reprimo paralogs are mainly expressed in the nervous system. Different vertebrate lineages have retained different reprimo paralogs, and even in species that have retained multiple copies, only one of them is heavily expressed.

Selection on different genes with equivalent functions: the convergence story told by Hox genes along the evolution of aquatic mammalian lineages

Background

Convergent evolution has been a challenging topic for decades, being cetaceans, pinnipeds and sirenians textbook examples of three independent origins of equivalent phenotypes. These mammalian lineages acquired similar anatomical features correlated to an aquatic life, and remarkably differ from their terrestrial counterparts. Whether their molecular evolutionary history also involved similar genetic mechanisms underlying such morphological convergence nevertheless remained unknown. To test for the existence of convergent molecular signatures, we studied the molecular evolution of Hox genes in these three aquatic mammalian lineages, comparing their patterns to terrestrial mammals. Hox genes are transcription factors that play a pivotal role in specifying embryonic regional identity of nearly any bilateral animal, and are recognized major agents for diversification of body plans.

Results

We detected few signatures of positive selection on Hox genes across the three aquatic mammalian lineages and verified that purifying selection prevails in these sequences, as expected for pleiotropic genes. Genes found as being positively selected differ across the aquatic mammalian lineages, but we identified a substantial overlap of their developmental functions. Such pattern likely resides on the duplication history of Hox genes, which probably provided different possible evolutionary routes for achieving the same phenotypic solution.

Conclusions

Our results indicate that convergence occurred at a functional level of Hox genes along three independent origins of aquatic mammals. This conclusion reinforces the idea that different changes in developmental genes may lead to similar phenotypes, probably due to the redundancy provided by the participation of Hox paralogous genes in several developmental functions.

Electronic supplementary material

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

Hox genes and evolution

Hox proteins are a deeply conserved group of transcription factors originally defined for their critical roles in governing segmental identity along the antero-posterior (AP) axis in Drosophila. Over the last 30 years, numerous data generated in evolutionarily diverse taxa have clearly shown that changes in the expression patterns of these genes are closely associated with the regionalization of the AP axis, suggesting that Hox genes have played a critical role in the evolution of novel body plans within Bilateria. Despite this deep functional conservation and the importance of these genes in AP patterning, key questions remain regarding many aspects of Hox biology. In this commentary, we highlight recent reports that have provided novel insight into the origins of the mammalian Hox cluster, the role of Hox genes in the generation of a limbless body plan, and a novel putative mechanism in which Hox genes may encode specificity along the AP axis. Although the data discussed here offer a fresh perspective, it is clear that there is still much to learn about Hox biology and the roles it has played in the evolution of the Bilaterian body plan.

HOX gene complement and expression in the planarian Schmidtea mediterranea

Background

Freshwater planarians are well known for their regenerative abilities. Less well known is how planarians maintain spatial patterning in long-lived adult animals or how they re-pattern tissues during regeneration. HOX genes are good candidates to regulate planarian spatial patterning, yet the full complement or genomic clustering of planarian HOX genes has not yet been described, primarily because only a few have been detectable by in situ hybridization, and none have given morphological phenotypes when knocked down by RNAi.

Results

Because the planarian Schmidteamediterranea (S. mediterranea) is unsegmented, appendage less, and morphologically simple, it has been proposed that it may have a simplified HOX gene complement. Here, we argue against this hypothesis and show that S. mediterranea has a total of 13 HOX genes, which represent homologs to all major axial categories, and can be detected by whole-mount in situ hybridization using a highly sensitive method. In addition, we show that planarian HOX genes do not cluster in the genome, yet 5/13 have retained aspects of axially restricted expression. Finally, we confirm HOX gene axial expression by RNA deep-sequencing 6 anterior–posterior “zones” of the animal, which we provide as a dataset to the community to discover other axially restricted transcripts.

Conclusions

Freshwater planarians have an unappreciated HOX gene complexity, with all major axial categories represented. However, we conclude based on adult expression patterns that planarians have a derived body plan and their asexual lifestyle may have allowed for large changes in HOX expression from the last common ancestor between arthropods, flatworms, and vertebrates. Using our in situ method and axial zone RNAseq data, it should be possible to further understand the pathways that pattern the anterior–posterior axis of adult planarians.

Electronic supplementary material

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

The vertebrate Hox gene regulatory network for hindbrain segmentation: Evolution and diversification: Coupling of a Hox gene regulatory network to hindbrain segmentation is an ancient trait originating at the base of vertebrates

Hindbrain development is orchestrated by a vertebrate gene regulatory network that generates segmental patterning along the anterior-posterior axis via Hox genes. Here, we review analyses of vertebrate and invertebrate chordate models that inform upon the evolutionary origin and diversification of this network. Evidence from the sea lamprey reveals that the hindbrain regulatory network generates rhombomeric compartments with segmental Hox expression and an underlying Hox code. We infer that this basal feature was present in ancestral vertebrates and, as an evolutionarily constrained developmental state, is fundamentally important for patterning of the vertebrate hindbrain across diverse lineages. Despite the common ground plan, vertebrates exhibit neuroanatomical diversity in lineage-specific patterns, with different vertebrates revealing variations of Hox expression in the hindbrain that could underlie this diversification. Invertebrate chordates lack hindbrain segmentation but exhibit some conserved aspects of this network, with retinoic acid signaling playing a role in establishing nested domains of Hox expression.

Origin and evolution of GATA2a and GATA2b in teleosts: insights from tongue sole, Cynoglossus semilaevis

Background. Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication occurred in the lineage of teleost fish and is considered to be responsible for much of the biological diversification within the lineage. GATA2, a member of GATA family of transcription factors, is an important regulator of gene expression in hematopoietic cell in mammals, yet the role of this gene or its putative paralogs in ray-finned fishes remains relatively unknown.

Methods. In this study, we attempted to identify GATA2 sequences from the transcriptomes and genomes of multiple teleosts using the bioinformatic tools MrBayes, MEME, and PAML. Following identification, comparative analysis of genome structure, molecular evolution rate, and expression by real-time qPCR were used to predict functional divergence of GATA2 paralogs and their relative transcription in organs of female and male tongue soles (Cynoglossus semilaevis).

Results. Two teleost GATA2 genes were identified in the transcriptomes of tongue sole and Japanese flounder (Paralichthysolivaceus). Synteny and phylogenetic analysis confirmed that the two genes likely originated from the teleost-specific genome duplication . Additionally, selection pressure analysis predicted these gene duplicates to have undergone purifying selection and possible divergent new functions. This was supported by differential expression pattern of GATA2a and GATA2b observed in organs of female and male tongue soles.

Discussion. Our results indicate that two GATA2 genes originating from the first teleost-specific genome duplication have remained transcriptionally active in some fish species and have likely undergone neofunctionalization. This knowledge provides novel insights into the evolution of the teleost GATA2 genes and constituted important groundwork for further research on the GATA gene family.

HOX transcription factors are potential targets and markers in malignant mesothelioma

BACKGROUND

The HOX genes are a family of homeodomain-containing transcription factors that determine cellular identity during development and which are dys-regulated in some cancers. In this study we examined the expression and oncogenic function of HOX genes in mesothelioma, a cancer arising from the pleura or peritoneum which is associated with exposure to asbestos.

METHODS

We tested the sensitivity of the mesothelioma-derived lines MSTO-211H, NCI-H28, NCI-H2052, and NCI-H226 to HXR9, a peptide antagonist of HOX protein binding to its PBX co-factor. Apoptosis was measured using a FACS-based assay with Annexin, and HOX gene expression profiles were established using RT-QPCR on RNA extracted from cell lines and primary mesotheliomas. The in vivo efficacy of HXR9 was tested in a mouse MSTO-211H flank tumor xenograft model.

RESULTS

We show that HOX genes are significantly dysregulated in malignant mesothelioma. Targeting HOX genes with HXR9 caused apoptotic cell death in all of the mesothelioma-derived cell lines, and prevented the growth of mesothelioma tumors in a mouse xenograft model. Furthermore, the sensitivity of these lines to HXR9 correlated with the relative expression of HOX genes that have either an oncogenic or tumor suppressive function in cancer. The analysis of HOX expression in primary mesothelioma tumors indicated that these cells could also be sensitive to the disruption of HOX activity by HXR9, and that the expression of HOXB4 is strongly associated with overall survival.

CONCLUSION

HOX genes are a potential therapeutic target in mesothelioma, and HOXB4 expression correlates with overall survival.

Ancestral whole-genome duplication in the marine chelicerate horseshoe crabs

Whole-genome duplication (WGD) results in new genomic resources that can be exploited by evolution for rewiring genetic regulatory networks in organisms. In metazoans, WGD occurred before the last common ancestor of vertebrates, and has been postulated as a major evolutionary force that contributed to their speciation and diversification of morphological structures. Here, we have sequenced genomes from three of the four extant species of horseshoe crabs-Carcinoscorpius rotundicauda, Limulus polyphemus and Tachypleus tridentatus. Phylogenetic and sequence analyses of their Hox and other homeobox genes, which encode crucial transcription factors and have been used as indicators of WGD in animals, strongly suggests that WGD happened before the last common ancestor of these marine chelicerates >135 million years ago. Signatures of subfunctionalisation of paralogues of Hox genes are revealed in the appendages of two species of horseshoe crabs. Further, residual homeobox pseudogenes are observed in the three lineages. The existence of WGD in the horseshoe crabs, noted for relative morphological stasis over geological time, suggests that genomic diversity need not always be reflected phenotypically, in contrast to the suggested situation in vertebrates. This study provides evidence of ancient WGD in the ecdysozoan lineage, and reveals new opportunities for studying genomic and regulatory evolution after WGD in the Metazoa.

Distinct functions of two olfactory marker protein genes derived from teleost-specific whole genome duplication

BACKGROUND

Whole genome duplications (WGDs) have been proposed to have made a significant impact on vertebrate evolution. Two rounds of WGD (1R and 2R) occurred in the common ancestor of Gnathostomata and Cyclostomata, followed by the third-round WGD (3R) in a common ancestor of all modern teleosts. The 3R-derived paralogs are good models for understanding the evolution of genes after WGD, which have the potential to facilitate phenotypic diversification. However, the recent studies of 3R-derived paralogs tend to be based on in silico analyses. Here we analyzed the paralogs encoding teleost olfactory marker protein (OMP), which was shown to be specifically expressed in mature olfactory sensory neurons and is expected to be involved in olfactory transduction.

RESULTS

Our genome database search identified two OMPs (OMP1 and OMP2) in teleosts, whereas only one was present in other vertebrates. Phylogenetic and synteny analyses suggested that OMP1 and 2 were derived from 3R. Both OMPs showed distinct expression patterns in zebrafish; OMP1 was expressed in the deep layer of the olfactory epithelium (OE), which is consistent with previous studies of mice and zebrafish, whereas OMP2 was sporadically expressed in the superficial layer. Interestingly, OMP2 was expressed in a very restricted region of the retina as well as in the OE. In addition, the analysis of transcriptome data of spotted gar, a non-teleost fish, revealed that single OMP gene was expressed in the eyes.

CONCLUSION

We found distinct expression patterns of zebrafish OMP1 and 2 at the tissue and cellular level. These differences in expression patterns may be explained by subfunctionalization as the model of molecular evolution. Namely, single OMP gene was speculated to be originally expressed in the OE and the eyes in the common ancestor of all Osteichthyes (bony fish including tetrapods). Then, two OMP gene paralogs derived from 3R-WGD reduced and specialized the expression patterns. This study provides a good example for analyzing a functional subdivision of the teleost OE and eyes as revealed by 3R-derived paralogs of OMPs.

Molecular characterization, origin, and evolution of teleost p68 gene family: Insights from Japanese flounder, Paralichthys olivaceus

Two rounds of whole-genome duplication occurred in the common ancestor of vertebrates. Later, a third round genome duplication occurred in the teleost fishes. As a prototype member of DEAD-box RNA helicases, the function of p68 helicase in development has been well investigated in human, however, limited information is available regarding the regulatory function of this gene in the development of teleosts. In this study, being an important farmed fish in North China, Japanese flounder (Paralichthys olivaceus) was used as model fish to investigate the role of p68 gene in teleost development. Two p68 genes were first identified from Japanese flounder. Molecular characterization of them was performed by analyzing the exon-intron boundaries. Then, we confirmed that such two teleost p68 genes originated from teleost-specific genome duplication through phylogenetic and synteny analyses. Additionally, comparative analyses of amino acid sequences, variation in selective pressure, and expression profiles of p68 genes revealed probable sub-functionalization fate of teleost p68 genes after the duplication. Therefore, this study supplements the evolutionary properties of teleost p68 gene family and provides the groundwork for further studying the regulatory function of p68 genes in the development of teleosts.

Comparative Evolution of Duplicated Ddx3 Genes in Teleosts: Insights from Japanese Flounder, Paralichthys olivaceus

Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication event occurred in the stem lineage of ray-finned fishes. This teleost-specific genome duplication is thought to be responsible for the biological diversification of ray-finned fishes. DEAD-box polypeptide 3 (DDX3) belongs to the DEAD-box RNA helicase family. Although their functions in humans have been well studied, limited information is available regarding their function in teleosts. In this study, two teleost Ddx3 genes were first identified in the transcriptome of Japanese flounder (Paralichthys olivaceus). We confirmed that the two genes originated from teleost-specific genome duplication through synteny and phylogenetic analysis. Additionally, comparative analysis of genome structure, molecular evolution rate, and expression pattern of the two genes in Japanese flounder revealed evidence of subfunctionalization of the duplicated Ddx3 genes in teleosts. Thus, the results of this study reveal novel insights into the evolution of the teleost Ddx3 genes and constitute important groundwork for further research on this gene family.

What to compare and how: Comparative transcriptomics for Evo-Devo

Evolutionary developmental biology has grown historically from the capacity to relate patterns of evolution in anatomy to patterns of evolution of expression of specific genes, whether between very distantly related species, or very closely related species or populations. Scaling up such studies by taking advantage of modern transcriptomics brings promising improvements, allowing us to estimate the overall impact and molecular mechanisms of convergence, constraint or innovation in anatomy and development. But it also presents major challenges, including the computational definitions of anatomical homology and of organ function, the criteria for the comparison of developmental stages, the annotation of transcriptomics data to proper anatomical and developmental terms, and the statistical methods to compare transcriptomic data between species to highlight significant conservation or changes. In this article, we review these challenges, and the ongoing efforts to address them, which are emerging from bioinformatics work on ontologies, evolutionary statistics, and data curation, with a focus on their implementation in the context of the development of our database Bgee (http://bgee.org). J. Exp. Zool. (Mol. Dev. Evol.) 324B: 372–382, 2015. © 2015 The Authors. J. Exp. Zool. (Mol. Dev. Evol.) published by Wiley Periodicals, Inc.

Transcriptional dynamics of homeobox C11 gene in water buffalo bubalus bubalis

The Hox complex contains 39 genes clustered into four groups involved in cell differentiation and development. We cloned full-length sequence of Hoxc11 gene from water buffalo Bubalus bubalis, assessed its copy number, localized the same onto the chromosome 5, and studied its evolutionary conservation across the species. Northern hybridization of Hoxc11 showed a 2.2 kb band in the tissues analyzed. Real-Time PCR showed highest expression of Hoxc11 gene in lung followed by spleen, spermatozoa, and testis. Six interacting partners of this gene showed higher expression in spleen, lung, testis, and spermatozoa. During the early stages of development, Hoxc11 and its interacting partners both showed lower expression, which then became prominent during the age of 1-3 years, regressed drastically thereafter, and remained so until the animal's life time (∼ 20 years). The high expression of Hoxc11 and its interacting partners in spermatozoa and testis during the onset of puberty suggests its likely role in the differentiation of gonads and subsequent reproductive activities. Additional work on Hoxc11 especially, in the context of respiratory, immunological, and in/fertility in other species, including humans would be useful for establishing its broader biological significance towards the enrichment of functional and comparative genomics.

The timing of Timezyme diversification in vertebrates

All biological functions in vertebrates are synchronized with daily and seasonal changes in the environment by the time keeping hormone melatonin. Its nocturnal surge is primarily due to the rhythmic activity of the arylalkylamine N-acetyl transferase AANAT, which thus became the focus of many investigations regarding its evolution and function. Various vertebrate isoforms have been reported from cartilaginous fish to mammals but their origin has not been clearly established. Using phylogeny and synteny, we took advantage of the increasing number of available genomes in order to test whether the various rounds of vertebrate whole genome duplications were responsible for the diversification of AANAT. We highlight a gene secondary loss of the AANAT2 in the Sarcopterygii, revealing for the first time that the AAANAT1/2 duplication occurred before the divergence between Actinopterygii (bony fish) and Sarcopterygii (tetrapods, lobe-finned fish, and lungfish). We hypothesize the teleost-specific whole genome duplication (WDG) generated the appearance of the AANAT1a/1b and the AANAT2/2′paralogs, the 2′ isoform being rapidly lost in the teleost common ancestor (ray-finned fish). We also demonstrate the secondary loss of the AANAT1a in a Paracantopterygii (Atlantic cod) and of the 1b in some Ostariophysi (zebrafish and cave fish). Salmonids present an even more diverse set of AANATs that may be due to their specific WGD followed by secondary losses. We propose that vertebrate AANAT diversity resulted from 3 rounds of WGD followed by previously uncharacterized secondary losses. Extant isoforms show subfunctionalized localizations, enzyme activities and affinities that have increased with time since their emergence.

Molecular characterization and functional divergence of two Gadd45g homologs in sex determination in half-smooth tongue sole (Cynoglossus semilaevis)

The growth arrest and DNA-damage-inducible protein 45 gamma (Gadd45g) is known to play a major role in embryonic development and sex determination. In this study, two Gadd45g genes were isolated from half-smooth tongue sole (Cynoglossus semilaevis). Using chromosomal fluorescence in situ hybridization (FISH), Gadd45g1 and Gadd45g2 were located on the W and Z chromosomes, respectively. The full-length cDNA sequences of Gadd45g1 (1270bp) and Gadd45g2 (1181bp) were predicted to contain a 480-bp coding sequence that could encode a protein of 159 amino acids residues. A phylogenetic tree showed that the predicted Gadd45g1 and Gadd45g2 amino acid sequences clustered closely in one branch. It is proposed that Gadd45g1 and Gadd45g2 are paralogous genes derived from the divergence of the sex chromosome. Ka/Ks ratios indicated that Gadd45g1 and Gadd45g2 may have undergone a high number of mutations and have a divergence time of only about 68,000years, although Gadd45g homologs are highly conserved. The qRT-PCR demonstrated that Gadd45g1 and Gadd45g2 were highly expressed in ovary, and negligibly expressed in testis of male and neo-male. During development of the ovary (from 80 to 150days), the expression levels of both genes reached high levels. Gadd45g1 was also highly expressed at 50days, the stage just before gonad differentiation in C. semilaevis. All these findings imply functional divergence of the two Gadd45g homologs; Gadd45g1 may be necessary for sex differentiation in the early stage of gonad development, and then Gadd45g1 and Gadd45g2 maintain ovary development and the female character of half-smooth tongue sole.