Convergent evolution of clustering of Iroquois homeobox genes across metazoans

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.

Iroquois Family Genes in Gastric Carcinogenesis: A Comprehensive Review

Gastric cancer (GC) is the fifth leading cause of cancer-associated death worldwide, accounting for 768,793 related deaths and 1,089,103 new cases in 2020. Despite diagnostic advances, GC is often detected in late stages. Through a systematic literature search, this study focuses on the associations between the Iroquois gene family and GC. Accumulating evidence indicates that Iroquois genes are involved in the regulation of various physiological and pathological processes, including cancer. To date, information about Iroquois genes in GC is very limited. In recent years, the expression and function of Iroquois genes examined in different models have suggested that they play important roles in cell and cancer biology, since they were identified to be related to important signaling pathways, such as wingless, hedgehog, mitogen-activated proteins, fibroblast growth factor, TGFβ, and the PI3K/Akt and NF-kB pathways. In cancer, depending on the tumor, Iroquois genes can act as oncogenes or tumor suppressor genes. However, in GC, they seem to mostly act as tumor suppressor genes and can be regulated by several mechanisms, including methylation, microRNAs and important GC-related pathogens. In this review, we provide an up-to-date review of the current knowledge regarding Iroquois family genes in GC.

Genome of the ramshorn snail Biomphalaria straminea-an obligate intermediate host of schistosomiasis

BACKGROUND

Schistosomiasis, or bilharzia, is a parasitic disease caused by trematode flatworms of the genus Schistosoma. Infection by Schistosoma mansoni in humans results when cercariae emerge into water from freshwater snails in the genus Biomphalaria and seek out and penetrate human skin. The snail Biomphalaria straminea is native to South America and is now also present in Central America and China, and represents a potential vector host for spreading schistosomiasis. To date, genomic information for the genus is restricted to the neotropical species Biomphalaria glabrata. This limits understanding of the biology and management of other schistosomiasis vectors, such as B. straminea.

FINDINGS

Using a combination of Illumina short-read, 10X Genomics linked-read, and Hi-C sequencing data, our 1.005 Gb B. straminea genome assembly is of high contiguity, with a scaffold N50 of 25.3 Mb. Transcriptomes from adults were also obtained. Developmental homeobox genes, hormonal genes, and stress-response genes were identified, and repeat content was annotated (40.68% of genomic content). Comparisons with other mollusc genomes (including Gastropoda, Bivalvia, and Cephalopoda) revealed syntenic conservation, patterns of homeobox gene linkage indicative of evolutionary changes to gene clusters, expansion of heat shock protein genes, and the presence of sesquiterpenoid and cholesterol metabolic pathway genes in Gastropoda. In addition, hormone treatment together with RT-qPCR assay reveal a sesquiterpenoid hormone responsive system in B. straminea, illustrating that this renowned insect hormonal system is also present in the lophotrochozoan lineage.

CONCLUSION

This study provides the first genome assembly for the snail B. straminea and offers an unprecedented opportunity to address a variety of phenomena related to snail vectors of schistosomiasis, as well as evolutionary and genomics questions related to molluscs more widely.

A duplication on chromosome 16q12 affecting the IRXB gene cluster is associated with autosomal dominant cone dystrophy with early tritanopic color vision defect

Cone dystrophies are a rare subgroup of inherited retinal dystrophies and hallmarked by color vision defects, low or decreasing visual acuity and central vision loss, nystagmus and photophobia. Applying genome-wide linkage analysis and array comparative genome hybridization, we identified a locus for autosomal dominant cone dystrophy on chromosome 16q12 in four independent multigeneration families. The locus is defined by duplications of variable size with a smallest region of overlap of 608 kb affecting the IRXB gene cluster and encompasses the genes IRX5 and IRX6. IRX5 and IRX6 belong to the Iroquois (Iro) protein family of homeodomain-containing transcription factors involved in patterning and regionalization of embryonic tissue in vertebrates, including the eye and the retina. All patients presented with a unique progressive cone dystrophy phenotype hallmarked by early tritanopic color vision defects. We propose that the disease underlies a misregulation of the IRXB gene cluster on chromosome 16q12 and demonstrate that overexpression of Irx5a and Irx6a, the two orthologous genes in zebrafish, results in visual impairment in 5-day-old zebrafish larvae.

Expression patterns of Irx genes in the developing chick inner ear

