hsd17b1 is a key gene for ovarian differentiation of the Siberian sturgeon.

This is the first work using gonads from undifferentiated, genetically-sexed Siberian sturgeon describing expression changes in genes related to steroid synthesis and female and male sex differentiation. One factor identified as relevant for ovarian differentiation was the gene coding for the enzyme Hsd17b1, which converts estrone into estradiol-17ß. hsd17b1 was highly activated in female gonads at 2.5 months of age, around the onset of sex differentiation, preceding activation of two other genes involved in estrogen production (cyp19a1 and foxl2). hsd17b1 was also strongly repressed in males. Two known foxl2 paralogs are found in Siberian sturgeon-foxl2 and foxl2l-but only foxl2 appeared to be associated with ovarian differentiation. With regard to the male pathway, neither 11-oxygenated androgens nor classic male genes (amh, dmrt1, sox9, and dhh) were found to be involved in male sex differentiation, leaving open the question of which genes participate in early male gonad development in this ancient fish. Taken together, these results indicate an estrogen-dependence of female sex differentiation and 11-oxygenated androgen-independence of male sex differentiation.

Evidence for DNA Sequence Encoding of an Accessible Nucleosomal Array across Vertebrates.

Nucleosome-depleted regions around which nucleosomes order following the "statistical" positioning scenario were recently shown to be encoded in the DNA sequence in human. This intrinsic nucleosomal ordering strongly correlates with oscillations in the local GC content as well as with the interspecies and intraspecies mutation profiles, revealing the existence of both positive and negative selection. In this letter, we show that these predicted nucleosome inhibitory energy barriers (NIEBs) with compacted neighboring nucleosomes are indeed ubiquitous to all vertebrates tested. These 1 kb-sized chromatin patterns are widely distributed along vertebrate chromosomes, overall covering more than a third of the genome. We have previously observed in human deviations from neutral evolution at these genome-wide distributed regions, which we interpreted as a possible indication of the selection of an open, accessible, and dynamic nucleosomal array to constitutively facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner. As a first, very appealing observation supporting this hypothesis, we report evidence of a strong association between NIEB borders and the poly(A) tails of Alu sequences in human. These results suggest that NIEBs provide adequate chromatin patterns favorable to the integration of Alu retrotransposons and, more generally to various transposable elements in the genomes of primates and other vertebrates.

The evolutionary conservation of the A Disintegrin-like and Metalloproteinase domain with Thrombospondin-1 motif metzincins across vertebrate species and their expression in teleost zebrafish.

BACKGROUND: The A Disintegrin-like and Metalloproteinase domain with Thrombospondin-1 motifs (ADAMTS) enzymes comprise 19 mammalian zinc-dependent metalloproteinases (metzincins) with homologues in a wide range of invertebrates. ADAMTS enzymes have a broad range of functions in development and diseases due to their extracellular matrix remodelling activity. Here, we report a detailed characterisation of their evolutionary conservation across vertebrates. RESULTS: Using bioinformatics complemented with de novo sequencing, gene sequences for ADAMTS enzymes were obtained from a variety of organisms. Detailed evolutionary analyses revealed a high level of conservation across vertebrates with evidence of ADAMTS gene expansion during two rounds of whole genome duplication (WGD) in vertebrates, while tandem duplication events and gene loss were also apparent. However, the additional round of teleost-specific WGD did not have a significant effect on ADAMTS gene family members suggesting their conserved roles have remained constant in teleost fish. Quantitative reverse-transcriptase polymerase chain reaction analysis revealed dynamic expression of adamts genes throughout zebrafish embryonic development reflecting the key conserved roles they play in vertebrate embryogenesis. Notably, several adamts mRNAs were maternally expressed with a dramatic increase in mRNA levels coinciding with zygotic expression and organogenesis. Broad adamts mRNA expression was also demonstrated in several adult organs indicating potential roles in adult homeostasis. CONCLUSIONS: Our data highlight the evolution of the ADAMTS gene family through duplication processes across metazoans supplemented by a burst of amplification through vertebrate WGD events. It also strongly posits the zebrafish as a potential model species to further elucidate the function of ADAMTS enzymes during vertebrate development.

Estrogen-related receptor γ is an in vivo receptor of bisphenol A.

