Homeobox gene clusters and the human paralogy map

Identification and validation of DLX4 as a prognostic and diagnostic biomarker for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is a lethal cancer, and biomarkers for exact diagnosis and predicting prognosis are urgently needed. The present study aimed to determine the roles of distal-less homeobox (DLX) family genes in ccRCC. The clinicopathological and mRNA expression data of patients with ccRCC were derived from The Cancer Genome Atlas database. Kaplan-Meier curves, univariate and multivariate Cox hazard analyses, in addition to receiver operator characteristic curves were used to evaluate the prognostic and diagnostic values. A single-sample gene set enrichment analysis was used to quantify the infiltration levels of immune cells. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry were conducted to examine the expression levels of DLX4 in tumor and adjacent tissue; the results demonstrated that DLX4 was highly expressed in ccRCC tissues compared with normal renal tissues. Furthermore, DLX4 expression was associated with tumor stage and grade. High proportions of males, advanced pathological stage, higher tumor grade and T, N and M stage were also observed in the high DLX4 expression group. Patients with the high DLX4 expression levels tended to have lower overall survival and disease-free survival rates compared with those with low DLX4 expression. DLX4 expression also showed favorable diagnostic efficiency in ccRCC patients. Based on functional enrichment analysis, cell cycle related pathways, epithelial-mesenchymal transition, glycolysis and inflammatory response were associated with the expression levels of DLX4. Furthermore, DLX4 expression was revealed to be associated with tumor immunosuppressive microenvironment. Overall, the expression level of DLX4 may be considered a novel prognostic indicator in ccRCC and a specific diagnostic biomarker for patients with ccRCC.

Deeply conserved synteny and the evolution of metazoan chromosomes

Animal genomes show networks of deeply conserved gene linkages whose phylogenetic scope and chromosomal context remain unclear. Here, we report chromosome-scale conservation of synteny among bilaterians, cnidarians, and sponges and use comparative analysis to reconstruct ancestral chromosomes across major animal groups. Comparisons among diverse metazoans reveal the processes of chromosome evolution that produced contemporary karyotypes from their Precambrian progenitors. On the basis of these findings, we introduce a simple algebraic representation of chromosomal change and use it to establish a unified systematic framework for metazoan chromosome evolution. We find that fusion-with-mixing, a previously unappreciated mode of chromosome change, has played a central role. We find that relicts of several metazoan chromosomal units are preserved in unicellular eukaryotes. These conserved pre-metazoan linkages include the chromosomal unit that encodes the most diverse set of metazoan homeobox genes, suggesting a candidate genomic context for the early diversification of this key gene family.

Hoxd10 Is Required Systemically for Secretory Activation in Lactation and Interacts Genetically with Hoxd9

Targeted disruption of the murine Hoxd10 gene (ΔHoxd10) leads to a high frequency of localized (gland-to-gland or regionally within a gland) lactation impairment in homozygous mutant mice as a single gene mutation. The effect of Hoxd10 disruption was enhanced by simultaneous disruption of Hoxd9 (ΔHoxd9/d10), a mutation shown previously to have no effect on mammary function as a single gene alteration. Mammary glands of homozygous ΔHoxd10 and ΔHoxd9/d10 females were indistinguishable from those of wild type littermate and age-matched control mice in late pregnancy. However, in lactation, 47% of homozygous ΔHoxd10 females, and 100% of homozygous ΔHoxd9/d10 females, showed localized or complete failure of two or more glands to undergo lactation-associated morphological changes and to secrete milk. Affected regions of ΔHoxd10 and ΔHoxd9/d10 mutants showed reduced prolactin receptor expression, reduced signal transducer and activator transcription protein 5 (STAT5) phosphorylation, reduced expression of downstream milk proteins, mislocalized glucose transporter 1 (GLUT1), increased STAT3 expression and phosphorylation, recruitment of leukocytes, altered cell cycle status, and increased apoptosis relative to unaffected regions and wild type control glands. Despite these local effects on alveolar function, transplantation results and hormone analysis indicate that Hoxd10 primarily has systemic functions that confer attenuated STAT5 phosphorylation on both wild type and ΔHoxd10 transplants when placed in ΔHoxd10 hosts, thereby exacerbating an underlying propensity for lactation failure in C57Bl/6 mice.

BP1, a potential biomarker for breast cancer prognosis

Homeobox genes are critical in tumor development. An isoform protein of DLX4 called BP1 is expressed in 80% of invasive ductal breast carcinomas. BP1 overexpression is implicated in an aggressive phenotype and poor prognosis. BP1 upregulation is associated with estrogen receptor negativity so those tumors do not respond to antiestrogens. Breast cancer is the second leading cause of death in women. BP1 could serve as both a novel prognostic biomarker for breast cancer and a therapeutic target. In this review, we address the role of BP1 protein in tumorigenesis of breast cancer and four other malignancies. A number of functions of BP1 in cancer are also discussed.

The Case of the Missing Ancient Fungal Polyploids

Polyploidy-the increase in the number of whole chromosome sets-is an important evolutionary force in eukaryotes. Polyploidy is well recognized throughout the evolutionary history of plants and animals, where several ancient events have been hypothesized to be drivers of major evolutionary radiations. However, fungi provide a striking contrast: while numerous recent polyploids have been documented, ancient fungal polyploidy is virtually unknown. We present a survey of known fungal polyploids that confirms the absence of ancient fungal polyploidy events. Three hypotheses may explain this finding. First, ancient fungal polyploids are indeed rare, with unique aspects of fungal biology providing similar benefits without genome duplication. Second, fungal polyploids are not successful in the long term, leading to few extant species derived from ancient polyploidy events. Third, ancient fungal polyploids are difficult to detect, causing the real contribution of polyploidy to fungal evolution to be underappreciated. We consider each of these hypotheses in turn and propose that failure to detect ancient events is the most likely reason for the lack of observed ancient fungal polyploids. We examine whether existing data can provide evidence for previously unrecognized ancient fungal polyploidy events but discover that current resources are too limited. We contend that establishing whether unrecognized ancient fungal polyploidy events exist is important to ascertain whether polyploidy has played a key role in the evolution of the extensive complexity and diversity observed in fungi today and, thus, whether polyploidy is a driver of evolutionary diversifications across eukaryotes. Therefore, we conclude by suggesting ways to test the hypothesis that there are unrecognized polyploidy events in the deep evolutionary history of the fungi.

Increased copy number of the DLX4 homeobox gene in breast axillary lymph node metastasis

