Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

Fluid evolution at the Variscan front in the vicinity of the Aachen thrust

Quartz-carbonate-chlorite veins were studied in borehole samples of the RWTH-1 well in Aachen. Veins formed in Devonian rocks in the footwall of the Aachen thrust during Variscan deformation and associated fluid flow. Primary fluid inclusions indicate subsolvus unmixing of a homogenous H(2)O-CO(2)-CH(4)-(N(2))-Na-(K)-Cl fluid into a H(2)O-Na-(K)-Cl solution and a vapour-rich CO(2)-(H(2)O, CH(4), N(2)) fluid. The aqueous end-member composition resembles that of metamorphic fluids of the Variscan front zone with salinities ranging from 4 to 7% NaCl equiv. and maximum homogenisation temperatures of close to 400A degrees C. Pressure estimates indicate a burial depth between 4,500 and 8,000 m at geothermal gradients between 50 and 75A degrees C/26 MPa, but pressure decrease to sublithostatic conditions is also indicated, probably as a consequence of fracture opening during episodic seismic activity. A second fluid system, mainly preserved in pseudo-secondary and secondary fluid inclusions, is characterised by fluid temperatures between 200 and 250A degrees C and salinities of < 5% NaCl equiv. Bulk stable isotope analyses of fluids released from vein quartz, calcite, and dolomite by decrepitation yielded delta D(H2O) values from -89 to -113 aEuro degrees, delta(13)C(CH4) from -26.9 to -28.9aEuro degrees (VPDB) and delta(13)C(CO2) from -12.8 to -23.3aEuro degrees (VPDB). The low delta D and delta(13)C range of the fluids is considered to be due to interaction with cracked hydrocarbons. The second fluid influx caused partial isotope exchange and disequilibrium. It is envisaged that an initial short lived flux of hot metamorphic fluids expelled from the epizonal metamorphic domains of the Stavelot-Venn massif. The metamorphic fluid was focused along major thrust faults of the Variscan front zone such as the Aachen thrust. A second fluid influx was introduced from formation waters in the footwall of the Aachen thrust as a consequence of progressive deformation. Mixing of the cooler and lower salinity formation water with the hot metamorphic fluid during episodic fluid trapping resulted in an evolving range of physicochemical fluid inclusion characteristics.

Evolution of proteins after whole-genome duplication

Evolution of proteins after whole-genome duplicationGene and genome duplication are considered major mechanisms in the creation of newfunctions in genomes, or in the refinement of networks by the division of function amongmore genes. In animals, the best demonstrated whole genome duplication occurred at theorigin of Teleost fishes. This makes fishes an ideal model to study the consequences ofgenome duplication, particularly since we have a good sampling of genome sequences,abundant functional information, and a very well studied outgroup: the tetrapodes (includinghuman). More specifically, I studied the consequences of duplication on proteins usingevolutionary models to infer adaptive events. I analysed the influence of positive selection invertebrate genes, by contrasting singleton genes and duplicated genes. The conclusion of theanalyses was threefold: (i) positive selection affects diverse phylogenetic branches anddiverse gene categories during vertebrate evolution; (ii) it concerns only a small proportion ofsites (1%-5%); and (iii) whole genome duplication had no detectable impact on theprevalence of this positive selection.I also studied evolution at the amino acid level with different methods to detect functionalshifts (covarion process and constant-but-different process). As in my previous research, Ifound similar numbers of functional shifts between duplicates and between orthologs.The accepted framework for studies of molecular evolution is that orthologs share the samefunction, whereas the function of paralogs diverges. This framework gives a special place togene duplication in evolution, as the main mechanism for generating novelty. With myprevious results showing that duplication and speciation are not so different, we investigatedthe literature to question the evidence for similar or divergent evolution of gene function afterduplication relative to speciation genes. This led us to propose a more rigorous design offuture studies of gene duplication.Finally, based on my automated protocol, we built a database of positive selection invertebrates' genes, Selectome. This database is freely available on the web and will helpfuture evolutionary as well as biochemical studies.

Self-organized evolution in a socio-economic environment

We propose a general scenario to analyze technological changes in socio-economic environments. We illustrate the ideas with a model that incorporating the main trends is simple enough to extract analytical results and, at the same time, sufficiently complex to display a rich dynamic behavior. Our study shows that there exists a macroscopic observable that is maximized in a regime where the system is critical, in the sense that the distribution of events follow power laws. Computer simulations show that, in addition, the system always self-organizes to achieve the optimal performance in the stationary state.