Bisphenol A (BPA) is an endocrine disruptor that displays estrogenic activity. Several reports suggest that BPA may have estrogen receptor-independent effects. In zebrafish, 50 µM BPA exposure induces otic vesicle abnormalities, including otolith aggregation. The purpose of this study was to test if BPA action was mediated in vivo during zebrafish development by the orphan nuclear estrogen related receptor (ERR) γ. Combining pharmacological and functional approaches, we demonstrate that the zebrafish ERRγ mediates BPA-induced malformations in otoliths. Using different bisphenol derivatives, we show that different compounds can induce a similar otolith phenotype than BPA and that the binding affinity of these derivatives to the zebrafish ERRγ correlates with their ability to induce otolith malformations. Morpholino knockdown of ERRγ function suppresses the BPA effect on otoliths whereas overexpression of ERRγ led to a BPA-like otolith phenotype. Moreover, a subphenotypical dose of BPA (1 µM) combined with ERRγ overexpression led to a full-dose (50 µM) BPA otolith phenotype. We therefore conclude that ERRγ mediates the otic vesicle phenotype generated by BPA. Our results suggest that the range of pathways perturbed by this compound and its potential harmful effect are larger than expected.-Tohmé, M., Prud'homme, S. M., Boulahtouf, A., Samarut, E., Brunet, F., Bernard, L., Bourguet, W., Gibert, Y., Balaguer, P., Laudet, V. Estrogen-related receptor γ is an in vivo receptor of bisphenol A.

The N-terminal domains of TRF1 and TRF2 regulate their ability to condense telomeric DNA.

TRF1 and TRF2 are key proteins in human telomeres, which, despite their similarities, have different behaviors upon DNA binding. Previous work has shown that unlike TRF1, TRF2 condenses telomeric, thus creating consequential negative torsion on the adjacent DNA, a property that is thought to lead to the stimulation of single-strand invasion and was proposed to favor telomeric DNA looping. In this report, we show that these activities, originating from the central TRFH domain of TRF2, are also displayed by the TRFH domain of TRF1 but are repressed in the full-length protein by the presence of an acidic domain at the N-terminus. Strikingly, a similar repression is observed on TRF2 through the binding of a TERRA-like RNA molecule to the N-terminus of TRF2. Phylogenetic and biochemical studies suggest that the N-terminal domains of TRF proteins originate from a gradual extension of the coding sequences of a duplicated ancestral gene with a consequential progressive alteration of the biochemical properties of these proteins. Overall, these data suggest that the N-termini of TRF1 and TRF2 have evolved to finely regulate their ability to condense DNA.

Structural shifts of aldehyde dehydrogenase enzymes were instrumental for the early evolution of retinoid-dependent axial patterning in metazoans.

Aldehyde dehydrogenases (ALDHs) catabolize toxic aldehydes and process the vitamin A-derived retinaldehyde into retinoic acid (RA), a small diffusible molecule and a pivotal chordate morphogen. In this study, we combine phylogenetic, structural, genomic, and developmental gene expression analyses to examine the evolutionary origins of ALDH substrate preference. Structural modeling reveals that processing of small aldehydes, such as acetaldehyde, by ALDH2, versus large aldehydes, including retinaldehyde, by ALDH1A is associated with small versus large substrate entry channels (SECs), respectively. Moreover, we show that metazoan ALDH1s and ALDH2s are members of a single ALDH1/2 clade and that during evolution, eukaryote ALDH1/2s often switched between large and small SECs after gene duplication, transforming constricted channels into wide opened ones and vice versa. Ancestral sequence reconstructions suggest that during the evolutionary emergence of RA signaling, the ancestral, narrow-channeled metazoan ALDH1/2 gave rise to large ALDH1 channels capable of accommodating bulky aldehydes, such as retinaldehyde, supporting the view that retinoid-dependent signaling arose from ancestral cellular detoxification mechanisms. Our analyses also indicate that, on a more restricted evolutionary scale, ALDH1 duplicates from invertebrate chordates (amphioxus and ascidian tunicates) underwent switches to smaller and narrower SECs. When combined with alterations in gene expression, these switches led to neofunctionalization from ALDH1-like roles in embryonic patterning to systemic, ALDH2-like roles, suggesting functional shifts from signaling to detoxification.

Molecular adaptation and resilience of the insect's nuclear receptor USP.