DLX4 is a homeobox gene strongly implicated in breast tumor progression and invasion. Our main objective was to determine the DLX4 copy number status in sentinel lymph node (SLN) metastasis to assess its involvement in the initial stages of the axillary metastatic process. A total of 37 paired samples of SLN metastasis and primary breast tumors (PBT) were evaluated by fluorescence in situ hybridization, quantitative polymerase chain reaction and array comparative genomic hybridization assays. DLX4 increased copy number was observed in 21.6% of the PBT and 24.3% of the SLN metastasis; regression analysis demonstrated that the DLX4 alterations observed in the SLN metastasis were dependent on the ones in the PBT, indicating that they occur in the primary tumor cell populations and are maintained in the early axillary metastatic site. In addition, regression analysis demonstrated that DLX4 alterations (and other DLX and HOXB family members) occurred independently of the ones in the HER2/NEU gene, the main amplification driver on the 17q region. Additional studies evaluating DLX4 copy number in non-SLN axillary lymph nodes and/or distant breast cancer metastasis are necessary to determine if these alterations are carried on and maintained during more advanced stages of tumor progression and if could be used as a predictive marker for axillary involvement.

Evolutionary expression of glucose-dependent-insulinotropic polypeptide (GIP)

Glucose-dependent insulinotropic polypeptide (GIP) is a mammalian incretin hormone released into the circulation following nutrient ingestion. We examined the functional evolution of GIP and its relationship with insulin to delineate their respective roles in promoting nutrient efficiency. Expression patterns were examined in the sea lamprey (Petromyzon marinus), a basal vertebrate lacking a distinct pancreas, and in the zebrafish, Xenopus laevis, chicken, and mouse, organisms possessing extraintestinal pancreata. Although sea lamprey genomic analysis predicted a potential GIP-like gene, transcripts were not detected, and insulin expression was confined to the caudal pancreatic bud. GIP was detected in both the intestine and pancreas of the zebrafish and X. laevis. In contrast, GIP and insulin expression were limited to the intestine and pancreas, respectively, in chicken and mouse. Phylogenetic analysis of the glucagon-like ligands suggested proglucagon as the common ancestor, supporting the theory that GIP arose as a gene duplication of proglucagon. Insulin-secreting cells in the sea lamprey intestine may have obviated the need for an enteroinsular axis, and zebrafish may represent an evolutionary transition where GIP does not yet function as an incretin hormone. These observations are consistent with the hypothesis that GIP and insulin influence survival advantage by enhancing the efficiency of nutrient absorption and energy storage.

HOXA10 is overexpressed in human ovarian clear cell adenocarcinoma and correlates with poor survival

Human homeobox gene (HOX) A10 is a homeobox allotype gene of the HOXA family in the HOX family. Human homeobox gene A10 may play an important role in cancer development. However, the role of HOXA10 in the carcinogenesis of ovarian clear cell adenocarcinoma (OCCA) has not been established. We have evaluated the prognostic significance of HOXA10 expression for human OCCA and the effects of HOXA10 on proliferation, motility, and invasion of OCCA cells. We found that HOXA10 was not expressed in normal ovarian epithelium, ovarian endometrial cysts, and ovarian serous carcinomas, but 20 (68.9%) of the 29 OCCAs were positive for the expression of HOXA10. Human homeobox gene A10 expression was negatively correlated to the 5-year survival of OCCA patients (R = -0.442, P = 0.043). When a HOXA10 expression vector was stably transfected into a human OCCA cell line, ES-2, the proliferation rate of ES-2-HOXA10 was much higher than the vector control, the motility of ES-2-HOXA10 cells was significantly increased compared with the control (P < 0.05), and the invasion of ES-2-HOXA10 cells was also much higher than the vector control (P < 0.01 at 5 hours and 12 hours after scratching). In conclusion, HOXA10 was overexpressed in OCCA and was correlated with poor survival. HOXA10 promotes proliferation, migration, and invasion of OCCA cells. Human homeobox gene A10 could be a promising prognostic marker for OCCA.

Mutations of polycomb-associated gene ASXL1 in myelodysplastic syndromes and chronic myelomonocytic leukaemia

The myelodysplastic syndromes (MDSs) are a heterogeneous group of clonal haematological diseases characterized by ineffective haematopoiesis and predisposition to acute myeloid leukaemia (AML). The pathophysiology of MDSs remains unclear. A definition of the molecular biology of MDSs may lead to a better classification, new prognosis indicators and new treatments. We studied a series of 40 MDS/AML samples by high-density array-comparative genome hybridization (aCGH). The genome of MDSs displayed a few alterations that can point to candidate genes, which potentially regulate histone modifications and WNT pathways (e.g. ASXL1, ASXL2, UTX, CXXC4, CXXC5, TET2, TET3). To validate some of these candidates we studied the sequence of ASXL1. We found mutations in the ASXL1 gene in four out of 35 MDS patients (11%). To extend these results we searched for mutations of ASXL1 in a series of chronic myelomonocytic leukaemias, a disease classified as MDS/Myeloproliferative disorder, and found mutations in 17 out of 39 patients (43%). These results show that ASXL1 might play the role of a tumour suppressor in myeloid malignancies.

Domain duplication, divergence, and loss events in vertebrate Msx paralogs reveal phylogenomically informed disease markers

Background

Msx originated early in animal evolution and is implicated in human genetic disorders. To reconstruct the functional evolution of Msx and inform the study of human mutations, we analyzed the phylogeny and synteny of 46 metazoan Msx proteins and tracked the duplication, diversification and loss of conserved motifs.

Results

Vertebrate Msx sequences sort into distinct Msx1, Msx2 and Msx3 clades. The sister-group relationship between MSX1 and MSX2 reflects their derivation from the 4p/5q chromosomal paralogon, a derivative of the original "MetaHox" cluster. We demonstrate physical linkage between Msx and other MetaHox genes (Hmx, NK1, Emx) in a cnidarian. Seven conserved domains, including two Groucho repression domains (N- and C-terminal), were present in the ancestral Msx. In cnidarians, the Groucho domains are highly similar. In vertebrate Msx1, the N-terminal Groucho domain is conserved, while the C-terminal domain diverged substantially, implying a novel function. In vertebrate Msx2 and Msx3, the C-terminal domain was lost. MSX1 mutations associated with ectodermal dysplasia or orofacial clefting disorders map to conserved domains in a non-random fashion.

Conclusion