Evolution under monogamy feminizes gene expression in Drosophila melanogaster.

Many genes have evolved sexually dimorphic expression as a consequence of divergent selection on males and females. However, because the sexes share a genome, the extent to which evolution can shape gene expression independently in each sex is controversial. Here, we use experimental evolution to reveal suboptimal sex-specific expression for much of the genome. By enforcing a monogamous mating system in populations of Drosophila melanogaster for over 100 generations, we eliminated major components of selection on males: female choice and male-male competition. If gene expression is subject to sexually antagonistic selection, relaxed selection on males should cause evolution towards female optima. Monogamous males and females show this pattern of feminization in both the whole-body and head transcriptomes. Genes with male-biased expression patterns evolved decreased expression under monogamy, while genes with female-biased expression evolved increased expression, relative to polygamous populations. Our results demonstrate persistent and widespread evolutionary tension between male and female adaptation.

Evolution of anaesthesia care and related events between 1996 and 2010 in Switzerland.

BACKGROUND: Anaesthesia Databank Switzerland (ADS) is a voluntary data registry introduced in 1996. Its ultimate goal is to promote quality in anaesthesiology. METHODS: The ADS registry analyses routinely recorded adverse events and provides benchmark comparisons between anaesthesia departments. Data collection comprises a set of 31 variables organised into three modules, one mandatory and two optional. RESULTS: In 2010, the database included 2,158,735 anaesthetic procedures. Over time, the proportions of older patients have increased, the largest group being aged 50-64 years. The percentage of patients with American Society of Anesthesiologists (ASA) status 1 has decreased while the percentage of ASA status 2 or 3 patients has increased. The most frequent comorbidities recorded were hypertension (21%), smoking (16%), allergy (15%) and obesity (12%). Between 1996 and 2010, 125,579 adverse events were recorded, of which 34% were cardiovascular, 7% respiratory, 39% technical and 20% non-specific. The most severe events were resuscitation (50%), oliguria (22%), myocardial ischaemia (17%) and haemorrhage (10%). CONCLUSION: Routine ADS data collection contributes to the monitoring of trends in anaesthesia care in Switzerland. The ADS system has proved to be usable in daily practice, although this remains a constant challenge that is highly dependent on local quality management and quality culture. Nevertheless, success in developing routine regular feedback to users to initiate discussions about anaesthetic events would most likely help strengthen departmental culture regarding safety and quality of care.

Magmatic-dominated fluid evolution in the Jurassic Nambija gold skarn deposits (southeastern Ecuador)