BACKGROUND: The maintenance of biological systems requires plasticity and robustness. The function of the ecdysone receptor, a heterodimer composed of the nuclear receptors ECR (NR1H1) and USP (NR2B4), was maintained in insects despite a dramatic divergence that occurred during the emergence of Mecopterida. This receptor is therefore a good model to study the evolution of plasticity. We tested the hypothesis that selection has shaped the Ligand-Binding Domain (LBD) of USP during evolution of Mecopterida. RESULTS: We isolated usp and cox1 in several species of Drosophilidae, Tenebrionidae and Blattaria and estimated non-synonymous/synonymous rate ratios using maximum-likelihood methods and codon-based substitution models. Although the usp sequences were mainly under negative selection, we detected relaxation at residues located on the surface of the LBD within Mecopterida families. Using branch-site models, we also detected changes in selective constraints along three successive branches of the Mecopterida evolution. Residues located at the bottom of the ligand-binding pocket (LBP) underwent strong positive selection during the emergence of Mecopterida. This change is correlated with the acquisition of a large LBP filled by phospholipids that probably allowed the stabilisation of the new Mecopterida structure. Later, when the two subgroups of Mecopterida (Amphiesmenoptera: Lepidoptera, Trichoptera; Antliophora: Diptera, Mecoptera, Siphonaptera) diverged, the same positions became under purifying selection. Similarly, several positions of the heterodimerisation interface experienced positive selection during the emergence of Mecopterida, rapidly followed by a phase of constrained evolution. An enlargement of the heterodimerisation surface is specific for Mecopterida and was associated with a reinforcement of the obligatory partnership between ECR and USP, at the expense of homodimerisation. CONCLUSIONS: In order to explain the episodic mode of evolution of USP, we propose a model in which the molecular adaptation of this protein is seen as a process of resilience for the maintenance of the ecdysone receptor functionality.

Expression and phylogeny of candidate genes for sex differentiation in a primitive fish species, the Siberian sturgeon, Acipenser baerii.

The molecular mechanisms underlying testis differentiation in basal actinopterygian fish remains poorly understood. The sex differentiation period was investigated in the Siberian sturgeon, Acipenser baerii, by expression profiling of Sertoli cell transcription factors (dmrt1, sox9) that control testis differentiation in vertebrates; Leydig cell factors (cyp17a1, star) affecting androgen production; the androgen receptor (ar); a growth factor controlling testis development (igf1); and a gene coding for a gonadotropin hormone (lh). Two genes were characterised for the first time in the Siberian sturgeon (dmrt1, cyp17a1), while the others came from public databases. Sturgeon gonad development is very slow, with a late sexual differentiation time during their juvenile stage, and are still immature at 3 years of age. Immature fish showed a sex-dimorphic pattern; all the genes studied displayed a higher expression level in male gonads. We took advantage of the presence of juvenile fish with pre- and post-differentiated gonads (16 and 18 months old) to characterise them at the molecular level. The post-differentiated fish displayed a sex dimorphism of gene expression in their gonads for all genes studied, with the exception of sox9. The trends in undifferentiated fish lead us to propose that sturgeons undergoing male differentiation express high levels of Sertoli cell factors (dmrt1, sox9) and of genes involved in the production and receptivity of androgens (cyp17a1, star and ar) together with lh. Expression profiles and phylogenetic studies suggest that these genes are potential regulators of testis development in the Siberian sturgeon.

Amphioxus postembryonic development reveals the homology of chordate metamorphosis.

Most studies in evolution are centered on how homologous genes, structures, and/or processes appeared and diverged. Although historical homology is well defined as a concept, in practice its establishment can be problematic, especially for some morphological traits or developmental processes. Metamorphosis in chordates is such an enigmatic character. Defined as a spectacular postembryonic larva-to-adult transition, it shows a wide morphological diversity between the different chordate lineages, suggesting that it might have appeared several times independently. In vertebrates, metamorphosis is triggered by binding of the thyroid hormones (THs) T(4) and T(3) to thyroid-hormone receptors (TRs). Here we show that a TH derivative, triiodothyroacetic acid (TRIAC), induces metamorphosis in the cephalochordate amphioxus. The amphioxus TR (amphiTR) mediates spontaneous and TRIAC-induced metamorphosis because it strongly binds to TRIAC, and a specific TR antagonist, NH3, inhibits both spontaneous and TRIAC-induced metamorphosis. Moreover, as in amphibians, amphiTR expression levels increase around metamorphosis and are enhanced by THs. Therefore, TH-regulated metamorphosis, mediated by TR, is an ancestral feature of all chordates. This conservation of a regulatory network supports the homology of metamorphosis in the chordate lineage.

Repetitive element-mediated recombination as a mechanism for new gene origination in Drosophila.

