Hypermutability of HoxA13A and functional divergence from its paralog are associated with the origin of a novel developmental feature in zebrafish and related taxa (Cypriniformes).

Gene duplication is widely regarded as the predominant mechanism by which genes with new functions and associated phenotypic novelties arise. A whole genome duplication occurred shortly before the most recent common ancestor of teleosts, the most diverse chordate group, resulting in duplication and retention of many Hox cluster genes. Because they play a key role in determination of body plan morphology, it has been widely assumed that Hox genes play a key role in the evolution of diverse metazoan body plans. However, it is not clear whether certain aspects of molecular evolution, such as asymmetric divergence and neofunctionalization, contribute to the initial retention of paralogs. We investigate the molecular evolution and functional divergence of the duplicated HoxA13 paralogs in zebrafish to determine when asymmetric divergence and functional divergence occurred after the duplication event. Our findings demonstrate the contribution of gene duplication to the evolution of novel features through evolutionary mechanisms other than those traditionally investigated, such as positive selection occurring immediately after gene duplication. Rather, we find a latent build up of molecular changes in a gene associated with the development of a novel feature in a very diverse group of fishes.

Adaptive evolution of Hox-gene homeodomains after cluster duplications.

BACKGROUND: Hox genes code for homeodomain-containing transcription factors that function in cell fate determination and embryonic development. Hox genes are arranged in clusters with up to 14 genes. This archetypical chordate cluster has duplicated several times in vertebrates, once at the origin of vertebrates and once at the origin of gnathostoms, an additional duplication event is associated with the origin of teleosts and the agnanths, suggesting that duplicated Hox cluster genes are involved in the genetic mechanisms behind the diversification of vertebrate body plans, and the origin of morphological novelties. Preservation of duplicate genes is promoted by functional divergence of paralogs, either by subfunction partitioning among paralogs or the acquisition of a novel function by one paralog. But for Hox genes the mechanisms of paralog divergence is unknown, leaving open the role of Hox gene duplication in morphological evolution. RESULTS: Here, we use several complementary methods, including branch-specific dN/dS ratio tests, branch-site dN/dS ratio tests, clade level amino acid conservation/variation patterns, and relative rate ratio tests, to show that the homeodomain of Hox genes was under positive Darwinian selection after cluster duplications. CONCLUSION: Our results suggest that positive selection acted on the homeodomain immediately after Hox clusters duplications. The location of sites under positive selection in the homeodomain suggests that they are involved in protein-protein interactions. These results further suggest that adaptive evolution actively contributed to Hox-gene homeodomain functions.

Repressor domain and nuclear localization signal of the murine Hoxa-11 protein are located in the homeodomain: no evidence for role of poly alanine stretches in transcriptional repression.

Hoxa-11 is a member of the homeodomain class of transcription factors, which play important roles in metazoan development. Hoxa-11 is particularly interesting because it is involved in a major mammalian innovation, uterus development and gestation. We are interested in the molecular changes underlying this evolutionary innovation. Although phenotypes resulting from loss of functions are well investigated (e.g., female sterility), little is known about the domains contributing to Hoxa-11 protein function. We therefore mapped the domains mediating two essential transcription factor functions, nuclear localization and transcriptional activity in the mouse Hoxa-11 protein. Our results show that the mammal-specific alanine repeat does not contribute to repressor activity, as has been hypothesized based on amino acid composition and analogy with other repressor domains. Interestingly, both the repressor domain as well as the nuclear localization signal (NLS) are located within the homeodomain, adding to the growing evidence that the homeodomain is a multifunctional domain which fulfills essential transcription factor functions beyond DNA binding. It is proposed that the high degree of conservation of the homeodomain is due to the multiple functional constraints that result from the various conserved functions accommodated in the homeodomain.

Adaptive evolution of HoxA-11 and HoxA-13 at the origin of the uterus in mammals.

The evolution of morphological characters is mediated by the evolution of developmental genes. Evolutionary changes can either affect cis-regulatory elements, leading to differences in their temporal and spatial regulation, or affect the coding region. Although there is ample evidence for the importance of cis-regulatory evolution, it has only recently been shown that transcription factors do not remain functionally equivalent during evolution. These results suggest that the evolution of transcription factors may play an active role in the evolution of development. To test this idea we investigated the molecular evolution of two genes essential for the development and function of the mammalian female reproductive organs, HoxA-11 and HoxA-13. We predicted that if coding-region evolution plays an active role in developmental evolution, then these genes should have experienced adaptive evolution at the origin of the mammalian female reproductive system. We report the sequences of HoxA-11 from basal mammalian and amniote taxa and analyse HoxA-11 and HoxA-13 for signatures of adaptive molecular evolution. The data demonstrate that these genes were under strong positive (directional) selection in the stem lineage of therian and eutherian mammals, coincident with the evolution of the uterus and vagina. These results support the idea that adaptive evolution of transcription factors can be an integral part in the evolution of novel structures.

GTF2IRD2 is located in the Williams-Beuren syndrome critical region 7q11.23 and encodes a protein with two TFII-I-like helix-loop-helix repeats.

Williams-Beuren syndrome (also known as Williams syndrome) is caused by a deletion of a 1.55- to 1.84-megabase region from chromosome band 7q11.23. GTF2IRD1 and GTF2I, located within this critical region, encode proteins of the TFII-I family with multiple helix-loop-helix domains known as I repeats. In the present work, we characterize a third member, GTF2IRD2, which has sequence and structural similarity to the GTF2I and GTF2IRD1 paralogs. The ORF encodes a protein with several features characteristic of regulatory factors, including two I repeats, two leucine zippers, and a single Cys-2/His-2 zinc finger. The genomic organization of human, baboon, rat, and mouse genes is well conserved. Our exon-by-exon comparison has revealed that GTF2IRD2 is more closely related to GTF2I than to GTF2IRD1 and apparently is derived from the GTF2I sequence. The comparison of GTF2I and GTF2IRD2 genes revealed two distinct regions of homology, indicating that the helix-loop-helix domain structure of the GTF2IRD2 gene has been generated by two independent genomic duplications. We speculate that GTF2I is derived from GTF2IRD1 as a result of local duplication and the further evolution of its structure was associated with its functional specialization. Comparison of genomic sequences surrounding GTF2IRD2 genes in mice and humans allows refinement of the centromeric breakpoint position of the primate-specific inversion within the Williams-Beuren syndrome critical region.