Gnathostome phylogenomics utilizing lungfish EST sequences

A detailed 3D MRI brain atlas of the African lungfish Protopterus annectens

The study of the brain by magnetic resonance imaging (MRI) in evolutionary analyses is still in its incipient stage, however, it is particularly useful as it allows us to analyze detailed anatomical images and compare brains of rare or otherwise inaccessible species, evolutionarily contextualizing possible differences, while at the same time being non-invasive. A good example is the lungfishes, sarcopterygians that are the closest living relatives of tetrapods and thus have an interesting phylogenetic position in the evolutionary conquest of the terrestrial environment. In the present study, we have developed a three-dimensional representation of the brain of the lungfish Protopterus annectens together with a rostrocaudal anatomical atlas. This methodological approach provides a clear delineation of the major brain subdivisions of this model and allows to measure both brain and ventricular volumes. Our results confirm that lungfish show neuroanatomical patterns reminiscent of those of extant basal sarcopterygians, with an evaginated telencephalon, and distinctive characters like a small optic tectum. These and additional characters uncover lungfish as a remarkable model to understand the origins of tetrapod diversity, indicating that their brain may contain significant clues to the characters of the brain of ancestral tetrapods.

Functional and Evolutionary Diversification of Otx2 and Crx in Vertebrate Retinal Photoreceptor and Bipolar Cell Development

Otx family homeoproteins Otx2 and Crx are expressed in photoreceptor precursor cells and bind to the common DNA-binding consensus sequence, but these two proteins have distinct functions in retinal development. To examine the functional substitutability of Otx2 and Crx, we generate knockin mouse lines in which Crx is replaced by Otx2 and vice versa. We find that Otx2 and Crx cannot be substituted in photoreceptor development. Subsequently, we investigate the function of Otx2 in photoreceptor and bipolar cell development. High Otx2 levels induce photoreceptor cell fate but not bipolar cell fate, whereas reduced Otx2 expression impairs bipolar cell maturation and survival. Furthermore, we identify Otx2 and Crx in the lamprey genome by using synteny analysis, suggesting that the last common ancestor of vertebrates possesses both Otx2 and Crx. We find that the retinal Otx2 expression pattern is different between lampreys and mice, suggesting that neofunctionalization of Otx2 occurred in the jawed vertebrate lineage.

Multimodal hypothalamo-hypophysial communication in the vertebrates

The vertebrate pituitary is arguably one of the most complex endocrine glands from the evolutionary, anatomical and functional perspectives. The pituitary plays a master role in endocrine physiology for the control of growth, metabolism, reproduction, water balance, and the stress response, among many other key processes. The synthesis and secretion of pituitary hormones are under the control of neurohormones produced by the hypothalamus. Under this conceptual framework, the communication between the hypophysiotropic brain and the pituitary gland is at the foundation of our understanding of endocrinology. The anatomy of the connections between the hypothalamus and the pituitary gland has been described in different vertebrate classes, revealing diverse modes of communication together with varying degrees of complexity. In this context, the evolution and variation in the neuronal, neurohemal, endocrine and paracrine modes will be reviewed in light of recent discoveries, and a re-evaluation of earlier observations. There appears to be three main hypothalamo-pituitary communication systems: 1. Diffusion, best exemplified by the agnathans; 2. Direct innervation of the adenohypophysis, which is most developed in teleost fish, and 3. The median eminence/portal blood vessel system, most conspicuously developed in tetrapods, showing also considerable variation between classes. Upon this basic classification, there exists various combinations possible, giving rise to taxon and species-specific, multimodal control over major physiological processes. Intrapituitary paracrine regulation and communication between folliculostellate cells and endocrine cells are additional processes of major importance. Thus, a more complex evolutionary picture of hypothalamo-hypophysial communication is emerging. There is currently little direct evidence to suggest which neuroendocrine genes may control the evolution of one communication system versus another. However, studies at the developmental and intergenerational timescales implicate several genes in the angiogenesis and axonal guidance pathways that may be important.

Aestivation Induces Changes in the mRNA Expression Levels and Protein Abundance of Two Isoforms of Urea Transporters in the Gills of the African Lungfish, Protopterus annectens

The African lungfish, Protopterus annectens, is ammonotelic in water despite being ureogenic. When it aestivates in mucus cocoon on land, ammonia is detoxified to urea. During the maintenance phase of aestivation, urea accumulates in the body, which is subsequently excreted upon arousal. Urea excretion involves urea transporters (UT/Ut). This study aimed to clone and sequence the ut isoforms from the gills of P. annectens, and to test the hypothesis that the mRNA and/or protein expression levels of ut/Ut isoforms could vary in the gills of P. annectens during the induction, maintenance, and arousal phases of aestivation. Two isoforms of ut, ut-a2a and ut-a2b, were obtained from the gills of P. annectens. ut-a2a consisted of 1227 bp and coded for 408 amino acids with an estimated molecular mass of 44.7 kDa, while ut-a2b consisted of 1392 bp and coded for 464 amino acids with an estimated molecular mass of 51.2 kDa. Ut-a2a and Ut-a2b of P. annectens had a closer phylogenetic relationship with Ut/UT of tetrapods than Ut of fishes. While the mRNA expression pattern of ut-a2a and ut-a2b across various tissues of P. annectens differed, the transcript levels of ut-a2a and ut-a2b in the gills were comparable, indicating that they might be equally important for branchial urea excretion during the initial arousal phase of aestivation. During the maintenance phase of aestivation, the transcript level of ut-a2a increased significantly, but the protein abundance of Ut-a2a remained unchanged in the gills of P. annectens. This could be an adaptive feature to prepare for an increase in the production of Ut-a2a upon arousal. Indeed, arousal led to a significant increase in the branchial Ut-a2a protein abundance. Although the transcript level of ut-a2b remained unchanged, there were significant increases in the protein abundance of Ut-a2b in the gills of P. annectens throughout the three phases of aestivation. The increase in the protein abundance of Ut-a2b during the maintenance phase could also be an adaptive feature to prepare for efficient urea excretion when water becomes available.