The vertebrate inner ear is a complex three-dimensional sensorial structure with auditory and vestibular functions. The molecular patterning of the developing otic epithelium creates various positional identities, consequently leading to the stereotyped specification of each neurosensory and non-sensory element of the membranous labyrinth. The Iroquois (Iro/Irx) genes, clustered in two groups (A: Irx1, Irx2, and Irx4; and B: Irx3, Irx5, and Irx6), encode for transcriptional factors involved directly in numerous patterning processes of embryonic tissues in many phyla. This work presents a detailed study of the expression patterns of these six Irx genes during chick inner ear development, paying particular attention to the axial specification of the otic anlagen. The Irx genes seem to play different roles at different embryonic periods. At the otic vesicle stage (HH18), all the genes of each cluster are expressed identically. Both clusters A and B seem involved in the specification of the lateral and posterior portions of the otic anlagen. Cluster B seems to regulate a larger area than cluster A, including the presumptive territory of the endolymphatic apparatus. Both clusters seem also to be involved in neurogenic events. At stages HH24/25-HH27, combinations of IrxA and IrxB genes participate in the specification of most sensory patches and some non-sensory components of the otic epithelium. At stage HH34, the six Irx genes show divergent patterns of expression, leading to the final specification of the membranous labyrinth, as well as to cell differentiation.

Favorable genomic environments for cis-regulatory evolution: A novel theoretical framework

Cis-regulatory changes are arguably the primary evolutionary source of animal morphological diversity. With the recent explosion of genome-wide comparisons of the cis-regulatory content in different animal species is now possible to infer general principles underlying enhancer evolution. However, these studies have also revealed numerous discrepancies and paradoxes, suggesting that the mechanistic causes and modes of cis-regulatory evolution are still not well understood and are probably much more complex than generally appreciated. Here, we argue that the mutational mechanisms and genomic regions generating new regulatory activities must comply with the constraints imposed by the molecular properties of cis-regulatory elements (CREs) and the organizational features of long-range chromatin interactions. Accordingly, we propose a new integrative evolutionary framework for cis-regulatory evolution based on two major premises for the origin of novel enhancer activity: (i) an accessible chromatin environment and (ii) compatibility with the 3D structure and interactions of pre-existing CREs. Mechanisms and DNA sequences not fulfilling these premises, will be less likely to have a measurable impact on gene expression and as such, will have a minor contribution to the evolution of gene regulation. Finally, we discuss current comparative cis-regulatory data under the light of this new evolutionary model, and propose that the two most prominent mechanisms for the evolution of cis-regulatory changes are the overprinting of ancestral CREs and the exaptation of transposable elements.

Reinforcing the egg-timer: recruitment of novel lophotrochozoa homeobox genes to early and late development in the pacific oyster

The metazoan superclade Lophotrochozoa includes mollusks, annelids, and several other animal phyla. It is reasonable to assume that this organismal diversity may be traced, in part, to changes in developmentally important genes, such as the homeobox genes. Although most comparative studies have focussed on ancient homeobox gene families conserved across bilaterians, there are also "novel" homeobox genes that have arisen more recently in evolution, presumably by duplication followed by radical divergence and functional change. We classify 136 homeobox genes in the genome sequence of the Pacific oyster, Crassostrea gigas. The genome shows an unusually low degree of homeobox gene clustering, with disruption of the NK, Hox, and ParaHox gene clusters. Among the oyster genes, 31 do not fall into ancient metazoan or bilaterian homeobox gene families; we deduce that they originated in the lophotrochozoan clade. We compared eight lophotrochozoan genomes to trace the pattern of homeobox gene evolution across this clade, allowing us to define 19 new lophotrochozoan-specific clades within the ANTP, PRD, TALE, ZF, SIX, and CUT classes. Using transcriptome data, we compared temporal expression of each homeobox gene in oyster development, and discovered that the lophotrochozoan-specific homeobox genes have peak expression either in early development (egg to gastrula) or in late development (after the trochophore larval stage), but rarely in between. This finding is consistent with the egg-timer, hourglass or phylotypic stage model of developmental evolution, in which there is a conserved central phase of development, but more evolutionarily labile early and late phases.

Genomics of ecological adaptation in cactophilic Drosophila

Cactophilic Drosophila species provide a valuable model to study gene–environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii–D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii–D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

The evolution of fungal metabolic pathways

