Comparative genomics of the vertebrate insulin/TOR signal transduction pathway: A network-level analysis of selective pressures

Complexity of biological function relies on large networks of interacting molecules. However, the evolutionary properties of these networks are not fully understood. It has been shown that selective pressures depend on the position of genes in the network. We have previously shown that in the Drosophila insulin/target of rapamycin (TOR) signal transduction pathway there is a correlation between the pathway position and the strength of purifying selection, with the downstream genes being most constrained. In this study, we investigated the evolutionary dynamics of this well-characterized pathway in vertebrates. More specifically, we determined the impact of natural selection on the evolution of 72 genes of this pathway. We found that in vertebrates there is a similar gradient of selective constraint in the insulin/TOR pathway to that found in Drosophila. This feature is neither the result of a polarity in the impact of positive selection nor of a series of factors affecting selective constraint levels (gene expression level and breadth, codon bias, protein length, and connectivity). We also found that pathway genes encoding physically interacting proteins tend to evolve under similar selective constraints. The results indicate that the architecture of the vertebrate insulin/TOR pathway constrains the molecular evolution of its components. Therefore, the polarity detected in Drosophila is neither specific nor incidental of this genus. Hence, although the underlying biological mechanisms remain unclear, these may be similar in both vertebrates and Drosophila.

Potential hosts number in cuckoo bees (Psithyrus subgen.) increases toward higher elevations

In severe and variable conditions, specialized resource selection strategies should be less fre‐ quent because extinction risks increase for species that depend on a single and unstable resource. Psithyrus (Bombus subgenus Psithyrus) are bumblebee parasites that usurp Bombus nests and display inter‐specific variation in the number of hosts they parasitize. Using a phylogenetic comparative frame‐ work, we show that Psithyrus species at higher elevations display a higher number of hosts species com‐ pared with species restricted to lower elevations. Species inhabiting high elevations also cover a larger temperature range, suggesting that species able to occur in colder conditions may benefit from recruit‐ ment from populations occurring in warmer conditions. Our results provide evidence for an 'altitudinal niche breadth hypothesis' in parasitic species, showing a decrease in the parasites' specialization along the elevational gradient, and also suggesting that Rapoport's rule might apply to Psithyrus.

Molecular signaling in zebrafish development and the vertebrate phylotypic period

During development vertebrate embryos pass through a stage where their morphology is most conserved between species, the phylotypic period (approximately the pharyngula). To explain the resistance to evolutionary changes of this period, one hypothesis suggests that it is characterized by a high level of interactions. Based on this hypothesis, we examined protein-protein interactions, signal transduction cascades and miRNAs over the course of zebrafish development, and the conservation of expression of these genes in mouse development. We also investigated the characteristics of genes highly expressed before or during the presumed phylotypic period. We show that while there is a high diversity of interactions during the phylotypic period (protein-DNA, RNA-RNA, cell-cell, and between tissues), which is well conserved with mouse, there is no clear difference with later, more morphologically divergent, stages. We propose that the phylotypic period may rather be the expression at the morphological level of strong conservation of molecular processes earlier in development.

Plants and tortoises: mutations in the Arabidopsis jasmonate pathway increase feeding in a vertebrate herbivore.

Photosynthetic tissues, the major food source of many invertebrates and vertebrates, are well defended. Many defence traits in leaves are controlled via the jasmonate signalling pathway in which jasmonate acts as a hormone by binding to a receptor to activate responses that lead to increased resistance to invertebrate folivores. We predicted that mutations in jasmonate synthesis might also increase the vulnerability of leaves to vertebrate folivores and tested this hypothesis using the Eastern Hermann's tortoise (Eurotestudo boettgeri) and an Arabidopsis thaliana (Brassicaceae) allene oxide synthase (aos) mutant unable to synthesize jasmonate. Tortoises preferred the aos mutant over the wild type (WT). Based on these results, we then investigated the effect of mutating jasmonate perception using a segregating population of the recessive A. thaliana jasmonate receptor mutant coronatine insensitive1-1 (coi1-1). Genotyping of these plants after tortoise feeding revealed that the homozygous coi1-1 receptor mutant was consumed more readily than the heterozygous mutant or the WT. Therefore, the plant's ability to synthesize or perceive jasmonate reduces feeding by a vertebrate herbivore. We also tested whether or not tortoise feeding behaviour was influenced by glucosinolates, the principal defence chemicals in Arabidopsis leaves with known roles in defence against many generalist insects. However, in contrast to what has been observed with such insects, leaves in which the levels of these compounds were reduced genetically were consumed at a similar rate to those of the WT.