Molecular Characterization of Aquaporin 1 and Aquaporin 3 from the Gills of the African Lungfish, Protopterus annectens, and Changes in Their Branchial mRNA Expression Levels and Protein Abundance during Three Phases of Aestivation

African lungfishes can undergo long periods of aestivation on land during drought. During aestivation, lungfishes are confronted with desiccation and dehydration, and their gills become non-functional and covered with a thick layer of dried mucus. Aquaporins (Aqps) are a superfamily of integral membrane proteins which generally facilitate the permeation of water through plasma membranes. This study aimed to obtain the complete cDNA coding sequences of aqp1 and aqp3 from the gills of Protopterus annectens, and to determine their branchial mRNA and protein expression levels during the induction, maintenance and arousal phases of aestivation. Dendrogramic analyses of the deduced Aqp1 and Aqp3 amino acid sequences of P. annectens revealed their close relationships with those of Latimeria chalumnae and tetrapods. During the induction phase, there were significant decreases in the transcript levels of aqp1 and aqp3 in the gills of P. annectens, but the branchial Aqp1 and Aqp3 protein abundance remained unchanged. As changes in transcription might precede changes in translation, this could be regarded as an adaptive response to decrease the protein abundance of Aqp1 and Aqp3 in the subsequent maintenance phase of aestivation. As expected, the branchial transcript levels and protein abundance of aqp1/Aqp1 and aqp3/Aqp3 were significantly down-regulated during the maintenance phase, probably attributable to the shutdown of branchial functions and the cessation of volume regulation of branchial epithelial cells. Additionally, these changes could reduce the loss of water through branchial epithelial surfaces, supplementing the anti-desiccating property of the dried mucus. Upon arousal, it was essential for the lungfish to restore branchial functions. Indeed, the protein abundance of Aqp1 recovered partially, with complete recovery of mRNA expression level and protein abundance of Aqp3, in the gills of P. annectens after 3 days of arousal. These results provide insights into how P. annectens regulates branchial Aqp expression to cope with desiccation and rehydration during different phases of aestivation.

Genetic signatures of adaptation revealed from transcriptome sequencing of Arctic and red foxes

Background

The genus Vulpes (true foxes) comprises numerous species that inhabit a wide range of habitats and climatic conditions, including one species, the Arctic fox (Vulpes lagopus) which is adapted to the arctic region. A close relative to the Arctic fox, the red fox (Vulpes vulpes), occurs in subarctic to subtropical habitats. To study the genetic basis of their adaptations to different environments, transcriptome sequences from two Arctic foxes and one red fox individual were generated and analyzed for signatures of positive selection. In addition, the data allowed for a phylogenetic analysis and divergence time estimate between the two fox species.

Results

The de novo assembly of reads resulted in more than 160,000 contigs/transcripts per individual. Approximately 17,000 homologous genes were identified using human and the non-redundant databases. Positive selection analyses revealed several genes involved in various metabolic and molecular processes such as energy metabolism, cardiac gene regulation, apoptosis and blood coagulation to be under positive selection in foxes. Branch site tests identified four genes to be under positive selection in the Arctic fox transcriptome, two of which are fat metabolism genes. In the red fox transcriptome eight genes are under positive selection, including molecular process genes, notably genes involved in ATP metabolism. Analysis of the three transcriptomes and five Sanger re-sequenced genes in additional individuals identified a lower genetic variability within Arctic foxes compared to red foxes, which is consistent with distribution range differences and demographic responses to past climatic fluctuations. A phylogenomic analysis estimated that the Arctic and red fox lineages diverged about three million years ago.

Conclusions

Transcriptome data are an economic way to generate genomic resources for evolutionary studies. Despite not representing an entire genome, this transcriptome analysis identified numerous genes that are relevant to arctic adaptation in foxes. Similar to polar bears, fat metabolism seems to play a central role in adaptation of Arctic foxes to the cold climate, as has been identified in the polar bear, another arctic specialist.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1724-9) contains supplementary material, which is available to authorized users.

Molecular characterization of betaine-homocysteine methyltransferase 1 from the liver, and effects of aestivation on its expressions and homocysteine concentrations in the liver, kidney and muscle, of the African lungfish, Protopterus annectens

Homocysteine accumulation has numerous deleterious effects, and betaine-homocysteine S-methyltransferase (BHMT) catalyses the synthesis of methionine from homocysteine and betaine. This study aimed to determine homocysteine concentrations, and mRNA expression levels and protein abundances of bhmt1/Bhmt1 in the liver, kidney and muscle of the African lungfish, Protopterus annectens, during the induction (6 days), maintenance (6 months) or arousal (3 days after arousal) phase of aestivation. The homocysteine concentration decreased significantly in the liver of P. annectens after 6 days or 6 months of aestivation, but it returned to the control level upon arousal. By contrast, homocysteine concentrations in the kidney and muscle remained unchanged during the three phases of aestivation. The complete coding cDNA sequence of bhmt1 from P. annectens consisted of 1236 bp, coding for 412 amino acids. The Bhmt1 from P. annectens had a close phylogenetic relationship with those from tetrapods and Callorhinchus milii. The expression of bhmt1 was detected in multiple organs/tissues of P. annectens, and this is the first report on the expression of bhmt1/Bhmt1 in animal skeletal muscle. The mRNA and protein expression levels of bhmt1/Bhmt1 were up-regulated in the liver of P. annectens during the induction and maintenance phases of aestivation, possibly to regulate the hepatic homocysteine concentration. The significant increase in hepatic Bhmt1 protein abundance during the arousal phase could be a response to increased cellular methylation for the purpose of tissue reconstruction. Unlike the liver, Bhmt1 expression in the kidney and muscle of P. annectens was regulated translationally, and its up-regulation could be crucial to prevent homocysteine accumulation.