Msx originated from a MetaHox ancestor that also gave rise to Tlx, Demox, NK, and possibly EHGbox, Hox and ParaHox genes. Duplication, divergence or loss of domains played a central role in the functional evolution of Msx. Duplicated domains allow pleiotropically expressed proteins to evolve new functions without disrupting existing interaction networks. Human missense sequence variants reside within evolutionarily conserved domains, likely disrupting protein function. This phylogenomic evaluation of candidate disease markers will inform clinical and functional studies.

Multiple chromosomal rearrangements structured the ancestral vertebrate Hox-bearing protochromosomes

While the proposal that large-scale genome expansions occurred early in vertebrate evolution is widely accepted, the exact mechanisms of the expansion--such as a single or multiple rounds of whole genome duplication, bloc chromosome duplications, large-scale individual gene duplications, or some combination of these--is unclear. Gene families with a single invertebrate member but four vertebrate members, such as the Hox clusters, provided early support for Ohno's hypothesis that two rounds of genome duplication (the 2R-model) occurred in the stem lineage of extant vertebrates. However, despite extensive study, the duplication history of the Hox clusters has remained unclear, calling into question its usefulness in resolving the role of large-scale gene or genome duplications in early vertebrates. Here, we present a phylogenetic analysis of the vertebrate Hox clusters and several linked genes (the Hox "paralogon") and show that different phylogenies are obtained for Dlx and Col genes than for Hox and ErbB genes. We show that these results are robust to errors in phylogenetic inference and suggest that these competing phylogenies can be resolved if two chromosomal crossover events occurred in the ancestral vertebrate. These results resolve conflicting data on the order of Hox gene duplications and the role of genome duplication in vertebrate evolution and suggest that a period of genome reorganization occurred after genome duplications in early vertebrates.

Phylogenetic and chromosomal analyses of multiple gene families syntenic with vertebrate Hox clusters

Background

Ever since the theory about two rounds of genome duplication (2R) in the vertebrate lineage was proposed, the Hox gene clusters have served as the prime example of quadruplicate paralogy in mammalian genomes. In teleost fishes, the observation of additional Hox clusters absent in other vertebrate lineages suggested a third tetraploidization (3R). Because the Hox clusters occupy a quite limited part of each chromosome, and are special in having position-dependent regulation within the multi-gene cluster, studies of syntenic gene families are needed to determine the extent of the duplicated chromosome segments. We have analyzed in detail 14 gene families that are syntenic with the Hox clusters to see if their phylogenies are compatible with the Hox duplications and the 2R/3R scenario. Our starting point was the gene family for the NPY family of peptides located near the Hox clusters in the pufferfish Takifugu rubripes, the zebrafish Danio rerio, and human.

Results

Seven of the gene families have members on at least three of the human Hox chromosomes and two families are present on all four. Using both neighbor-joining and quartet-puzzling maximum likelihood methods we found that 13 families have a phylogeny that supports duplications coinciding with the Hox cluster duplications. One additional family also has a topology consistent with 2R but due to lack of urochordate or cephalocordate sequences the time window when these duplications could have occurred is wider. All but two gene families also show teleost-specific duplicates.

Conclusion

Based on this analysis we conclude that the Hox cluster duplications involved a large number of adjacent gene families, supporting expansion of these families in the 2R, as well as in the teleost 3R tetraploidization. The gene duplicates presumably provided raw material in early vertebrate evolution for neofunctionalization and subfunctionalization.

The early ANTP gene repertoire: insights from the placozoan genome

The evolution of ANTP genes in the Metazoa has been the subject of conflicting hypotheses derived from full or partial gene sequences and genomic organization in higher animals. Whole genome sequences have recently filled in some crucial gaps for the basal metazoan phyla Cnidaria and Porifera. Here we analyze the complete genome of Trichoplax adhaerens, representing the basal metazoan phylum Placozoa, for its set of ANTP class genes. The Trichoplax genome encodes representatives of Hox/ParaHox-like, NKL, and extended Hox genes. This repertoire possibly mirrors the condition of a hypothetical cnidarian-bilaterian ancestor. The evolution of the cnidarian and bilaterian ANTP gene repertoires can be deduced by a limited number of cis-duplications of NKL and “extended Hox” genes and the presence of a single ancestral “ProtoHox” gene.

Early vertebrate chromosome duplications and the evolution of the neuropeptide Y receptor gene regions

Background

One of the many gene families that expanded in early vertebrate evolution is the neuropeptide (NPY) receptor family of G-protein coupled receptors. Earlier work by our lab suggested that several of the NPY receptor genes found in extant vertebrates resulted from two genome duplications before the origin of jawed vertebrates (gnathostomes) and one additional genome duplication in the actinopterygian lineage, based on their location on chromosomes sharing several gene families. In this study we have investigated, in five vertebrate genomes, 45 gene families with members close to the NPY receptor genes in the compact genomes of the teleost fishes Tetraodon nigroviridis and Takifugu rubripes. These correspond to Homo sapiens chromosomes 4, 5, 8 and 10.

Results

Chromosome regions with conserved synteny were identified and confirmed by phylogenetic analyses in H. sapiens, M. musculus, D. rerio, T. rubripes and T. nigroviridis. 26 gene families, including the NPY receptor genes, (plus 3 described recently by other labs) showed a tree topology consistent with duplications in early vertebrate evolution and in the actinopterygian lineage, thereby supporting expansion through block duplications. Eight gene families had complications that precluded analysis (such as short sequence length or variable number of repeated domains) and another eight families did not support block duplications (because the paralogs in these families seem to have originated in another time window than the proposed genome duplication events). RT-PCR carried out with several tissues in T. rubripes revealed that all five NPY receptors were expressed in the brain and subtypes Y2, Y4 and Y8 were also expressed in peripheral organs.

Conclusion

We conclude that the phylogenetic analyses and chromosomal locations of these gene families support duplications of large blocks of genes or even entire chromosomes. Thus, these results are consistent with two early vertebrate tetraploidizations forming a paralogon comprising human chromosomes 4, 5, 8 and 10 and one teleost tetraploidization. The combination of positional and phylogenetic data further strengthens the identification of orthologs and paralogs in the NPY receptor family.

Expansion of signaling genes for adaptive immune system evolution in early vertebrates

Background

The adaptive immune system (AIS) of jawed vertebrates is a sophisticated system mediated by numerous genes in specialized cells. Phylogenetic analysis indicates that emergence of the AIS followed the occurrence of two rounds of whole-genome duplication (2R-WGD) in early vertebrates, but little direct evidence linking these two events is available.

Results