Simultaneous Mendelian and clonal genome transmission in a sexually reproducing, all-triploid vertebrate.

Meiosis in triploids faces the seemingly insuperable difficulty of dividing an odd number of chromosome sets by two. Triploid vertebrates usually circumvent this problem through either asexuality or some forms of hybridogenesis, including meiotic hybridogenesis that involve a reproductive community of different ploidy levels and genome composition. Batura toads (Bufo baturae; 3n = 33 chromosomes), however, present an all-triploid sexual reproduction. This hybrid species has two genome copies carrying a nucleolus-organizing region (NOR+) on chromosome 6, and a third copy without it (NOR-). Males only produce haploid NOR+ sperm, while ova are diploid, containing one NOR+ and one NOR- set. Here, we conduct sibship analyses with co-dominant microsatellite markers so as (i) to confirm the purely clonal and maternal transmission of the NOR- set, and (ii) to demonstrate Mendelian segregation and recombination of the NOR+ sets in both sexes. This new reproductive mode in vertebrates ('pre-equalizing hybrid meiosis') offers an ideal opportunity to study the evolution of non-recombining genomes. Elucidating the mechanisms that allow simultaneous transmission of two genomes, one of Mendelian, the other of clonal inheritance, might shed light on the general processes that regulate meiosis in vertebrates.

High prevalence and lineage diversity of avian malaria in wild populations of great tits (Parus major) and mosquitoes (Culex pipiens).

Avian malaria studies have taken a prominent place in different aspects of evolutionary ecology. Despite a recent interest in the role of vectors within the complex interaction system of the malaria parasite, they have largely been ignored in most epidemiological studies. Epidemiology of the disease is however strongly related to the vector's ecology and behaviour, and there is a need for basic investigations to obtain a better picture of the natural associations between Plasmodium lineages, vector species and bird hosts. The aim of the present study was to identify the mosquito species involved in the transmission of the haemosporidian parasites Plasmodium spp. in two wild populations of breeding great tits (Parus major) in western Switzerland. Additionally, we compared Plasmodium lineages, based on mitochondrial DNA cytochrome b sequences, between the vertebrate and dipteran hosts, and evaluated the prevalence of the parasite in the mosquito populations. Plasmodium spp. were detected in Culex pipiens only, with an overall 6.6% prevalence. Among the six cytochrome b lineages of Plasmodium identified in the mosquitoes, three were also present in great tits. The results provide evidence for the first time that C. pipiens can act as a natural vector of avian malaria in Europe and yield baseline data for future research on the epidemiology of avian malaria in European countries.

Geochemical study of vertebrate fossils from the Upper Cretaceous (Santonian) Csehbanya Formation (Hungary): Evidence for a freshwater habitat of mosasaurs and pycnodont fish

The diverse vertebrate remains from the Upper Cretaceous freshwater settings at Iharkut, Hungary, contain two fossil groups, Pycnodontiformes fish and Mosasauridae that are almost exclusively known from marine palaeo-environments. Hence, their appearance in alluvial sediments is very unusual. Trace element and isotope compositions of the remains have been analyzed to investigate the taphonomy and the ecological differences among the different fossil groups present at Iharkut. All examined fossils have undergone post-depositional diagenetic alteration, which resulted in high concentrations of REE, U, and Fe, together with almost complete homogenization of delta(18)O(CO3) values. Similar REE patterns in different fossils suggest a common origin for all remains, hence the discovered species most likely lived in the same local ecosystem. Despite partial diagenetic overprinting, the delta(18)O(PO4) values of the fossils indicate sufficient taxon-specific isotopic diversity to permit some broad conclusions on the palaeo-environment of the fossils. In particular, it is apparent that the isotopic composition of the Pycnodontiformes fish and Mosasauridae remains is most compatible with a freshwater palaeo-habitat and incompatible with a marine palaeo-environment. In addition, the Sr concentration and isotope data indicate that the Pycnodontiformes and Mosasauridae likely lived predominantly in a freshwater environment and were not simply occasional visitors to the Iharkut river ecosystem. Regarding other fossil groups, high delta(18)O(PO4) values of Alligatoroidea and Iharkutosuchus teeth suggest that these small crocodile species might have inhabited swamps and ponds where the water was relatively rich in (18)O due to evaporation. (c) 2009 Elsevier B.V. All rights reserved.

Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution.

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.