Aestivation induces changes in transcription and translation of coagulation factor II and fibrinogen gamma chain in the liver of the African lungfish, Protopterus annectens

This study aimed to sequence and characterize two pro-coagulant genes, coagulation factor II (f2) and fibrinogen gamma chain (fgg), from the liver of the African lungfish Protopterus annectens, and to determine their hepatic mRNA expression levels during three phases of aestivation. The protein abundances of F2 and Fgg in the liver and plasma were determined by immunoblotting. Results indicated that F2 and Fgg of P. annectens were phylogenetically closer to those of amphibians than those of teleosts. Three days of aestivation resulted in an up-regulation in the hepatic fgg mRNA expression level, while 6 days of aestivation led to a significant increase (3-fold) in the protein abundance of Fgg in the plasma. Hence, there could be an increase in the blood clotting ability in P. annectens during the induction phase of aestivation. By contrast, the blood clotting ability in P. annectens might be reduced in response to decreased blood flow and increased possibility of thrombosis during the maintenance phase of aestivation, as 6 months of aestivation led to significant decreases in mRNA expression levels of f2 and fgg in the liver. There could also be a decrease in the export of F2 and Fgg from the liver to the plasma so as to avert thrombosis. Upon 3-6 days of arousal from 6 months of aestivation, the protein abundances of F2 and Fgg recovered partially in the plasma of P. annectens, and a complete recovery of the transcription and translation of f2/F2 in the liver might occur only after refeeding.

Conserved localization of Pax6 and Pax7 transcripts in the brain of representatives of sarcopterygian vertebrates during development supports homologous brain regionalization

Many of the genes involved in brain patterning during development are highly conserved in vertebrates and similarities in their expression patterns help to recognize homologous cell types or brain regions. Among these genes, Pax6 and Pax7 are expressed in regionally restricted patterns in the brain and are essential for its development. In the present immunohistochemical study we analyzed the distribution of Pax6 and Pax7 cells in the brain of six representative species of tetrapods and lungfishes, the closest living relatives of tetrapods, at several developmental stages. The distribution patterns of these transcription factors were largely comparable across species. In all species only Pax6 was expressed in the telencephalon, including the olfactory bulbs, septum, striatum, and amygdaloid complex. In the diencephalon, Pax6 and Pax7 were distinct in the alar and basal parts, mainly in prosomeres 1 and 3. Pax7 specifically labeled cells in the optic tectum (superior colliculus) and Pax6, but not Pax7, cells were found in the tegmentum. Pax6 was found in most granule cells of the cerebellum and Pax7 labeling was detected in cells of the ventricular zone of the rostral alar plate and in migrated cells in the basal plate, including the griseum centrale and the interpeduncular nucleus. Caudally, Pax6 cells formed a column, whereas the ventricular zone of the alar plate expressed Pax7. Since the observed Pax6 and Pax7 expression patterns are largely conserved they can be used to identify subdivisions in the brain across vertebrates that are not clearly discernible with classical techniques.

Brain Na+/K+-ATPase α-subunit isoforms and aestivation in the African lungfish, Protopterus annectens

This study aimed to clone and sequence Na (+) / K (+)-ATPase (nka) α-subunit isoforms from, and to determine their mRNA expression levels and protein abundance in the brain of the African lungfish, Protopterus annectens during the induction, maintenance and arousal phases of aestivation in air. We obtained the full cDNA sequences of nkaα1, nkaα2 and nkaα3 from the brain of P. annectens. Phylogenetic analysis of their deduced amino acid sequences revealed that they are closer to the corresponding NKA α-subunits of tetrapods than to those of fishes. The mRNA expression of these three nkaα isoforms showed differential changes in the brain of P. annectens during the three phases of aestivation. After 12 days of aestivation, there was a significant increase in the protein abundance of Nkaα1 in the brain of P. annectens. This could be an important response to maintain cellular Na(+) and K(+) concentrations and regulate cell volume during the early maintenance phase of aestivation. On the other hand, the mRNA expression of nkaα2 decreased significantly in the brain of P. annectens after 6 months of aestivation, which could be a result of a suppression of transcriptional activities to reduce energy expenditure. The down-regulation of mRNA expression of nkaα1, nkaα2 and nkaα3 and the overall protein abundance of Nka α-subunit isoforms in the brain of P. annectens after 1 day of arousal from 6 months of aestivation were novel observations, and it could be an adaptive response to restore blood pressure and/or to prevent brain oedema.

Intraclass diversification of immunoglobulin heavy chain genes in the African lungfish