Previous studies of repetitive elements (REs) have implicated a mechanistic role in generating new chimerical genes. Such examples are consistent with the classic model for exon shuffling, which relies on non-homologous recombination. However, recent data for chromosomal aberrations in model organisms suggest that ectopic homology-dependent recombination may also be important. Lack of a dataset comprising experimentally verified young duplicates has hampered an effective examination of these models as well as an investigation of sequence features that mediate the rearrangements. Here we use approximately 7,000 cDNA probes (approximately 112,000 primary images) to screen eight species within the Drosophila melanogaster subgroup and identify 17 duplicates that were generated through ectopic recombination within the last 12 mys. Most of these are functional and have evolved divergent expression patterns and novel chimeric structures. Examination of their flanking sequences revealed an excess of repetitive sequences, with the majority belonging to the transposable element DNAREP1 family, associated with the new genes. Our dataset strongly suggests an important role for REs in the generation of chimeric genes within these species.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution.

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype.

Tetraodon nigroviridis is a freshwater puffer fish with the smallest known vertebrate genome. Here, we report a draft genome sequence with long-range linkage and substantial anchoring to the 21 Tetraodon chromosomes. Genome analysis provides a greatly improved fish gene catalogue, including identifying key genes previously thought to be absent in fish. Comparison with other vertebrates and a urochordate indicates that fish proteins have diverged markedly faster than their mammalian homologues. Comparison with the human genome suggests approximately 900 previously unannotated human genes. Analysis of the Tetraodon and human genomes shows that whole-genome duplication occurred in the teleost fish lineage, subsequent to its divergence from mammals. The analysis also makes it possible to infer the basic structure of the ancestral bony vertebrate genome, which was composed of 12 chromosomes, and to reconstruct much of the evolutionary history of ancient and recent chromosome rearrangements leading to the modern human karyotype.

Sequencing the Genome of Indian Flying Fox, Natural Reservoir of Nipah Virus, Using Hybrid Assembly and Conservative Secondary Scaffolding.

Indian fruit bats, flying fox Pteropus medius was identified as an asymptomatic natural host of recently emerged Nipah virus, which is known to induce a severe infectious disease in humans. The absence of P. medius genome sequence presents an important obstacle for further studies of virus-host interactions and better understanding of mechanisms of zoonotic viral emergence. Generation of the high-quality genome sequence is often linked to a considerable effort associated to elevated costs. Although secondary scaffolding methods have reduced sequencing expenses, they imply the development of new tools for the integration of different data sources to achieve more reliable sequencing results. We initially sequenced the P. medius genome using the combination of Illumina paired-end and Nanopore sequencing, with a depth of 57.4x and 6.1x, respectively. Then, we introduced the novel scaff2link software to integrate multiple sources of information for secondary scaffolding, allowing to remove the association with discordant information among two sources. Different quality metrics were next produced to validate the benefits from secondary scaffolding. The P. medius genome, assembled by this method, has a length of 1,985 Mb and consists of 33,613 contigs and 16,113 scaffolds with an NG50 of 19 Mb. At least 22.5% of the assembled sequences is covered by interspersed repeats already described in other species and 19,823 coding genes are annotated. Phylogenetic analysis demonstrated the clustering of P. medius genome with two other Pteropus bat species, P. alecto and P. vampyrus, for which genome sequences are currently available. SARS-CoV entry receptor ACE2 sequence of P. medius was 82.7% identical with ACE2 of Rhinolophus sinicus bats, thought to be the natural host of SARS-CoV. Altogether, our results confirm that a lower depth of sequencing is enough to obtain a valuable genome sequence, using secondary scaffolding approaches and demonstrate the benefits of the scaff2link application. The genome sequence is now available to the scientific community to (i) proceed with further genomic analysis of P. medius, (ii) to characterize the underlying mechanism allowing Nipah virus maintenance and perpetuation in its bat host, and (iii) to monitor their evolutionary pathways toward a better understanding of bats' ability to control viral infections.

Origin and evolution of the Notch signalling pathway: an overview from eukaryotic genomes.