Fungi contain a remarkable range of metabolic pathways, sometimes encoded by gene clusters, enabling them to digest most organic matter and synthesize an array of potent small molecules. Although metabolism is fundamental to the fungal lifestyle, we still know little about how major evolutionary processes, such as gene duplication (GD) and horizontal gene transfer (HGT), have interacted with clustered and non-clustered fungal metabolic pathways to give rise to this metabolic versatility. We examined the synteny and evolutionary history of 247,202 fungal genes encoding enzymes that catalyze 875 distinct metabolic reactions from 130 pathways in 208 diverse genomes. We found that gene clustering varied greatly with respect to metabolic category and lineage; for example, clustered genes in Saccharomycotina yeasts were overrepresented in nucleotide metabolism, whereas clustered genes in Pezizomycotina were more common in lipid and amino acid metabolism. The effects of both GD and HGT were more pronounced in clustered genes than in their non-clustered counterparts and were differentially distributed across fungal lineages; specifically, GD, which was an order of magnitude more abundant than HGT, was most frequently observed in Agaricomycetes, whereas HGT was much more prevalent in Pezizomycotina. The effect of HGT in some Pezizomycotina was particularly strong; for example, we identified 111 HGT events associated with the 15 Aspergillus genomes, which sharply contrasts with the 60 HGT events detected for the 48 genomes from the entire Saccharomycotina subphylum. Finally, the impact of GD within a metabolic category was typically consistent across all fungal lineages, whereas the impact of HGT was variable. These results indicate that GD is the dominant process underlying fungal metabolic diversity, whereas HGT is episodic and acts in a category- or lineage-specific manner. Both processes have a greater impact on clustered genes, suggesting that metabolic gene clusters represent hotspots for the generation of fungal metabolic diversity.

Deep conservation of cis-regulatory elements in metazoans

Despite the vast morphological variation observed across phyla, animals share multiple basic developmental processes orchestrated by a common ancestral gene toolkit. These genes interact with each other building complex gene regulatory networks (GRNs), which are encoded in the genome by cis-regulatory elements (CREs) that serve as computational units of the network. Although GRN subcircuits involved in ancient developmental processes are expected to be at least partially conserved, identification of CREs that are conserved across phyla has remained elusive. Here, we review recent studies that revealed such deeply conserved CREs do exist, discuss the difficulties associated with their identification and describe new approaches that will facilitate this search.

AmphiEST: Enabling comparative analysis of ESTs from five developmental stages of amphioxus

Amphioxus has been an important model for understanding the evolution of chordates and origin of vertebrates. Comparative transcriptome analysis can facilitate delineation of gene expression patterns of amphioxus at different developmental stages. So far, however, few such analyses have been performed. Here we have systematically compared amphioxus ESTs from five developmental stages. For the egg, gastrula, neurula, larva and adult stages, amphioxus ESTs were assembled, respectively, into 3364, 3230, 10,299, 4052 and 3866 contigs, and 193, 183, 933, 178 and 151 singlets. 25,796 ORFs were identified, of which 6,529 predicted ORFs (egg, 922; gastrula, 997; neurula, 2,159; larva, 1,014; adult, 1,437) have found matches from the UniProt database. For all five stages, more GO terms were assigned to stage-specific contigs than to repertoire ones. Respectively, 2, 5, 11, 4 and 2 miRNAs were matched with ESTs from the five stages, and 12,232, 10,147, 43,170, 8,049 and 12,811 miRNA target sites were predicted from the 3' UTR sequences. A web interface was developed for visualization of the contigs, miRNAs and annotations (http://www.rich.yunda.org/test/amphioxusest/). The constructed AmphiEST database may serve as an essential resource for future studies of the development and evolution of amphioxus.

Extensive conservation of ancient microsynteny across metazoans due to cis-regulatory constraints

The order of genes in eukaryotic genomes has generally been assumed to be neutral, since gene order is largely scrambled over evolutionary time. Only a handful of exceptional examples are known, typically involving deeply conserved clusters of tandemly duplicated genes (e.g., Hox genes and histones). Here we report the first systematic survey of microsynteny conservation across metazoans, utilizing 17 genome sequences. We identified nearly 600 pairs of unrelated genes that have remained tightly physically linked in diverse lineages across over 600 million years of evolution. Integrating sequence conservation, gene expression data, gene function, epigenetic marks, and other genomic features, we provide extensive evidence that many conserved ancient linkages involve (1) the coordinated transcription of neighboring genes, or (2) genomic regulatory blocks (GRBs) in which transcriptional enhancers controlling developmental genes are contained within nearby bystander genes. In addition, we generated ChIP-seq data for key histone modifications in zebrafish embryos, which provided further evidence of putative GRBs in embryonic development. Finally, using chromosome conformation capture (3C) assays and stable transgenic experiments, we demonstrate that enhancers within bystander genes drive the expression of genes such as Otx and Islet, critical regulators of central nervous system development across bilaterians. These results suggest that ancient genomic functional associations are far more common than previously thought—involving ∼12% of the ancestral bilaterian genome—and that cis-regulatory constraints are crucial in determining metazoan genome architecture.