Lungfish (Dipnoi) are the closest living relatives to tetrapods, and they represent the transition from water to land during vertebrate evolution. Lungfish are armed with immunoglobulins (Igs), one of the hallmarks of the adaptive immune system of jawed vertebrates, but only three Ig forms have been characterized in Dipnoi to date. We report here a new diversity of Ig molecules in two African lungfish species (Protopterus dolloi and Protopterus annectens). The African lungfish Igs consist of three IgMs, two IgWs, three IgNs, and an IgQ, where both IgN and IgQ originated evidently from the IgW lineage. Our data also suggest that the IgH genes in the lungfish are organized in a transiting form from clusters (IgH loci in cartilaginous fish) to a translocon configuration (IgH locus in tetrapods). We propose that the intraclass diversification of the two primordial gnathostome Ig classes (IgM and IgW) as well as acquisition of new isotypes (IgN and IgQ) has allowed lungfish to acquire a complex and functionally diverse Ig repertoire to fight a variety of microorganisms. Furthermore, our results support the idea that "tetrapod-specific" Ig classes did not evolve until the vertebrate adaptation to land was completed ~360 million years ago.

Elephant shark genome provides unique insights into gnathostome evolution

The emergence of jawed vertebrates (gnathostomes) from jawless vertebrates was accompanied by major morphological and physiological innovations, such as hinged jaws, paired fins and immunoglobulin-based adaptive immunity. Gnathostomes subsequently diverged into two groups, the cartilaginous fishes and the bony vertebrates. Here we report the whole-genome analysis of a cartilaginous fish, the elephant shark (Callorhinchus milii). We find that the C. milii genome is the slowest evolving of all known vertebrates, including the ‘living fossil’ coelacanth, and features extensive synteny conservation with tetrapod genomes, making it a good model for comparative analyses of gnathostome genomes. Our functional studies suggest that the lack of genes encoding secreted calcium-binding phosphoproteins in cartilaginous fishes explains the absence of bone in their endoskeleton. Furthermore, the adaptive immune system of cartilaginous fishes is unusual: it lacks the canonical CD4 co-receptor and most transcription factors, cytokines and cytokine receptors related to the CD4 lineage, despite the presence of polymorphic major histocompatibility complex class II molecules. It thus presents a new model for understanding the origin of adaptive immunity.

Whole-genome analysis of the elephant shark, a cartilaginous fish, shows that it is the slowest evolving of all known vertebrates, lacks critical bone formation genes and has an unusual adaptive immune system.

The elephant shark (Callorhinchus milii) is a cartilaginous fish native to the temperate waters off southern Australia and New Zealand, living at depths of 200 to 500 metres and migrating into shallow waters during spring for breeding. The genome sequence is published in this issue of Nature. Comparison with other vertebrate genomes shows that it is the slowest evolving genome of all known vertebrates — coelacanth included. Genome analysis points to an unusual adaptive immune system lacking the CD4 receptor and some associated cytokines, indicating that cartilaginous fishes possess a primordial gnathostome adaptive immune system. Also absent are genes encoding secreted calcium-binding phosphoproteins, in line with the absence of bone in cartilaginous fish.

Drastic neofunctionalization associated with evolution of the timezyme AANAT 500 Mya

Melatonin (N-acetyl-5-methoxytrypamine) is the vertebrate hormone of the night: circulating levels at night are markedly higher than day levels. This increase is driven by precisely regulated increases in acetylation of serotonin in the pineal gland by arylalkylamine N-acetyltransferase (AANAT), the penultimate enzyme in the synthesis of melatonin. This unique essential role of AANAT in vertebrate timekeeping is recognized by the moniker the timezyme. AANAT is also found in the retina, where melatonin is thought to play a paracrine role. Here, we focused on the evolution of AANAT in early vertebrates. AANATs from Agnathans (lamprey) and Chondrichthyes (catshark and elephant shark) were cloned, and it was found that pineal glands and retinas from these groups express a form of AANAT that is compositionally, biochemically, and kinetically similar to AANATs found in bony vertebrates (VT-AANAT). Examination of the available genomes indicates that VT-AANAT is absent from other forms of life, including the Cephalochordate amphioxus. Phylogenetic analysis and evolutionary rate estimation indicate that VT-AANAT evolved from the nonvertebrate form of AANAT after the Cephalochordate-Vertebrate split over one-half billion years ago. The emergence of VT-AANAT apparently involved a dramatic acceleration of evolution that accompanied neofunctionalization after a duplication of the nonvertebrate AANAT gene. This scenario is consistent with the hypotheses that the advent of VT-AANAT contributed to the evolution of the pineal gland and lateral eyes from a common ancestral photodetector and that it was not a posthoc recruitment.

On the measurement of ecological novelty: scale-eating pupfish are separated by 168 my from other scale-eating fishes

The colonization of new adaptive zones is widely recognized as one of the hallmarks of adaptive radiation. However, the adoption of novel resources during this process is rarely distinguished from phenotypic change because morphology is a common proxy for ecology. How can we quantify ecological novelty independent of phenotype? Our study is split into two parts: we first document a remarkable example of ecological novelty, scale-eating (lepidophagy), within a rapidly-evolving adaptive radiation of Cyprinodon pupfishes on San Salvador Island, Bahamas. This specialized predatory niche is known in several other fish groups, but is not found elsewhere among the 1,500 species of atherinomorphs. Second, we quantify this ecological novelty by measuring the time-calibrated phylogenetic distance in years to the most closely-related species with convergent ecology. We find that scale-eating pupfish are separated by 168 million years of evolution from the nearest scale-eating fish. We apply this approach to a variety of examples and highlight the frequent decoupling of ecological novelty from phenotypic divergence. We observe that novel ecology is not always tightly correlated with rates of phenotypic or species diversification, particularly within recent adaptive radiations, necessitating the use of additional measures of ecological novelty independent of phenotype.