We examined the relationship between 2R-WGD and the gain of AIS-related functions by numerous genes. To analyze the evolution of the many genes related to signal transduction in the AIS (defined as AIS genes), we identified groups of genes (defined as AIS subfamilies) that included at least one human AIS gene, its paralogs (if any), and its Drosophila ortholog(s). Genomic mapping revealed that numerous pairs of AIS genes and their paralogs were part of paralogons – series of paralogous regions that derive from a common ancestor – throughout the human genome, indicating that the genes were retained as duplicates after 2R-WGD. Outgroup comparison analysis revealed that subfamilies in which human and fly genes shared a nervous system-related function were significantly enriched among AIS subfamilies, as compared with the overall incidence of shared nervous system-related functions among all subfamilies in bilaterians. This finding statistically supports the hypothesis that AIS-related signaling genes were ancestrally involved in the nervous system of urbilaterians.

Conclusion

The current results suggest that 2R-WGD played a major role in the duplication of many signaling genes, ancestrally used in nervous system development and function, that were later co-opted for new functions during evolution of the AIS.

Evolution after gene duplication: models, mechanisms, sequences, systems, and organisms

Gene duplication is postulated to have played a major role in the evolution of biological novelty. Here, gene duplication is examined across levels of biological organization in an attempt to create a unified picture of the mechanistic process by which gene duplication can have played a role in generating biodiversity. Neofunctionalization and subfunctionalization have been proposed as important processes driving the retention of duplicate genes. These models have foundations in population genetic theory, which is now being refined by explicit consideration of the structural constraints placed upon genes encoding proteins through physical chemistry. Further, such models can be examined in the context of comparative genomics, where an integration of gene-level evolution and species-level evolution allows an assessment of the frequency of duplication and the fate of duplicate genes. This process, of course, is dependent upon the biochemical role that duplicated genes play in biological systems, which is in turn dependent upon the mechanism of duplication: whole genome duplication involving a co-duplication of interacting partners vs. single gene duplication. Lastly, the role that these processes may have played in driving speciation is examined.

Ancient complexity of the non-Hox ANTP gene complement in the anthozoan Nematostella vectensis: implications for the evolution of the ANTP superclass

The origin and evolution of ANTP superclass genes has raised controversial discussions. While recent evidence suggests that a true Hox cluster emerged after the cnidarian bilaterian split, the origin of the ANTP superclass as a whole remains unclear. Based on analyses of bilaterian genomes, it seems very likely that clustering has once been a characteristic of all ANTP homeobox genes and that their ancestors have emerged through several series of cis-duplications from the same genomic region. Since the diploblastic Cnidaria possess orthologs of some non-Hox ANTP genes, at least some steps of the expansion of this hypothetical homeobox gene array must have occurred in the last common ancestor of both lineages--but it is unknown to what extent. By screening the unassembled Nematostella genome, we have identified unambiguous orthologs to almost all non-Hox ANTP genes which are present in Bilateria--with the exception of En, Tlx and (possibly) Vax. Furthermore, Nematostella possesses ANTP genes that are missing in some bilaterian lineages, like the rough gene or NK7. In addition, several ANTP homeobox gene families have been independently duplicated in Nematostella. We conclude that the last cnidarian/bilaterian ancestor already harboured the almost full complement of non-Hox ANTP genes before the Hox system evolved.

Pre-bilaterian origins of the Hox cluster and the Hox code: evidence from the sea anemone, Nematostella vectensis

BACKGROUND

Hox genes were critical to many morphological innovations of bilaterian animals. However, early Hox evolution remains obscure. Phylogenetic, developmental, and genomic analyses on the cnidarian sea anemone Nematostella vectensis challenge recent claims that the Hox code is a bilaterian invention and that no "true" Hox genes exist in the phylum Cnidaria.

METHODOLOGY/PRINCIPAL FINDINGS

Phylogenetic analyses of 18 Hox-related genes from Nematostella identify putative Hox1, Hox2, and Hox9+ genes. Statistical comparisons among competing hypotheses bolster these findings, including an explicit consideration of the gene losses implied by alternate topologies. In situ hybridization studies of 20 Hox-related genes reveal that multiple Hox genes are expressed in distinct regions along the primary body axis, supporting the existence of a pre-bilaterian Hox code. Additionally, several Hox genes are expressed in nested domains along the secondary body axis, suggesting a role in "dorsoventral" patterning.

CONCLUSIONS/SIGNIFICANCE

A cluster of anterior and posterior Hox genes, as well as ParaHox cluster of genes evolved prior to the cnidarian-bilaterian split. There is evidence to suggest that these clusters were formed from a series of tandem gene duplication events and played a role in patterning both the primary and secondary body axes in a bilaterally symmetrical common ancestor. Cnidarians and bilaterians shared a common ancestor some 570 to 700 million years ago, and as such, are derived from a common body plan. Our work reveals several conserved genetic components that are found in both of these diverse lineages. This finding is consistent with the hypothesis that a set of developmental rules established in the common ancestor of cnidarians and bilaterians is still at work today.

Ancient linkage of a POU class 6 and an anterior hox-like gene in cnidaria: implications for the evolution of homeobox genes

Linkage analyses in metazoan genomes suggest two ancestral arrays for the majority of homeobox genes. The related homeobox genes and chromosomal regions that are dispersed in extant species derived possibly from only two single common ancestor regions. One proposed ancestral array, designated as ANTP mega-array, contains most of the ANTP class homeobox genes; the second, named the contraHox super-paralogon, would consist of the classes PRD, POU, LIM, CUT, prospero, TALE and SIX. Here, we report the tight linkage of a POU class 6 gene to an anterior Hox-like gene in the hydrozoan Eleutheria dichotoma and discuss its possible significance for the evolution of homeobox genes. POU class 6 genes also seem to be ancestrally linked to the HoxC and A clusters in vertebrates, despite POU homeobox genes belonging to the contraHox paralogon. Hence, the much tighter linkage of a POU class 6 gene to an anterior Hox-like gene in a cnidarian is possibly the evolutionary echo of an ancestral genomic region from which most metazoan homeobox classes emerged.

Evolution of the vertebrate glucose-dependent insulinotropic polypeptide (GIP) gene