The Jurassic (approximately 145 Ma) Nambija oxidized gold skarns are hosted by the Triassic volcanosedimentary Piuntza unit in the sub-Andean zone of southeastern Ecuador. The skarns consist dominantly of granditic garnet (Ad(20-98)) with subordinate pyroxene (Di(46-92)Hd(17-42)Jo(0-19)) and epidote and are spatially associated with porphyritic quartz-diorite to granodiorite intrusions. Endoskarn is developed at the intrusion margins and grades inwards into a potassic alteration zone. Exoskarn has an outer K- and Na-enriched zone in the volcanosedimentary unit. Gold mineralization is associated with the weakly developed retrograde alteration of the exoskarn and occurs mainly in sulfide-poor vugs and milky quartz veins and veinlets in association with hematite. Fluid inclusion data for the main part of the prograde stage indicate the coexistence of high-temperature (500A degrees C to > 600A degrees C), high-salinity (up to 65 wt.% eq. NaCl), and moderate- to low-salinity aqueous-carbonic fluids interpreted to have been trapped at pressures around 100-120 MPa, corresponding to about 4-km depth. Lower-temperature (510-300A degrees C) and moderate- to low-salinity (23-2 wt.% eq. NaCl) aqueous fluids are recorded in garnet and epidote of the end of the prograde stage. The microthermometric data (Th from 513A degrees C to 318A degrees C and salinity from 1.0 to 23 wt.% eq. NaCl) and delta(18)O values between 6.2aEuro degrees and 11.5aEuro degrees for gold-bearing milky quartz from the retrograde stage suggest that the ore-forming fluid was dominantly magmatic. Pressures during the early retrograde stage were in the range of 50-100 MPa, in line with the evidence for CO(2) effervescence and probable local boiling. The dominance of magmatic low-saline to moderately saline oxidizing fluids during the retrograde stage is consistent with the depth of the skarn system, which could have delayed the ingression of external fluids until relatively low temperatures were reached. The resulting low water-to-rock ratios explain the weak retrograde alteration and the compositional variability of chlorite, essentially controlled by host rock compositions. Gold was precipitated at this stage as a result of cooling and pH increase related to CO(2) effervescence, which both result in destabilization of gold-bearing chloride complexes. Significant ingression of external fluids took place after gold deposition only, as recorded by delta(18)O values of 0.4aEuro degrees to 6.2aEuro degrees for fluids depositing quartz (below 350A degrees C) in sulfide-rich barren veins. Low-temperature (< 300A degrees C) meteoric fluids (delta(18)O(water) between -10.0aEuro degrees and -2.0aEuro degrees) are responsible for the precipitation of late comb quartz and calcite in cavities and veins and indicate mixing with cooler fluids of higher salinities (about 100A degrees C and 25 wt.% eq. NaCl). The latter are similar to low-temperature fluids (202-74.5A degrees C) with delta(18)O values of -0.5aEuro degrees to 3.1aEuro degrees and salinities in the range of 21.1 to 17.3 wt.% eq. CaCl(2), trapped in calcite of late veins and interpreted as basinal brines. Nambija represents a deep equivalent of the oxidized gold skarn class, the presence of CO(2) in the fluids being partly a consequence of the relatively deep setting at about 4-km depth. As in other Au-bearing skarn deposits, not only the prograde stage but also the gold-precipitating retrograde stage is dominated by fluids of magmatic origin.

Exact temporal evolution for some non-linear diffusion process

Exact solutions to FokkerPlanck equations with nonlinear drift are considered. Applications of these exact solutions for concrete models are studied. We arrive at the conclusion that for certain drifts we obtain divergent moments (and infinite relaxation time) if the diffusion process can be extended without any obstacle to the whole space. But if we introduce a potential barrier that limits the diffusion process, moments converge with a finite relaxation time.

Evolution of the Pelagonian carbonate platform complex and the adjacent oceanic realm in response to plate tectonic forcing (Late Triassic and Jurassic), Evvoia, Greece

The Late Triassic and Jurassic platform and the oceanic complexes in Evvoia, Greece, share a complementary plate-tectonic evolution. Shallow marine carbonate deposition responded to changing rates of subsidence and uplift, whilst the adjacent ocean underwent spreading, and then convergence, collision and finally obduction over the platform complex. Late Triassic ocean spreading correlated with platform subsidence and the formation of a long-persisting peritidal passive-margin platform. Incipient drowning occurred from the Sinemurian to the late Middle Jurassic. This subsidence correlated with intra-oceanic subduction and plate convergence that led to supra-subduction calc-alkaline magmatism and the formation of a primitive volcanic arc. During the Middle Jurassic, plate collision caused arc uplift above the carbonate compensation depth (CCD) in the oceanic realm, and related thrust-faulting, on the platform, led to sub-aerial exposures. Patch-reefs developed there during the Late Oxfordian to Kimmeridgian. Advanced oceanic nappe-loading caused platform drowning below the CCD during the Tithonian, which is documented by intercalations of reefal turbidites with non-carbonate radiolarites. Radiolarites and bypass-turbidites, consisting of siliciclastic greywacke, terminate the platform succession beneath the emplaced oceanic nappe during late Tithonian to Valanginian time.

RPS4Y gene family evolution in primates

Background: The RPS4 gene codifies for ribosomal protein S4, a very well-conserved protein present in all kingdoms. In primates, RPS4 is codified by two functional genes located on both sex chromosomes: the RPS4X and RPS4Y genes. In humans, RPS4Y is duplicated and the Y chromosome therefore carries a third functional paralog: RPS4Y2, which presents a testis-specific expression pattern. Results: DNA sequence analysis of the intronic and cDNA regions of RPS4Y genes from species covering the entire primate phylogeny showed that the duplication event leading to the second Y-linked copy occurred after the divergence of New World monkeys, about 35 million years ago. Maximum likelihood analyses of the synonymous and non-synonymous substitutions revealed that positive selection was acting on RPS4Y2 gene in the human lineage, which represents the first evidence of positive selection on a ribosomal protein gene. Putative positive amino acid replacements affected the three domains of the protein: one of these changes is located in the KOW protein domain and affects the unique invariable position of this motif, and might thus have a dramatic effect on the protein function.Conclusion: Here, we shed new light on the evolutionary history of RPS4Y gene family, especially on that of RPS4Y2. The results point that the RPS4Y1 gene might be maintained to compensate gene dosage between sexes, while RPS4Y2 might have acquired a new function, at least in the lineage leading to humans.