Early vertebrate origin of melanocortin 2 receptor accessory proteins (MRAPs)

The melanocortin 2 receptor (MC2R) accessory proteins, MRAP, along with its homolog, MRAP2, are two among a growing number of G protein-coupled receptor accessory proteins that have been identified in recent years. These proteins interact directly with MC2R and are essential for trafficking of this receptor from the endoplasmic reticulum to the cell surface, where it mediates the effects of ACTH. lthough earlier studies have identified MRAP and MRAP2 subtypes in distant species, an overall evolutionary analysis of these families is still missing. Here, we performed a comprehensive evolutionary analysis of the MRAP and MRAP2 homologs based on whole genome sequences. We systematically mined and analyzed the genomes of metazoans to identify these genes. Overall, we identified 70 sequences of MRAP and MRAP2 from 44 species belonging to several vertebrate lineages, including at least 40 new sequences previously not reported in the literature. Herein, we provide evidence that MRAP2 is likely to be the ancestor of the MRAP family because MRAP2-like protein, but not MRAP, was identified in Petromyzon marinus (sea lamprey), which belong to an ancient basal vertebrate lineage. Later in vertebrate evolution, MRAP2 duplicated and gave rise to MRAP in an event before the emergence of actinopterygii (ray-finned fishes). However, we observed losses of MRAP in sarcopterygii (lobe-finned fish), amphibians and reptiles while both subtypes are present in chicken and most mammals studied. Synteny analysis showed a conserved synteny within same lineages and an inversion of gene order between lineages. An evolutionary rate shift analysis indicated that these genes were under high purifying selection. Overall, this study provides a comprehensive analysis of the evolution and gene repertoire of MRAP and MRAP2.

Basal jawed vertebrate phylogenomics using transcriptomic data from Solexa sequencing

The traditionally accepted relationships among basal jawed vertebrates have been challenged by some molecular phylogenetic analyses based on mitochondrial sequences. Those studies split extant gnathostomes into two monophyletic groups: tetrapods and piscine branch, including Chondrichthyes, Actinopterygii and sarcopterygian fishes. Lungfish and bichir are found in a basal position on the piscine branch. Based on transcriptomes of an armored bichir (Polypterus delhezi) and an African lungfish (Protopterus sp.) we generated, expressed sequences and whole genome sequences available from public databases, we obtained 111 genes to reconstruct the phylogenetic tree of basal jawed vertebrates and estimated their times of divergence. Our phylogenomic study supports the traditional relationship. We found that gnathostomes are divided into Chondrichthyes and the Osteichthyes, both with 100% support values (posterior probabilities and bootstrap values). Chimaeras were found to have a basal position among cartilaginous fishes with a 100% support value. Osteichthyes were divided into Actinopterygii and Sarcopterygii with 100% support value. Lungfish and tetrapods form a monophyletic group with 100% posterior probability. Bichir and two teleost species form a monophyletic group with 100% support value. The previous tree, based on mitochondrial data, was significantly rejected by an approximately unbiased test (AU test, p = 0). The time of divergence between lungfish and tetrapods was estimated to be 391.8 Ma and the divergence of bichir from pufferfish and medaka was estimated to be 330.6 Ma. These estimates closely match the fossil record. In conclusion, our phylogenomic study successfully resolved the relationship of basal jawed vertebrates based on transtriptomes, EST and whole genome sequences.

Characterization of the bed nucleus of the stria terminalis in the forebrain of anuran amphibians

Major common features have been reported for the organization of the basal telencephalon in amniotes, and most characteristics were thought to be acquired in the transition from anamniotes to amniotes. However, gene expression, neurochemical, and hodological data obtained for the basal ganglia and septal and amygdaloid complexes in amphibians (anamniotic tetrapods) have strengthened the idea of a conserved organization in tetrapods. A poorly characterized region in the forebrain of amniotes has been the bed nucleus of the stria terminalis (BST), but numerous recent investigations have characterized it as a member of the extended amygdala. Our study analyzes the main features of the BST in anuran amphibians to establish putative homologies with amniotes. Gene expression patterns during development identified the anuran BST as a subpallial, nonstriatal territory. The BST shows Nkx2.1 and Lhx7 expression and contains an Islet1-positive cell subpopulation derived from the lateral ganglionic eminence. Immunohistochemistry for diverse peptides and neurotransmitters revealed that the distinct chemoarchitecture of the BST is strongly conserved among tetrapods. In vitro tracing techniques with dextran amines revealed important connections between the BST and the central and medial amygdala, septal territories, medial pallium, preoptic area, lateral hypothalamus, thalamus, and prethalamus. The BST receives dopaminergic projections from the ventral tegmental area and is connected with the laterodorsal tegmental nucleus and the rostral raphe in the brainstem. All these data suggest that the anuran BST shares many features with its counterpart in amniotes and belongs to a basal continuum, likely controlling similar reflexes, reponses, and behaviors in tetrapods.

The oxytocin/vasopressin receptor family has at least five members in the gnathostome lineage, inclucing two distinct V2 subtypes