Nd and Sr isotope compositions in modern and fossil bones - Proxies for vertebrate provenance and taphonomy

Rare earth elements (REE), while not essential for the physiologic functions of animals, are ingested and incorporated in ppb concentrations in bones and teeth. Nd isotope compositions of modern bones of animals from isotopically distinct habitats demonstrate that the (143)Nd/(144)Nd of the apatite can be used as a fingerprint for bedrock geology or ambient water mass. This potentially allows the provenance and migration of extant vertebrates to be traced, similar to the use of Sr isotopes. Although REE may be enriched by up to 5 orders of magnitude during diagenesis and recrystallization of bone apatite, in vivo (143)Nd/(144)Nd may be preserved in the inner cortex of fossil bones or enamel. However, tracking the provenance of ancient or extinct vertebrates is possible only for well-preserved archeological and paleontological skeletal remains with in vivo-like Nd contents at the ppb-level. Intra-bone and -tooth REE analysis can be used to screen for appropriate areas. Large intra-bone Nd concentration gradients of 10(1)-10(3) are often measured. Nd concentrations in the inner bone cortex increase over timescales of millions of years, while bone rims may be enriched over millenial timescales. Nevertheless, epsilon(Nd) values are often similar within one epsilon(Nd) unit within a single bone. Larger intra-bone differences in specimens may either reflect a partial preservation of in vivo values or changing epsilon(Nd) values of the diagenetic fluid during fossilization. However, most fossil specimens and the outer rims of bones will record taphonomic (143)Nd/(144)Nd incorporated post mortem during diagenesis. Unlike REE patterns, (143)Nd/(144)Nd are not biased by fractionation processes during REE-uptake into the apatite crystal lattice, hence the epsilon(Nd) value is an important tracer for taphonomy and reworking. Bones and teeth from autochthonous fossil assemblages have small variations of +/- 1 epsilon(Nd) unit only. In contrast, fossil bones and teeth from over 20 different marine and terrestrial fossil sites have a total range of epsilon(Nd) values from -13.0 to 4.9 (n = 80), often matching the composition of the embedding sediment. This implies that the surrounding sediment is the source of Nd in the fossil bones and that the specimens of this study seem not to have been reworked. Differences in epsilon(Nd) values between skeletal remains and embedding sediment may either indicate reworking of fossils and/or a REE-uptake from a diagenetic fluid with non-sediment derived epsilon(Nd) values. The latter often applies to fossil shark teeth, which may preserve paleo-seawater values. Complementary to epsilon(Nd) values, (87)Sr/(86)Sr can help to further constrain the fossil provenance and reworking. (C) 2011 Elsevier Ltd. All rights reserved.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Vertebrate and nematode genes coding for yolk proteins are derived from a common ancestor.

One of the most obvious characteristics of the egg cells of oviparous animals is their large size resulting to a major extent from the deposition of nutritional reserves, mainly constituted of yolk proteins. In general, these are derived from a precursor called vitellogenin, which undergoes posttranslational modifications during secretion and during transport into and storage within the oocytes. Comparative analysis of the structural organization of the vitellogenin gene and of its product in different species shows that the vitellogenin gene is very ancient and that in vertebrates the gene may have more resemblance to the earliest gene than in invertebrates.

Mites as biological tags of their hosts.

Movements and spatial distribution of host populations are expected to shape the genetic structure of their parasite populations. Comparing the genetic patterns of both interacting species may improve our understanding of their evolutionary history. Moreover, genetic analyses of parasites with horizontal transmission may serve as indicators of historical events or current demographic processes that are not apparent in the genetic signature of their hosts. Here, we compared mitochondrial variation in populations of the ectoparasitic mite Spinturnix myoti with the genetic pattern of its host, the Maghrebian bat Myotis punicus in North Africa and in the islands of Corsica and Sardinia. Mite mitochondrial differentiation among populations was correlated with both host mitochondrial and nuclear differentiation, suggesting spatial co-differentiation of the lineages of the two interacting species. Therefore our results suggest that parasite dispersal is exclusively mediated by host movements, with open water between landmasses as a main barrier for host and parasite dispersal. Surprisingly the unique presence of a continental European mite lineage in Corsica was inconsistent with host phylogeographical history and strongly suggests the former presence of European mouse-eared bats on this island. Parasites may thus act as biological tags to reveal the presence of their now locally extinct host.