The non-LTR retrotransposons in Ciona intestinalis: new insights into the evolution of chordate genomes

Background: Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes. Results: Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families. Conclusions: The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.

Genesis and evolution of the Evx and Mox genes and the extended Hox and ParaHox gene clusters

Background: Hox and ParaHox gene clusters are thought to have resulted from the duplication of a ProtoHox gene cluster early in metazoan evolution. However, the origin and evolution of the other genes belonging to the extended Hox group of homeobox-containing genes, that is, Mox and Evx, remains obscure. We constructed phylogenetic trees with mouse, amphioxus and Drosophila extended Hox and other related Antennapedia-type homeobox gene sequences and analyzed the linkage data available for such genes.Results: We claim that neither Mox nor Evx is a Hox or ParaHox gene. We propose a scenariothat reconciles phylogeny with linkage data, in which an Evx/Mox ancestor gene linked to aProtoHox cluster was involved in a segmental tandem duplication event that generated an arrayof all Hox-like genes, referred to as the `coupled¿ cluster. A chromosomal breakage within thiscluster explains the current composition of the extended Hox cluster (with Evx, Hox and Moxgenes) and the ParaHox cluster.Conclusions: Most studies dealing with the origin and evolution of Hox and ParaHox clustershave not included the Hox-related genes Mox and Evx. Our phylogenetic analyses and theavailable linkage data in mammalian genomes support an evolutionary scenario in which anancestor of Evx and Mox was linked to the ProtoHox cluster, and that a tandem duplication of alarge genomic region early in metazoan evolution generated the Hox and ParaHox clusters, plusthe cluster-neighbors Evx and Mox. The large `coupled¿ Hox-like cluster EvxHox/MoxParaHox wassubsequently broken, thus grouping the Mox and Evx genes to the Hox clusters, and isolating theParaHox cluster.

Thermal evolution of gene expression profiles in Drosophila subobscura

Background: Despite its pervasiveness, the genetic basis of adaptation resulting in variation directly or indirectly related to temperature (climatic) gradients is poorly understood. By using 3-fold replicated laboratory thermal stocks covering much of the physiologically tolerable temperature range for the temperate (i.e., cold tolerant) species Drosophila subobscura we have assessed whole-genome transcriptional responses after three years of thermal adaptation, when the populations had already diverged for inversion frequencies, pre-adult life history components, and morphological traits. Total mRNA from each population was compared to a reference pool mRNA in a standard, highly replicated two-colour competitive hybridization experiment using cDNA microarrays.Results: A total of 306 (6.6%) cDNA clones were identified as 'differentially expressed' (following a false discovery rate correction) after contrasting the two furthest apart thermal selection regimes (i.e., 13°C vs . 22°C), also including four previously reported candidate genes for thermotolerance in Drosophila (Hsp26, Hsp68, Fst, and Treh). On the other hand, correlated patterns of gene expression were similar in cold- and warm-adapted populations. Analysis of functional categories defined by the Gene Ontology project point to an overrepresentation of genes involved in carbohydrate metabolism, nucleic acids metabolism and regulation of transcription among other categories. Although the location of differently expressed genes was approximately at random with respect to chromosomes, a physical mapping of 88 probes to the polytene chromosomes of D. subobscura has shown that a larger than expected number mapped inside inverted chromosomal segments.Conclusion: Our data suggest that a sizeable number of genes appear to be involved in thermal adaptation in Drosophila, with a substantial fraction implicated in metabolism. This apparently illustrates the formidable challenge to understanding the adaptive evolution of complex trait variation. Furthermore, some clustering of genes within inverted chromosomal sections was detected. Disentangling the effects of inversions will be obviously required in any future approach if we want to identify the relevant candidate genes.