Widespread recurrent evolution of genomic features

The recent explosion of genome sequences from all major phylogenetic groups has unveiled an unexpected wealth of cases of recurrent evolution of strikingly similar genomic features in different lineages. Here, we review the diverse known types of recurrent evolution in eukaryotic genomes, with a special focus on metazoans, ranging from reductive genome evolution to origins of splice-leader trans-splicing, from tandem exon duplications to gene family expansions. We first propose a general classification scheme for evolutionary recurrence at the genomic level, based on the type of driving force-mutation or selection-and the environmental and genomic circumstances underlying these forces. We then discuss various cases of recurrent genomic evolution under this scheme. Finally, we provide a broader context for repeated genomic evolution, including the unique relationship of genomic recurrence with the genotype-phenotype map, and the ways in which the study of recurrent genomic evolution can be used to understand fundamental evolutionary processes.

An ancient genomic regulatory block conserved across bilaterians and its dismantling in tetrapods by retrogene replacement

Developmental genes are regulated by complex, distantly located cis-regulatory modules (CRMs), often forming genomic regulatory blocks (GRBs) that are conserved among vertebrates and among insects. We have investigated GRBs associated with Iroquois homeobox genes in 39 metazoans. Despite 600 million years of independent evolution, Iroquois genes are linked to ankyrin-repeat-containing Sowah genes in nearly all studied bilaterians. We show that Iroquois-specific CRMs populate the Sowah locus, suggesting that regulatory constraints underlie the maintenance of the Iroquois-Sowah syntenic block. Surprisingly, tetrapod Sowah orthologs are intronless and not associated with Iroquois; however, teleost and elephant shark data demonstrate that this is a derived feature, and that many Iroquois-CRMs were ancestrally located within Sowah introns. Retroposition, gene, and genome duplication have allowed selective elimination of Sowah exons from the Iroquois regulatory landscape while keeping associated CRMs, resulting in large associated gene deserts. These results highlight the importance of CRMs in imposing constraints to genome architecture, even across large phylogenetic distances, and of gene duplication-mediated genetic redundancy to disentangle these constraints, increasing genomic plasticity.

An evolutionarily conserved three-dimensional structure in the vertebrate Irx clusters facilitates enhancer sharing and coregulation

Developmental gene clusters are paradigms for the study of gene regulation; however, the mechanisms that mediate phenomena such as coregulation and enhancer sharing remain largely elusive. Here we address this issue by analysing the vertebrate Irx clusters. We first present a deep enhancer screen of a 2-Mbp span covering the IrxA cluster. Using chromosome conformation capture, we show that enhancer sharing is widespread within the cluster, explaining its evolutionarily conserved organization. We also identify a three-dimensional architecture, probably formed through interactions with CCCTC-binding factor, which is present within both Irx clusters of mouse, Xenopus and zebrafish. This architecture brings the promoters of the first two genes together in the same chromatin landscape. We propose that this unique and evolutionarily conserved genomic architecture of the vertebrate Irx clusters is essential for the coregulation of the first two genes and simultaneously maintains the third gene in a partially independent regulatory landscape.

Annelid Distal-less/Dlx duplications reveal varied post-duplication fates

Background

Dlx (Distal-less) genes have various developmental roles and are widespread throughout the animal kingdom, usually occurring as single copy genes in non-chordates and as multiple copies in most chordate genomes. While the genomic arrangement and function of these genes is well known in vertebrates and arthropods, information about Dlx genes in other organisms is scarce. We investigate the presence of Dlx genes in several annelid species and examine Dlx gene expression in the polychaete Pomatoceros lamarckii.

Results