The glucose-dependent insulinotropic polypeptide (GIP) gene is believed to have originated from a gene duplication event very early in vertebrate evolution that also produced the proglucagon gene, yet so far GIP has only been described within mammals. Here we report the identification of GIP genes in chicken, frogs, and zebrafish. The chicken and frog genes are organized in a similar fashion to mammalian GIP genes and contain 6 exons and 5 introns in homologous locations. These genes can also potentially be proteolytically processed in identical patterns as observed in the mammalian sequences that would yield a GIP hormone that is only one amino shorter than the mammalian sequences due to the removal of an extra basic residue by carboxypeptidase E. The zebrafish GIP gene and precursor protein is shorter than other vertebrate GIP genes and is missing exon 5. The predicted zebrafish GIP hormone is also shorter, being only 31 amino acids in length. The zebrafish GIP hormone is similar in length to the proglucagon-derived peptide hormones, peptides encoded from the gene most closely related to GIP. We suggest that the structure of zebrafish GIP is more similar to the ancestral gene, and that tetrapod GIP has been extended. The mammalian GIP hormone has also undergone a period of rapid sequence evolution early in mammalian evolution. The discovery of a conserved GIP in diverse vertebrate suggests that it has an essential role in physiology in diverse vertebrates, although it may have only recently evolved a role as an incretin hormone.

A low diversity of ANTP class homeobox genes in Placozoa

Homeobox genes of the ANTP and PRD classes play important roles in body patterning of metazoans, and a large diversity of these genes have been described in bilaterian animals and cnidarians. Trichoplax adhaerens (Phylum Placozoa) is a small multicellular marine animal with one of the simplest body organizations of all metazoans, showing no symmetry and a small number of distinct cell types. Only two ANTP class genes have been described from Trichoplax: the Hox/ParaHox gene Trox-2 and a gene related to the Not family. Here we report an extensive screen for ANTP class genes in Trichoplax, leading to isolation of three additional ANTP class genes. These can be assigned to the Dlx, Mnx and Hmx gene families. Sequencing approximately 12-20 kb around each gene indicates that none are part of tight gene clusters, and in situ hybridization reveals that at least two have spatially restricted expression around the periphery of the animal. The low diversity of ANTP class genes isolated in Trichoplax can be reconciled with the low anatomical complexity of this animal, although the finding that these genes are assignable to recognized gene families is intriguing.

Hox genes and their candidate downstream targets in the developing central nervous system

1. Homeobox (Hox) genes were originally discovered in the fruit fly Drosophila, where they function through a conserved homeodomain as transcriptional regulators to control embryonic morphogenesis. Since then over 1000 homeodomain proteins have been identified in several species. In vertebrates, 39 Hox genes have been identified as homologs of the original Drosophila complex, and like their Drosophila counterparts they are organized within chromosomal clusters. Vertebrate Hox genes have also been shown to play a critical role in embryonic development as transcriptional regulators. 2. Both the Drosophila and vertebrate Hox genes have been shown to interact with various cofactors, such as the TALE homeodomain proteins, in recognition of consensus sequences within regulatory elements of their target genes. These protein-protein interactions are believed to contribute to enhancing the specificity of target gene recognition in a cell-type or tissue- dependent manner. The regulatory activity of a particular Hox protein on a specific regulatory element is highly variable and dependent on its interacting partners within the transcriptional complex. 3. In vertebrates, Hox genes display spatially restricted patterns of expression within the developing CNS, both along the anterioposterior and dorsoventral axis of the embryo. Their restricted gene expression is suggestive of a regulatory role in patterning of the CNS, as well as in cell specification. Determining the precise function of individual Hox genes in CNS morphogenesis through classical mutational analyses is complicated due to functional redundancy between Hox genes. 4. Understanding the precise mechanisms through which Hox genes mediate embryonic morphogenesis requires the identification of their downstream target genes. Although Hox genes have been implicated in the regulation of several pathways, few target genes have been shown to be under their direct regulatory control. Development of methodologies used for the isolation of target genes and for the analysis of putative targets will be beneficial in establishing the genetic pathways controlled by Hox factors. 5. Within the developing CNS various cell adhesion molecules and signaling molecules have been identified as candidate downstream target genes of Hox proteins. These targets play a role in processes such as cell migration and differentiation, and are implicated in contributing to neuronal processes such as plasticity and/or specification. Hence, Hox genes not only play a role in patterning of the CNS during early development, but may also contribute to cell specification and identity.

An evolutionary history of the FGF superfamily

Fibroblast growth factors (FGF) are associated with multiple developmental and metabolic processes in triploblasts, and perhaps also in diploblasts. The evolution of the FGF superfamily has accompanied the major morphological and functional innovations of metazoan species. The study of FGFs throughout species shows that the FGF superfamily can be subdivided in eight families in present-day organisms and has evolved through phases of gene duplications and gene losses. At least two major expansions of the superfamily can be recognized: a first expansion increased the number of FGFs from one or few archeo-FGFs to eight proto-FGFs, prototypic of the eight families. A second expansion, which took place during euchordate evolution, is associated with genome duplications. It increased the number of members in the families. Subsequent losses reduced that number to the present-day figures.

Two rounds of whole genome duplication in the ancestral vertebrate

The hypothesis that the relatively large and complex vertebrate genome was created by two ancient, whole genome duplications has been hotly debated, but remains unresolved. We reconstructed the evolutionary relationships of all gene families from the complete gene sets of a tunicate, fish, mouse, and human, and then determined when each gene duplicated relative to the evolutionary tree of the organisms. We confirmed the results of earlier studies that there remains little signal of these events in numbers of duplicated genes, gene tree topology, or the number of genes per multigene family. However, when we plotted the genomic map positions of only the subset of paralogous genes that were duplicated prior to the fish-tetrapod split, their global physical organization provides unmistakable evidence of two distinct genome duplication events early in vertebrate evolution indicated by clear patterns of four-way paralogous regions covering a large part of the human genome. Our results highlight the potential for these large-scale genomic events to have driven the evolutionary success of the vertebrate lineage.

Conserved synteny between the Ciona genome and human paralogons identifies large duplication events in the molecular evolution of the insulin-relaxin gene family

The aims of the study were to outline the sequence of events that gave rise to the vertebrate insulin-relaxin gene family and the chromosomal regions in which they reside. We analyzed the gene content surrounding the human insulin/relaxin genes with respect to what family they belonged to and if the duplication history of investigated families parallels the evolution of the insulin-relaxin family members. Markov Clustering and phylogenetic analysis were used to determine family identity. More than 15% of the genes belonged to families that have paralogs in the regions, defining two sets of quadruplicate paralogy regions. Thereby, the localization of insulin/relaxin genes in humans is in accordance with those regions on human chromosomes 1, 11, 12, 19q (insulin/insulin-like growth factors) and 1, 6p/15q, 9/5, 19p (insulin-like factors/relaxins) were formed during two genome duplications. We compared the human genome with that of Ciona intestinalis, a species that split from the vertebrate lineage before the two suggested genome duplications. Two insulin-like orthologs were discovered in addition to the already described Ci-insulin gene. Conserved synteny between the Ciona regions hosting the insulin-like genes and the two sets of human paralogons implies their common origin. Linkage of the two human paralogons, as seen in human chromosome 1, as well as the two regions hosting the Ciona insulin-like genes suggests that a segmental duplication gave rise to the region prior to the genome doublings. Thus, preserved gene content provides support that genome duplication(s) in addition to segmental and single-gene duplications shaped the genomes of extant vertebrates.