The vertebrate oxytocin and vasopressin receptors form a family of G-protein-coupled receptors (GPCRs) that mediate a large variety of functions, including social behavior and the regulation of blood pressure, water balance and reproduction. In mammals four family members have been identified, three of which respond to vasopressin (VP) named V1A, V1B and V2, and one of which is activated by oxytocin (OT), called the OT receptor. Four receptors have been identified in chicken as well, but these have received different names. Until recently only V1-type receptors have been described in several species of teleost fishes. We have identified family members in several gnathostome genomes and performed phylogenetic analyses to classify OT/VP-receptors across species and determine orthology relationships. Our phylogenetic tree identifies five distinct ancestral gnathostome receptor subtypes in the OT/VP receptor family: V1A, V1B, V2A, V2B and OT receptors. The existence of distinct V2A and V2B receptors has not been previously recognized. We have found these two subtypes in all examined teleost genomes as well as in available frog and lizard genomes and conclude that the V2A-type is orthologous to mammalian V2 receptors whereas the V2B-type is orthologous to avian V2 receptors. Some teleost fishes have acquired additional and more recent gene duplicates with up to eight receptor family members. Thus, this analysis reveals an unprecedented complexity in the gnathostome repertoire of OT/VP receptors, opening interesting research avenues regarding functions such as regulation of water balance, reproduction and behavior, particularly in reptiles, amphibians, teleost fishes and cartilaginous fishes.

Overview of the transcriptome profiles identified in hagfish, shark, and bichir: current issues arising from some nonmodel vertebrate taxa

Because of their crucial phylogenetic positions, hagfishes, sharks, and bichirs are recognized as key taxa in our understanding of vertebrate evolution. The expression patterns of the regulatory genes involved in developmental patterning have been analyzed in the context of evolutionary developmental studies. However, in a survey of public sequence databases, we found that the large-scale sequence data for these taxa are still limited. To address this deficit, we used conventional Sanger DNA sequencing and a next-generation sequencing technology based on 454 GS FLX sequencing to obtain expressed sequence tags (ESTs) of the Japanese inshore hagfish (Eptatretus burgeri; 161,482 ESTs), cloudy catshark (Scyliorhinus torazame; 165,819 ESTs), and gray bichir (Polypterus senegalus; 34,336 ESTs). We deposited the ESTs in a newly constructed database, designated the "Vertebrate TimeCapsule." The ESTs include sequences from genes that can be effectively used in evolutionary developmental studies; for instance, several encode cartilaginous extracellular matrix proteins, which are central to an understanding of the ways in which evolutionary processes affected the skeletal elements, whereas others encode regulatory genes involved in craniofacial development and early embryogenesis. Here, we discuss how hagfishes, sharks, and bichirs contribute to our understanding of vertebrate evolution, we review the current status of the publicly available sequence data for these three taxa, and we introduce our EST projects and newly developed database.

Phylogenetic analysis of kindlins suggests subfunctionalization of an ancestral unduplicated kindlin into three paralogs in vertebrates

Kindlin proteins represent a newly discovered family of evolutionarily conserved FERM domain-containing proteins. This family includes three highly conserved proteins: Kindlin-1, Kindlin-2 and Kindlin-3. All three Kindlin proteins are associated with focal adhesions and are involved in integrin activation. The FERM domain of each Kindlin is bipartite and plays a key role in integrin activation. We herein explore for the first time the evolutionary history of these proteins. The phylogeny of the Kindlins suggests a single ancestral Kindlin protein present in even the earliest metazoan ie, hydra. This protein then underwent duplication events in insects and also experienced genome duplication in vertebrates, leading to the Kindlin family. A comparative study of the Kindlin paralogs showed that Kindlin-2 is the slowest evolving protein among the three family members. The analysis of synonymous and non-synonymous substitutions in orthologous Kindlin sequences in different species showed that all three Kindlins have been evolving under the influence of purifying selection. The expression pattern of Kindlins along with phylogenetic studies supports the subfunctionalization model of gene duplication.

The non-evaginated secondary prosencephalon of vertebrates

The secondary prosencephalon (telencephalon plus hypothalamus) is probably the most complex area of the brain, with complicated patterning specifications. As yet, no prosomeric subdivisions have been reported and only distinct histogenetic territories have been recognized. In the present comparative study we analyzed cross-correlated expression maps in the non-evaginated territories of the secondary prosencephalon in different vertebrates throughout development, to assess the existence of comparable divisions and subdivisions in the different groups. Each division is characterized by expression of a unique combination of developmental regulatory genes, and each appears to represent a self-regulated and topologically constant histogenetic brain compartment that gives rise to a specific cell group. The non-evaginated area of the telencephalon corresponds to the preoptic region, whereas the hypothalamus, topologically rostral to the diencephalic prethalamus, includes basal (mammillary and tuberal) and alar (paraventricular and suprachiasmatic) parts. This complex area is specified by a cascade of transcription factors, among which the Dlx family members and Nkx2.1 are essential for the correct development. The only exception is found in the subdivision named termed the supraoptoparaventricular area, in which the transcription factor Orthopedia is essential in restricting the fate of multiple categories of neuroendocrine neurons, in the absence of the Dlx/Nkx2.1 combination. Our analysis, based on own data and published results by other researchers, suggests that common features are shared at least by all tetrapods and, therefore, they most likely were present in the stem tetrapods. The available data for agnathans (lampreys) and other fish groups indicate that not all subdivisions of the secondary prosencephalon were present at the origin of vertebrates, raising important questions about their evolution.

Immunohistochemical localization of calbindin-D28k and calretinin in the spinal cord of lungfishes