Two Dlx genes are present in P. lamarckii, Capitella teleta and Helobdella robusta. The C. teleta Dlx genes are closely linked in an inverted tail-to-tail orientation, reminiscent of the arrangement of vertebrate Dlx pairs, and gene conversion appears to have had a role in their evolution. The H. robusta Dlx genes, however, are not on the same genomic scaffold and display divergent sequences, while, if the P. lamarckii genes are linked in a tail-to-tail orientation they are a minimum of 41 kilobases apart and show no sign of gene conversion. No expression in P. lamarckii appendage development has been observed, which conflicts with the supposed conserved role of these genes in animal appendage development. These Dlx duplications do not appear to be annelid-wide, as the polychaete Platynereis dumerilii likely possesses only one Dlx gene.

Conclusions

On the basis of the currently accepted annelid phylogeny, we hypothesise that one Dlx duplication occurred in the annelid lineage after the divergence of P. dumerilii from the other lineages and these duplicates then had varied evolutionary fates in different species. We also propose that the ancestral role of Dlx genes is not related to appendage development.

Contrasting 5' and 3' evolutionary histories and frequent evolutionary convergence in Meis/hth gene structures

Organisms show striking differences in genome structure; however, the functional implications and fundamental forces that govern these differences remain obscure. The intron–exon organization of nuclear genes is involved in a particularly large variety of structures and functional roles. We performed a 22-species study of Meis/hth genes, intron-rich homeodomain-containing transcription factors involved in a wide range of developmental processes. Our study revealed three surprising results that suggest important and very different functions for Meis intron–exon structures. First, we find unexpected conservation across species of intron positions and lengths along most of the Meis locus. This contrasts with the high degree of structural divergence found in genome-wide studies and may attest to conserved regulatory elements residing within these conserved introns. Second, we find very different evolutionary histories for the 5′ and 3′ regions of the gene. The 5′-most 10 exons, which encode the highly conserved Meis domain and homeodomain, show striking conservation. By contrast, the 3′ of the gene, which encodes several domains implicated in transcriptional activation and response to cell signaling, shows a remarkably active evolutionary history, with diverse isoforms and frequent creation and loss of new exons and splice sites. This region-specific diversity suggests evolutionary “tinkering,” with alternative splicing allowing for more subtle regulation of protein function. Third, we find a large number of cases of convergent evolution in the 3′ region, including 1) parallel losses of ancestral coding sequence, 2) parallel gains of external and internal splice sites, and 3) recurrent truncation of C-terminal coding regions. These results attest to the importance of locus-specific splicing functions in differences in structural evolution across genes, as well as to commonalities of forces shaping the evolution of individual genes along different lineages.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background

The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians.

Results

In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern.

Conclusions

The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

Multiple GAL pathway gene clusters evolved independently and by different mechanisms in fungi

A notable characteristic of fungal genomes is that genes involved in successive steps of a metabolic pathway are often physically linked or clustered. To investigate how such clusters of functionally related genes are assembled and maintained, we examined the evolution of gene sequences and order in the galactose utilization (GAL) pathway in whole-genome data from 80 diverse fungi. We found that GAL gene clusters originated independently and by different mechanisms in three unrelated yeast lineages. Specifically, the GAL cluster found in Saccharomyces and Candida yeasts originated through the relocation of native unclustered genes, whereas the GAL cluster of Schizosaccharomyces yeasts was acquired through horizontal gene transfer from a Candida yeast. In contrast, the GAL cluster of Cryptococcus yeasts was assembled independently from the Saccharomyces/Candida and Schizosaccharomyces GAL clusters and coexists in the Cryptococcus genome with unclustered GAL paralogs. These independently evolved GAL clusters represent a striking example of analogy at the genomic level. We also found that species with GAL clusters exhibited significantly higher rates of GAL pathway loss than species with unclustered GAL genes. These results suggest that clustering of metabolic genes might facilitate fungal adaptation to changing environments both through the acquisition and loss of metabolic capacities.

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome-deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)-(Nk4/tinman)-(Nk3/bagpipe)-(Lbx/ladybird)-(Tlx/c15)-(Nk7)-(Nk6/hgtx)-(Nk1/slouch)-(Nk5/Hmx). Of these genes, only NKX2.6-NKX3.1, LBX1-TLX1 and LBX2-TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have "missed" certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1-Tlx1, Lbx2-Tlx2 and Nkx3.1-Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2-Hmx3 genes are also tightly linked. Finally, we show that Lbx1-Tlx1 and Hmx2-Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2-Tlx2 and Nkx3.1-Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.