A pair as a minimum: the two fibroblast growth factors of the nematode Caenorhabditis elegans

Fibroblast growth factors (FGFs) regulate many important developmental and homeostatic physiological events. The FGF superfamily contains several families. In this review, we present recent findings on the two FGFs of the nematode Caenorhabditis elegans from both functional and phylogenic points of view. C. elegans has a single FGFR (EGL-15) with two functionally exclusive isoforms, and two FGFs (LET-756 and EGL-17), which play distinct roles: an essential function for the former, and guidance of the migrating sex myoblasts for the latter. Regulation of homeostasis by control of the fluid balance could be the basis for the essential function of LET-756. Phylogenetic and functional studies suggest that LET-756, like vertebrate FGF9, -16, and -20, belongs to the FGF9 family, whereas EGL-17, like vertebrate FGF8, -17, and -18, could be included in the FGF8 family.

The genomic environment around the Aromatase gene: evolutionary insights

Background

The cytochrome P450 aromatase (CYP19), catalyses the aromatisation of androgens to estrogens, a key mechanism in vertebrate reproductive physiology. A current evolutionary hypothesis suggests that CYP19 gene arose at the origin of vertebrates, given that it has not been found outside this clade. The human CYP19 gene is located in one of the proposed MHC-paralogon regions (HSA15q). At present it is unclear whether this genomic location is ancestral (which would suggest an invertebrate origin for CYP19) or derived (genomic location with no evolutionary meaning). The distinction between these possibilities should help to clarify the timing of the CYP19 emergence and which taxa should be investigated.

Results

Here we determine the "genomic environment" around CYP19 in three vertebrate species Homo sapiens, Tetraodon nigroviridis and Xenopus tropicalis. Paralogy studies and phylogenetic analysis of six gene families suggests that the CYP19 gene region was structured through "en bloc" genomic duplication (as part of the MHC-paralogon formation). Four gene families have specifically duplicated in the vertebrate lineage. Moreover, the mapping location of the different paralogues is consistent with a model of "en bloc" duplication. Furthermore, we also determine that this region has retained the same gene content since the divergence of Actinopterygii and Tetrapods. A single inversion in gene order has taken place, probably in the mammalian lineage. Finally, we describe the first invertebrate CYP19 sequence, from Branchiostoma floridae.

Conclusion

Contrary to previous suggestions, our data indicates an invertebrate origin for the aromatase gene, given the striking conservation pattern in both gene order and gene content, and the presence of aromatase in amphioxus. We propose that CYP19 duplicated in the vertebrate lineage to yield four paralogues, followed by the subsequent loss of all but one gene in vertebrate evolution. Finally, we suggest that agnathans and lophotrocozoan protostomes should be investigated for the presence of aromatase.

Inhibition of DLX4 promotes apoptosis in choriocarcinoma cell lines

Homeodomain (HDM) proteins encoded by homeobox (HBX) genes represent a large family of transcriptional factors that control differentiation and development in certain cell types. DLX4 is a member of Distal-less (DLX) family of HBX genes. Recent studies have demonstrated that abnormal expression of DLX4 is present in several types of human tumors, such as breast cancer, leukemia and colon cancer. In the present study, we investigated DLX4 mRNA and protein expression in both normal placental tissues and human choriocarcinoma cell lines. Also, using RNA interference (RNAi) technique, we knocked down the expression of DLX4 and examined apoptosis in JEG-3 cells. Our studies demonstrated that DLX4 RNAi inhibited DLX4 mRNA expression and decreased DLX4 protein mass specifically and effectively, potentially enhancing apoptosis. Moreover, we examined expression of caspase-3 and caspase-8, and found that both caspases were increased after DLX4 knockdown. However, DLX4 RNAi did not influence Bax expression in JEG-3 cells. In conclusion, this study suggests that DLX4 may be involved in the survival of human choriocarcinoma cells, which may be mediated by the inhibition of apoptosis. The detailed mechanism needs further investigation.

Extensive duplications of phototransduction genes in early vertebrate evolution correlate with block (chromosome) duplications

Many gene families in mammals have members that are expressed more or less uniquely in the retina or differentially in specific retinal cell types. We describe here analyses of nine such gene families with regard to phylogenetic relationships and chromosomal location. The families are opsins, G proteins (alpha, beta, and gamma subunits), phosphodiesterases type 6, cyclic nucleotide-gated channels, G-protein-coupled receptor kinases, arrestins, and recoverins. The results suggest that multiple new gene copies arose in all of these families very early in vertebrate evolution during a period with extensive gene duplications. Many of the new genes arose through duplications of large chromosome regions (blocks of genes) or even entire chromosomes, as shown by linkage with other gene families. Some of the phototransduction families belong to the same duplicated regions and were thus duplicated simultaneously. We conclude that gene duplications in early vertebrate evolution probably helped facilitate the specialization of the retina and the subspecialization of different retinal cell types.

Phylogenetic analysis of Ciona intestinalis gene superfamilies supports the hypothesis of successive gene expansions

Understanding the formation of metazoan multigene families is a good approach to reconstitute the evolution of the chordate genome. In this attempt, the analysis of the genome of selected species provides valuable information. Ciona intestinalis belongs to the urochordates, whose lineage separated from the chordate lineage that later gave birth to vertebrates. We have searched available sequences from the small marine ascidian C. intestinalis for orthologs of members of five vertebrate superfamilies, including tyrosine kinase receptors, ETS, FOX and SOX transcription factors, and WNT secreted regulatory factors, and conducted phylogenetic analyses. We have found that most vertebrate subfamilies have a single C. intestinalis ortholog. Our results support the hypothesis of a gene expansion prior the base of chordate ancestry followed by another gene expansion during vertebrate evolution. They also indicate that Ciona intestinalis genome will be a very valuable tool for evolutionary analyses.

Molecular evolution of NPY receptor subtypes