A common pattern of distribution of neurons and fibers containing the calcium-binding proteins calbindin-D28k (CB) and calretinin (CR) in the spinal cord of terrestrial vertebrates has been recently demonstrated. Lungfishes are considered the closest living relatives of tetrapods, but practically no experimental data exist on the organization of their spinal cord. By means of immunohistochemical techniques, the localization of CB and CR was investigated in the spinal cord of the African (Protopterus dolloi) and Australian (Neoceratodus forsteri) lungfishes. Abundant cell bodies and fibers immunoreactive for either CB or CR were widely distributed throughout the spinal cord. A large population of immunoreactive cells was found in the dorsal column of the gray matter in both species, and abundant cells were distributed in the lateral and ventral columns. Ventrolateral motoneurons and multipolar cells were only intensely CB and CR immunoreactive in Neoceratodus. For the most part, separate cell populations contained either CB or CR, but a small subset of dorsally located neurons contained both in the two lungfishes. Colocalization was found in motoneurons and in ventrolaterally located cells only in Neoceratodus. Fiber labeling showed a predominance of CR-containing axons in the lateral and ventral funiculi of presumed supraspinal origin. These results show that lung-fishes and tetrapods have many features in common, suggesting that primitive anatomical, and likely functional, organization of the spinal cord of tetrapods is present in lungfishes.

Lungfishes, like tetrapods, possess a vomeronasal system

The vomeronasal system (VNS) is an accessory olfactory system that in tetrapod vertebrates is composed of specific receptor neurons in the nasal organ and a set of centers in the forebrain that receive and relay the information consecutively towards the hypothalamus. Thus, only in tetrapods the VNS comprises a discrete vomeronasal (Jacobson's) organ, which contains receptor cells that are morphologically distinct from those of the olfactory epithelium and use different transduction mechanisms. The axons of the vomeronasal receptors in tetrapods project to the accessory olfactory bulb (AOB) in the rostral telencephalon. Secondary vomeronasal connections exist through the medial amygdala to the hypothalamus. Currently, the lungfishes are considered the closest living relatives of tetrapods. Here we show that the African lungfish, Protopterus dolloi, has epithelial crypts at the base of the lamellae of the olfactory epithelium that express markers of the vomeronasal receptors in tetrapods. The projections of these crypts allow us to identify an AOB on the lateral margin of the main olfactory bulb. The projections of this AOB reach a region that is topologically, hodologically, and immunohistochemically identical to the medial amygdala and could represent its homolog. Neurons of this putative medial amygdala were demonstrated to project to the lateral hypothalamus, as they do in tetrapods. All these features that lungfishes share with tetrapods indicate that lungfishes have the complete set of brain centers and connections involved in processing vomeronasal information and that these features were already present in the last common ancestor of lungfishes and tetrapods.

Sonic hedgehog expression during Xenopus laevis forebrain development

We have analyzed the developing expression pattern of x-Shh in the Xenopus forebrain, interpreting the results within the framework of the neuromeric model to assess evolutionary trends and clues. To achieve this goal, we have characterized phenotypically the developing x-Shh expressing forebrain subdivisions and neurons by means of the combination of in situ hybridization for x-Shh and immunohistochemistry for the detection of forebrain essential regulators and markers, such as the homeodomain transcription factors Islet 1, Orthopedia, NKX2.1 and NKX2.2 and tyrosine hydroxylase. Substantial evidence was found for x-Shh expression in the telencephalic commissural preoptic area and this is strongly correlated with the presence of a pallidum and/or a basal telencephalic cholinergic system. In the diencephalon, x-Shh was demonstrated in the zona limitans intrathalamica and the x-Shh expressing cells were extended into the prethalamus. Throughout development and in the adult hypothalamic x-Shh expression was strong in basal regions but, in addition, in the alar suprachiasmatic region. The findings obtained in the forebrain of Xenopus revealed a largely conserved pattern of Shh expression among tetrapods. However, interesting differences were also noted that could be related to evolutive changes in forebrain organization.

Systematics of the subfamily Danioninae (Teleostei: Cypriniformes: Cyprinidae)

The members of the cyprinid subfamily Danioninae form a diverse and scientifically important group of fishes, which includes the zebrafish, Danio rerio. The diversity of this assemblage has attracted much scientific interest but its monophyly and the relationships among its members are poorly understood. The phylogenetic relationships of the Danioninae are examined herein using sequence data from mitochondrial cytochrome b, mitochondrial cytochrome c oxidase I, nuclear opsin, and nuclear recombination activating gene 1. A combined data matrix of 4117 bp for 270 taxa was compiled and analyzed. The resulting topology supports some conclusions drawn by recent studies on the group and certain portions of the traditional classification, but our results also contradict key aspects of the traditional classification. The subfamily Danioninae is not monophyletic, with putative members scattered throughout Cyprinidae. Therefore, we restrict Danioninae to the monophyletic group that includes the following genera: Amblypharyngodon, Barilius, Cabdio, Chela, Chelaethiops, Danio, Danionella, Devario (including Inlecypris), Esomus, Horadandia, Laubuca, Leptocypris, Luciosoma, Malayochela, Microdevario, Microrasbora, Nematabramis, Neobola, Opsaridium, Opsarius, Paedocypris, Pectenocypris, Raiamas, Rasbora (including Boraras and Trigonostigma), Rasboroides, Salmostoma, Securicula, and Sundadanio. This Danioninae sensu stricto is divided into three major lineages, the tribes Chedrini, Danionini, and Rasborini, where Chedrini is sister to a Danionini-Rasborini clade. Each of these tribes is monophyletic, following the restriction of Danioninae. The tribe Chedrini includes a clade of exclusively African species and contains several genera of uncertain monophyly (Opsarius, Raiamas, Salmostoma). Within the tribe Rasborini, the species-rich genus Rasbora is rendered non-monophyletic by the placement of two monophyletic genera, Boraras and Trigonostigma, hence we synonymize those two genera with Rasbora. In the tribe Danionini, the miniature genus Danionella is recovered as the sister group of Danio, with D. nigrofasciatus sister to D. rerio.

Phylogenomic analyses reveal the evolutionary origin of the inhibin alpha-subunit, a unique TGFbeta superfamily antagonist

Transforming growth factor-beta (TGFβ) homologues form a diverse superfamily that arose early in animal evolution and control cellular function through membrane-spanning, conserved serine-threonine kinases (RII and RI receptors). Activin and inhibin are related dimers within the TGFβ superfamily that share a common β-subunit. The evolution of the inhibin α-subunit created the only antagonist within the TGFβ superfamily and the only member known to act as an endocrine hormone. This hormone introduced a new level of complexity and control to vertebrate reproductive function. The novel functions of the inhibin α-subunit appear to reflect specific insertion-deletion changes within the inhibin β-subunit that occurred during evolution. Using phylogenomic analysis, we correlated specific insertions with the acquisition of distinct functions that underlie the phenotypic complexity of vertebrate reproductive processes. This phylogenomic approach presents a new way of understanding the structure-function relationships between inhibin, activin, and the larger TGFβ superfamily.

Network dynamics of eukaryotic LTR retroelements beyond phylogenetic trees

BACKGROUND

Sequencing projects have allowed diverse retroviruses and LTR retrotransposons from different eukaryotic organisms to be characterized. It is known that retroviruses and other retro-transcribing viruses evolve from LTR retrotransposons and that this whole system clusters into five families: Ty3/Gypsy, Retroviridae, Ty1/Copia, Bel/Pao and Caulimoviridae. Phylogenetic analyses usually show that these split into multiple distinct lineages but what is yet to be understood is how deep evolution occurred in this system.

RESULTS

We combined phylogenetic and graph analyses to investigate the history of LTR retroelements both as a tree and as a network. We used 268 non-redundant LTR retroelements, many of them introduced for the first time in this work, to elucidate all possible LTR retroelement phylogenetic patterns. These were superimposed over the tree of eukaryotes to investigate the dynamics of the system, at distinct evolutionary times. Next, we investigated phenotypic features such as duplication and variability of amino acid motifs, and several differences in genomic ORF organization. Using this information we characterized eight reticulate evolution markers to construct phenotypic network models.

CONCLUSION

The evolutionary history of LTR retroelements can be traced as a time-evolving network that depends on phylogenetic patterns, epigenetic host-factors and phenotypic plasticity. The Ty1/Copia and the Ty3/Gypsy families represent the oldest patterns in this network that we found mimics eukaryotic macroevolution. The emergence of the Bel/Pao, Retroviridae and Caulimoviridae families in this network can be related with distinct inflations of the Ty3/Gypsy family, at distinct evolutionary times. This suggests that Ty3/Gypsy ancestors diversified much more than their Ty1/Copia counterparts, at distinct geological eras. Consistent with the principle of preferential attachment, the connectivities among phenotypic markers, taken as network-represented combinations, are power-law distributed. This evidences an inflationary mode of evolution where the system diversity; 1) expands continuously alternating vertical and gradual processes of phylogenetic divergence with episodes of modular, saltatory and reticulate evolution; 2) is governed by the intrinsic capability of distinct LTR retroelement host-communities to self-organize their phenotypes according to emergent laws characteristic of complex systems.