Internal and external paralogy in the evolution of tropomyosin genes in metazoans

Nature contains a tremendous diversity of forms both at the organismal and genomic levels. This diversity motivates the twin central questions of molecular evolution: what are the molecular mechanisms of adaptation, and what are the functional consequences of genomic diversity. We report a 22-species comparative analysis of tropomyosin (PPM) genes, which exist in a variety of forms and are implicated in the emergence of a wealth of cellular functions, including the novel muscle functions integral to the functional diversification of bilateral animals. TPM genes encode either or both of long-form [284 amino acid (aa)] and short-form (approximately 248 aa) proteins. Consistent with a role of TPM diversification in the origins and radiation of bilaterians, we find evidence that the muscle-specific long-form protein arose in proximal bilaterian ancestors (the bilaterian 'stem'). Duplication of the 5' end of the gene led to alternative promoters encoding long- and short-form transcripts with distinct functions. This dual-function gene then underwent strikingly parallel evolution in different bilaterian lineages. In each case, recurrent tandem exon duplication and mutually exclusive alternative splicing of the duplicates, with further association between these alternatively spliced exons along the gene, led to long- and short-form-specific exons, allowing for gradual emergence of alternative "internal paralogs" within the same gene. We term these Mutually exclusively Alternatively spliced Tandemly duplicated Exon sets "MATEs". This emergence of internal paralogs in various bilaterians has employed every single TPM exon in at least one lineage and reaches striking levels of divergence with up to 77% of long- and short-form transcripts being transcribed from different genomic regions. Interestingly, in some lineages, these internal alternatively spliced paralogs have subsequently been "externalized" by full gene duplication and reciprocal retention/loss of the two transcript isoforms, a particularly clear case of evolution by subfunctionalization. This parallel evolution of TPM genes in diverse metazoans attests to common selective forces driving divergence of different gene transcripts and represents a striking case of emergence of evolutionary novelty by alternative splicing.

It's a long way from amphioxus: descendants of the earliest chordate

The origin of chordates and the consequent genesis of vertebrates were major events in natural history. The amphioxus (lancelet) is now recognised as the closest extant relative to the stem chordate and is the only living invertebrate that retains a vertebrate-like development and body plan through its lifespan, despite more than 500 million years of independent evolution from the stem vertebrate. The inspiring data coming from its recently sequenced genome confirms that amphioxus has a prototypical chordate genome with respect to gene content and structure, and even chromosomal organisation. Pushed by joint efforts of amphioxus researchers, amphioxus is now entering a new era, namely its maturation as a laboratory model, through the availability of a large amount of molecular data and the advent of experimental manipulation of the embryo. These two facts may well serve to illuminate the hidden secrets of the genetic changes that generated, among other vertebrates, ourselves.

Evolutionary history of the iroquois/Irx genes in metazoans

BACKGROUND

The iroquois (iro/Irx) genes encode transcriptional regulators that belong to the TALE superclass of homeodomain proteins and have key functions during development in both vertebrates and insects. The Irx genes occur in one or two genomic clusters containing three genes each within the Drosophila and several vertebrate genomes, respectively. The similar genomic organization in Drosophila and vertebrates is widely considered as a result of convergent evolution, due to independent tandem gene duplications. In this study, we investigate the evolutionary history of the Irx genes at the scale of the whole metazoan kingdom.

RESULTS

We identified in silico the putative full complement of Irx genes in the sequenced genomes of 36 different species representative of the main metazoan lineages, including non bilaterian species, several arthropods, non vertebrate chordates, and a basal vertebrate, the sea lamprey. We performed extensive phylogenetic analyses of the identified Irx genes and defined their genomic organizations. We found that, in most species, there are several Irx genes, these genes form two to four gene clusters, and the Irx genes are physically linked to a structurally and functionally unrelated gene known as CG10632 in Drosophila.

CONCLUSION

Three main conclusions can be drawn from our study. First, an Irx cluster composed of two genes, araucan/caupolican and mirror, is ancestral to the crustaceans+insects clade and has been strongly conserved in this clade. Second, three Irx genes were probably present in the last common ancestor of vertebrates and the duplication that has given rise to the six genes organized into two clusters found in most vertebrates, likely occurred in the gnathostome lineage after its separation from sea lampreys. Third, the clustered organization of the Irx genes in various evolutionary lineages may represent an exceptional case of convergent evolution or may point to the existence of an Irx gene cluster ancestral to bilaterians.