The neuropeptide Y (NPY) system consists in mammals of three peptides and 4-5 G-protein-coupled receptors called Y receptors that are involved in a variety of physiological functions such as appetite regulation, circadian rhythm and anxiety. Both the receptor family and the peptide family display unexpected evolutionary complexity and flexibility as shown by information from different classes of vertebrates. The vertebrate ancestor most likely had a single peptide gene and three Y receptor genes, the progenitors of the Y1, Y2 and Y5 subfamilies. The receptor genes were probably located in the same chromosomal segment. Additional gene copies arose through the chromosome quadruplication that took place before the emergence of jawed vertebrates (gnathostomes) whereupon differential losses of the gene copies ensued. The inferred ancestral gnathostome gene repertoire most likely consisted of two peptide genes, NPY and PYY, and no less than seven Y receptor genes: four Y1-like (Y1, Y4/a, Y6, and Yb), two Y2-like (Y2 and Y7), and a single Y5 gene. Whereas additional peptide genes have arisen in various lineages, the most common trend among the Y receptor genes has been further losses. Mammals have lost Yb and Y7 (the latter still exists in frogs) and Y6 is a pseudogene in several mammalian species but appears to be still functional in some. One challenge is to find out if mammals have been deprived of any functions through these gene losses. Teleost fishes like zebrafish and pufferfish, on the other hand, have lost the two major appetite-stimulating receptors Y1 and Y5. Nevertheless, teleost fishes seem to respond to NPY with increased feeding why some other subtype probably mediates this effect. Another challenge is to deduce how Y2 and Y4 came to evolve an inhibitory effect on appetite. Changes in anatomical distribution of receptor expression may have played an important part in such functional switching along with changes in receptor structures and ligand preferences.

Natural history and functional divergence of protein tyrosine kinases

Cellular signaling is important for many biological processes including growth, differentiation, adhesion, motility and apoptosis. The protein tyrosine kinase (PTK) supergene family is the key mediator in cellular signaling in metazoans, directly associated with a variety of human diseases. All PTKs contain a highly conserved catalytic kinase domain, in spite of variable multi-domain structures. Within each PTK gene family, members exhibit functional divergence in substrate-specificity or temporal/tissue-specific expression, although their primary function is conserved. After conducting phylogenetic analysis on major PTK gene families, we found that the expanding of each PTK family was likely caused by gene or genome duplication event(s) that occurred before the emergence of teleosts but after the vertebrate-amphioxus split. We further investigated the evolutionary pattern of functional divergence after gene duplication in those gene families. Our results show that site-specific shifted evolutionary rate (altered functional constraint) is a common pattern in PTK gene family evolution.

Retinoic acid regulates the expression of PBX1, PBX2, and PBX3 in P19 cells both transcriptionally and post-translationally

Pre-B cell leukemia transcription factors (PBXs) are important co-factors for the transcriptional regulation mediated by a number of Hox proteins during embryonic development. It was previously shown that the expression of several Pbx genes is elevated in mouse embryo limb buds and embryonal carcinoma P19 cells upon retinoic acid (RA) treatment although the mechanism of this induction is not well understood. In this report, we demonstrate that PBX1a, PBX1b, PBX2, and PBX3 mRNAs and PBX1/2/3 proteins are induced during endodermal and neuronal differentiation of P19 cells in a RAR-dependent subtype-unspecific manner following RA treatment. The increases in both PBX1 mRNA and PBX3 mRNA levels are secondary responses to RA treatment requiring new proteins synthesis while the increase in PBX2 mRNA is a primary response. The RA-dependent increases in PBX1 mRNA, PBX2 mRNA, and PBX3 mRNA levels are likely to be transcriptionally regulated since the stability of these mRNAs does not change. In addition, the half-lives of PBX1/2/3 proteins are significantly extended by RA treatment. Two possible mechanisms could contribute to the stabilization of PBX proteins: PBX proteins associate with RA-dependent increased levels of MEIS proteins, and RA may decrease the proteasome dependent degradation of PBX proteins.

ParaDB: a tool for paralogy mapping in vertebrate genomes