BACKGROUND: Of the 20 or so signal transduction pathways that orchestrate cell-cell interactions in metazoans, seven are involved during development. One of these is the Notch signalling pathway which regulates cellular identity, proliferation, differentiation and apoptosis via the developmental processes of lateral inhibition and boundary induction. In light of this essential role played in metazoan development, we surveyed a wide range of eukaryotic genomes to determine the origin and evolution of the components and auxiliary factors that compose and modulate this pathway. RESULTS: We searched for 22 components of the Notch pathway in 35 different species that represent 8 major clades of eukaryotes, performed phylogenetic analyses and compared the domain compositions of the two fundamental molecules: the receptor Notch and its ligands Delta/Jagged. We confirm that a Notch pathway, with true receptors and ligands is specific to the Metazoa. This study also sheds light on the deep ancestry of a number of genes involved in this pathway, while other members are revealed to have a more recent origin. The origin of several components can be accounted for by the shuffling of pre-existing protein domains, or via lateral gene transfer. In addition, certain domains have appeared de novo more recently, and can be considered metazoan synapomorphies. CONCLUSION: The Notch signalling pathway emerged in Metazoa via a diversity of molecular mechanisms, incorporating both novel and ancient protein domains during eukaryote evolution. Thus, a functional Notch signalling pathway was probably present in Urmetazoa.

Isotopic characteristics of the Garonne River and its tributaries.

The Garonne is the largest river in the south-west of France, and its drainage basin stretches between the Pyrénées and the Massif Central mountains. Until now, no water stable isotope study has been performed on the whole Garonne river basin which is composed of different geological substrata, and where the water resources are limited during the dry summer period. This study focuses on the Garonne river and its tributaries from the Pyrénées foothill upstream to its confluence with the Lot River downstream. The aim of the study is to determine the origins of the surface waters using their chemical and stable isotopic compositions ((18)O, D and (13)C), to better understand their circulation within the drainage basin and to assess the anthropogenic influences. The Garonne displays a specific (18)O seasonal effect, and keeps its Pyrénean characteristics until its confluence with the Tarn River. The difference in the dissolved inorganic carbon (DIC) comes mainly from the change in lithology between the Pyrénées and the Massif Central mountains. Agriculture activity is only detected in the small tributaries.

OTX5 regulates pineal expression of the zebrafish REV-ERB alpha through a new DNA binding site.

The pineal gland plays a central role in the photoneuroendocrine system and acts as a photosensory organ in lower vertebrates. The orphan nuclear receptor Rev-erbalpha (NR1D1) has previously been shown to be expressed in the pineal and to be regulated with a robust circadian rhythm during zebrafish embryogenesis. This early pineal expression is under the control of the transcription factor Orthodenticle homeobox 5 (Otx5). In this paper, we show that Otx5 regulates the second zfRev-erbalpha promoter, ZfP2. Despite the absence of a classical Otx-binding site within ZfP2, this regulation depends on the integrity of the Otx5 homeodomain. Mapping experiments as well as EMSAs show that this interaction between Otx5 and ZfP2 depends on a noncanonical bipartite Otx-binding site (GANNCTTA and TAAA) that we called pineal expression related element (PERE). We showed that PERE is necessary for pineal expression in vivo by injecting zebrafish embryos with wild type and mutated versions of zfRev-erbalpha promoter fused to green fluorescent protein. Interestingly, PERE is found upstream of other genes expressed in the pineal gland, suggesting that it may play an important role in governing pineal expression. Our data establish that PERE is a novel cis-acting element contributing to pineal-specific gene expression and to Otx target gene regulation.

Highly variable rates of genome rearrangements between hemiascomycetous yeast lineages.

Hemiascomycete yeasts cover an evolutionary span comparable to that of the entire phylum of chordates. Since this group currently contains the largest number of complete genome sequences it presents unique opportunities to understand the evolution of genome organization in eukaryotes. We inferred rates of genome instability on all branches of a phylogenetic tree for 11 species and calculated species-specific rates of genome rearrangements. We characterized all inversion events that occurred within synteny blocks between six representatives of the different lineages. We show that the rates of macro- and microrearrangements of gene order are correlated within individual lineages but are highly variable across different lineages. The most unstable genomes correspond to the pathogenic yeasts Candida albicans and Candida glabrata. Chromosomal maps have been intensively shuffled by numerous interchromosomal rearrangements, even between species that have retained a very high physical fraction of their genomes within small synteny blocks. Despite this intensive reshuffling of gene positions, essential genes, which cluster in low recombination regions in the genome of Saccharomyces cerevisiae, tend to remain syntenic during evolution. This work reveals that the high plasticity of eukaryotic genomes results from rearrangement rates that vary between lineages but also at different evolutionary times of a given lineage.

The fate of the duplicated androgen receptor in fishes: a late neofunctionalization event?