REVIEWERS

This article was reviewed by Eugene V. Koonin, Eric Bapteste, and Enmanuelle Lerat (nominated by King Jordan).

Whole genome duplications and expansion of the vertebrate GATA transcription factor gene family

Background

GATA transcription factors influence many developmental processes, including the specification of embryonic germ layers. The GATA gene family has significantly expanded in many animal lineages: whereas diverse cnidarians have only one GATA transcription factor, six GATA genes have been identified in many vertebrates, five in many insects, and eleven to thirteen in Caenorhabditis nematodes. All bilaterian animal genomes have at least one member each of two classes, GATA123 and GATA456.

Results

We have identified one GATA123 gene and one GATA456 gene from the genomic sequence of two invertebrate deuterostomes, a cephalochordate (Branchiostoma floridae) and a hemichordate (Saccoglossus kowalevskii). We also have confirmed the presence of six GATA genes in all vertebrate genomes, as well as additional GATA genes in teleost fish. Analyses of conserved sequence motifs and of changes to the exon-intron structure, and molecular phylogenetic analyses of these deuterostome GATA genes support their origin from two ancestral deuterostome genes, one GATA 123 and one GATA456. Comparison of the conserved genomic organization across vertebrates identified eighteen paralogous gene families linked to multiple vertebrate GATA genes (GATA paralogons), providing the strongest evidence yet for expansion of vertebrate GATA gene families via genome duplication events.

Conclusion

From our analysis, we infer the evolutionary birth order and relationships among vertebrate GATA transcription factors, and define their expansion via multiple rounds of whole genome duplication events. As the genomes of four independent invertebrate deuterostome lineages contain single copy GATA123 and GATA456 genes, we infer that the 0R (pre-genome duplication) invertebrate deuterostome ancestor also had two GATA genes, one of each class. Synteny analyses identify duplications of paralogous chromosomal regions (paralogons), from single ancestral vertebrate GATA123 and GATA456 chromosomes to four paralogons after the first round of vertebrate genome duplication, to seven paralogons after the second round of vertebrate genome duplication, and to fourteen paralogons after the fish-specific 3R genome duplication. The evolutionary analysis of GATA gene origins and relationships may inform understanding vertebrate GATA factor redundancies and specializations.