We present ParaDB (http://abi.marseille.inserm.fr/paradb/), a new database for large-scale paralogy studies in vertebrate genomes. We intended to collect all information (sequence, mapping and phylogenetic data) needed to map and detect new paralogous regions, previously defined as Paralogons. The AceDB database software was used to generate graphical objects and to organize data. General data were automatically collated from public sources (Ensembl, GadFly and RefSeq). ParaDB provides access to data derived from whole genome sequences (Homo sapiens, Mus musculus and Drosophila melanogaster): cDNA and protein sequences, positional information, bibliographical links. In addition, we provide BLAST results for each protein sequence, InParanoid orthologs and 'In-Paralogs' data, previously established paralogy data, and, to compare vertebrates and Drosophila, orthology data.

Correlation of expression of BP1, a homeobox gene, with estrogen receptor status in breast cancer

BACKGROUND

BP1 is a novel homeobox gene cloned in our laboratory. Our previous studies in leukemia demonstrated that BP1 has oncogenic properties, including as a modulator of cell survival. Here BP1 expression was examined in breast cancer, and the relationship between BP1 expression and clinicopathological data was determined.

METHODS

Total RNA was isolated from cell lines, tumors, and matched normal adjacent tissue or tissue from autopsy. Reverse transcription polymerase chain reaction was performed to evaluate BP1 expression. Statistical analysis was accomplished with SAS.

RESULTS

Analysis of 46 invasive ductal breast tumors demonstrated BP1 expression in 80% of them, compared with a lack of expression in six normal breast tissues and low-level expression in one normal breast tissue. Remarkably, 100% of tumors that were negative for the estrogen receptor (ER) were BP1-positive, whereas 73% of ER-positive tumors expressed BP1 (P = 0.03). BP1 expression was also associated with race: 89% of the tumors of African American women were BP1-positive, whereas 57% of those from Caucasian women expressed BP1 (P = 0.04). However, there was no significant difference in BP1 expression between grades I, II, and III tumors. Interestingly, BP1 mRNA expression was correlated with the ability of malignant cell lines to cause breast cancer in mice.

CONCLUSION

Because BP1 is expressed abnormally in breast tumors, it could provide a useful target for therapy, particularly in patients with ER-negative tumors. The frequent expression of BP1 in all tumor grades suggests that activation of BP1 is an early event.

The human Hox-bearing chromosome regions did arise by block or chromosome (or even genome) duplications

Many chromosome regions in the human genome exist in four similar copies, suggesting that the entire genome was duplicated twice in early vertebrate evolution, a concept called the 2R hypothesis. Forty-two gene families on the four Hox-bearing chromosomes were recently analyzed by others, and 32 of these were reported to have evolutionary histories incompatible with duplications concomitant with the Hox clusters, thereby contradicting the 2R hypothesis. However, we show here that nine of the families have probably been translocated to the Hox-bearing chromosomes more recently, and that three of these belong to other chromosome quartets where they actually support the 2R hypothesis. We consider 13 families too complex to shed light on the chromosome duplication hypothesis. Among the remaining 20 families, 14 display phylogenies that support or are at least consistent with the Hox-cluster duplications. Only six families seem to have other phylogenies, but these trees are highly uncertain due to shortage of sequence information. We conclude that all relevant and analyzable families support or are consistent with block/chromosome duplications and that none clearly contradicts the 2R hypothesis.

Interaction between two ubiquitin-protein isopeptide ligases of different classes, CBLC and AIP4/ITCH

In metazoans, CBL proteins are RING finger type ubiquitin-protein isopeptide (E3) ligases involved in the down-regulation of epidermal growth factor tyrosine kinase receptors (EGFR). Among the three CBL proteins described in humans, CBLC (CBL3) remains poorly studied. By screening in parallel a human and a Caenorhabditis elegans library using the two-hybrid procedure in yeast, we found a novel interaction between Hsa-CBLC and Hsa-AIP4 or its C. elegans counterpart Cel-WWP1. Hsa-AIP4 and Cel-WWP1 are also ubiquitin E3 ligases. They contain a HECT (homologous to E6-AP C terminus) catalytic domain and four WW domains known to bind proline-rich regions. We confirmed the interaction between Hsa-CBLC and Hsa-AIP4 by a combination of glutathione S-transferase pull-down, co-immunoprecipitation, and colocalization experiments. We show that these two E3 ligases are involved in EGFR signaling because both become phosphorylated on tyrosine following epidermal growth factor stimulation. In addition, we observed that CBLC increases the ubiquitination of EGFR, and that coexpressing the WW domains of AIP4 exerts a dominant negative effect on EGFR ubiquitination. Finally, coexpressing CBLC and AIP4 induces a down-regulation of EGFR signaling. In conclusion, our data demonstrate that two E3 ligases of different classes can interact and cooperate to down-regulate EGFR signaling.

Carcinogenesis and translational controls: TACC1 is down-regulated in human cancers and associates with mRNA regulators

The three human TACC genes encode a family of proteins that are suspected to play a role in carcinogenesis. Their function is not precisely known; a Xenopus TACC protein called Maskin is involved in translational control, while the Drosophila D-TACC associates with microtubules and centrosomes. We have characterized the human TACC1 gene and its products. The TACC1 gene is located in region p12 of chromosome 8; its mRNA is ubiquitously expressed and encodes a protein with an apparent molecular mass of 125 kDa, which is cytoplasmic and mainly perinuclear. We show that TACC1 mRNA gene expression is down-regulated in various types of tumors. Using immunohistochemistry of tumor tissue-microarrays and sections, we confirm that the level of TACC1 protein is down-regulated in breast cancer. Finally, using the two-hybrid screen in yeast, GST pull-downs and co-immunoprecipitations, we identified two potential binding partners for TACC1, LSM7 and SmG. They constitute a conserved subfamily of Sm-like small proteins that associate with U6 snRNPs and play a role in several aspects of mRNA processing. We speculate that down-regulation of TACC1 may alter the control of mRNA homeostasis in polarized cells and participates in the oncogenic processes.

ORMDL proteins are a conserved new family of endoplasmic reticulum membrane proteins

BACKGROUND

Annotations of completely sequenced genomes reveal that nearly half of the genes identified are of unknown function, and that some belong to uncharacterized gene families. To help resolve such issues, information can be obtained from the comparative analysis of homologous genes in model organisms.

RESULTS

While characterizing genes from the retinitis pigmentosa locus RP26 at 2q31-q33, we have identified a new gene, ORMDL1, that belongs to a novel gene family comprising three genes in humans (ORMDL1, ORMDL2 and ORMDL3), and homologs in yeast, microsporidia, plants, Drosophila, urochordates and vertebrates. The human genes are expressed ubiquitously in adult and fetal tissues. The Drosophila ORMDL homolog is also expressed throughout embryonic and larval stages, particularly in ectodermally derived tissues. The ORMDL genes encode transmembrane proteins anchored in the endoplasmic reticulum (ER). Double knockout of the two Saccharomyces cerevisiae homologs leads to decreased growth rate and greater sensitivity to tunicamycin and dithiothreitol. Yeast mutants can be rescued by human ORMDL homologs.

CONCLUSIONS

From protein sequence comparisons we have defined a novel gene family, not previously recognized because of the absence of a characterized functional signature. The sequence conservation of this family from yeast to vertebrates, the maintenance of duplicate copies in different lineages, the ubiquitous pattern of expression in human and Drosophila, the partial functional redundancy of the yeast homologs and phenotypic rescue by the human homologs, strongly support functional conservation. Subcellular localization and the response of yeast mutants to specific agents point to the involvement of ORMDL in protein folding in the ER.

Ancient duplications of the human proglucagon gene

The human proglucagon gene (GCG) is encoded within a finished 576-kb DNA sequence generated by the Human Genome Project. GCG is flanked by 18 kb and 65 kb of DNA, 5' and 3', respectively, that do not encode genes. The genomic sequence that includes GCG was found to have a long history of gene duplication events. Some members of the glucagon-like family of genes, GCG on chromosome 2 and GIP on chromosome 17, may be products of ancient genome duplications on the early vertebrate lineage. A large genomic tandem duplication event that included DPP4-like and GCG genes occurred before the amphibian-mammal divergence, but one of the duplicated copies of GCG has been lost on the human lineage. Recently, a processed pseudogene of the X-chromosome-linked gene TIMM8A was inserted downstream of GCG. Some ancient duplicates of GCG may retain physiological functions in other vertebrates.

Coparalogy: physical and functional clusterings in the human genome

Two rounds of large-scale duplications are thought to have occurred in early vertebrate ancestry; this is now known as the "2R hypothesis." They have led to the constitution of subfamilies of paralogous genes. Chromosomal regions that contain present-day paralogs (paralogous regions or paralogons) have been identified in mammals. We show that sets of paralogons (PGs) can be assembled in a tentative "human genome paralogy map" that includes all autosomes and X. A total of 14 PGs, containing more than 1600 genes, were assembled in this paralogy map. Genes that belong to the same PG are coparalogs. We show that identification of coparalogy can be used (i) to broaden data on gene mapping, (ii) to identify physical gene clusters that derive from early cis-duplications, and (iii) to speculate on coevolution and coregulation of genes sharing a common structure or function (functional clusters). Thus, coparalogy analyses should parallel phylogenetic analyses and can help draw hypotheses on gene and genome evolution.