BACKGROUND: Based on the observation of an increased number of paralogous genes in teleost fishes compared with other vertebrates and on the conserved synteny between duplicated copies, it has been shown that a whole genome duplication (WGD) occurred during the evolution of Actinopterygian fish. Comparative phylogenetic dating of this duplication event suggests that it occurred early on, specifically in teleosts. It has been proposed that this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish, notably by allowing the sub- or neo-functionalization of many duplicated genes. RESULTS: In this paper, we studied in a wide range of Actinopterygians the duplication and fate of the androgen receptor (AR, NR3C4), a nuclear receptor known to play a key role in sex-determination in vertebrates. The pattern of AR gene duplication is consistent with an early WGD event: it has been duplicated into two genes AR-A and AR-B after the split of the Acipenseriformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. Genomic and syntenic analyses in addition to lack of PCR amplification show that one of the duplicated copies, AR-B, was lost in several basal Clupeocephala such as Cypriniformes (including the model species zebrafish), Siluriformes, Characiformes and Salmoniformes. Interestingly, we also found that, in basal teleost fish (Osteoglossiformes and Anguilliformes), the two copies remain very similar, whereas, specifically in Percomorphs, one of the copies, AR-B, has accumulated substitutions in both the ligand binding domain (LBD) and the DNA binding domain (DBD). CONCLUSION: The comparison of the mutations present in these divergent AR-B with those known in human to be implicated in complete, partial or mild androgen insensitivity syndrome suggests that the existence of two distinct AR duplicates may be correlated to specific functional differences that may be connected to the well-known plasticity of sex determination in fish. This suggests that three specific events have shaped the present diversity of ARs in Actinopterygians: (i) early WGD, (ii) parallel loss of one duplicate in several lineages and (iii) putative neofunctionalization of the same duplicate in percomorphs, which occurred a long time after the WGD.

The amphioxus genome enlightens the evolution of the thyroid hormone signaling pathway.

Thyroid hormones (THs) have pleiotropic effects on vertebrate development, with amphibian metamorphosis as the most spectacular example. However, developmental functions of THs in non-vertebrate chordates are largely hypothetical and even TH endogenous production has been poorly investigated. In order to get better insight into the evolution of the thyroid hormone signaling pathway in chordates, we have taken advantage of the recent release of the amphioxus genome. We found amphioxus homologous sequences to most of the genes encoding proteins involved in thyroid hormone signaling in vertebrates, except the fast-evolving thyroglobulin: sodium iodide symporter, thyroid peroxidase, deiodinases, thyroid hormone receptor, TBG, and CTHBP. As only some genes encoding proteins involved in TH synthesis regulation were retrieved (TRH, TSH receptor, and CRH receptor but not their corresponding receptors and ligands), there may be another mode of upstream regulation of TH synthesis in amphioxus. In accord with the notion that two whole genome duplications took place at the base of the vertebrate tree, one amphioxus gene often corresponded to several vertebrate homologs. However, some amphioxus specific duplications occurred, suggesting that several steps of the TH pathway were independently elaborated in the cephalochordate and vertebrate lineages. The present results therefore indicate that amphioxus is capable of producing THs. As several genes of the TH signaling pathway were also found in the sea urchin genome, we propose that the thyroid hormone signaling pathway is of ancestral origin in chordates, if not in deuterostomes, with specific elaborations in each lineage, including amphioxus.

Nuclear hormone receptor signaling in amphioxus.

The nuclear hormone receptors (NRs) form a superfamily of transcription factors unified by conserved protein structure and mode of function. While most members of this superfamily are activated by ligands, such as thyroid hormones, steroids, vitamin D or retinoic acid, other NRs are called orphan receptors because they have no known ligand. NR-dependent signaling is crucial for vertebrate development with the majority of receptors being expressed in the developing embryo. Due to massive gene duplications during vertebrate diversification, there are usually more NRs in vertebrates than in invertebrates. In this study, we examine the evolutionary diversification of the NR superfamily and of NR-dependent signaling in chordates (vertebrates, tunicates, and amphioxus). We take advantage of the unique features of the genome of the invertebrate amphioxus, which is characterized by a vertebrate-like gene content without having undergone massive duplications, to assess the NR signaling complement (NRs and NR coregulators) of the ancestral chordate. We find 33 NRs in amphioxus, which are more NRs than originally anticipated. This increase is mainly due to an amphioxus-specific duplication of genes encoding receptors of the NR1H group. In addition, there are three heterologous NRs in amphioxus that could not be placed within the framework of the NR superfamily. Apart from these exceptions, there is usually one amphioxus NR or NR signaling coregulator for each paralogous group of two, three, or four human receptors suggesting that the ancestral chordate had a set of 22 different NRs plus one copy of each NR coregulator.