"The 3/3 strategy": a successful multifaceted hospital wide hand hygiene intervention based on WHO and continuous quality improvement methodology.

BACKGROUND Only multifaceted hospital wide interventions have been successful in achieving sustained improvements in hand hygiene (HH) compliance. METHODOLOGY/PRINCIPAL FINDINGS Pre-post intervention study of HH performance at baseline (October 2007-December 2009) and during intervention, which included two phases. Phase 1 (2010) included multimodal WHO approach. Phase 2 (2011) added Continuous Quality Improvement (CQI) tools and was based on: a) Increase of alcohol hand rub (AHR) solution placement (from 0.57 dispensers/bed to 1.56); b) Increase in frequency of audits (three days every three weeks: "3/3 strategy"); c) Implementation of a standardized register form of HH corrective actions; d) Statistical Process Control (SPC) as time series analysis methodology through appropriate control charts. During the intervention period we performed 819 scheduled direct observation audits which provided data from 11,714 HH opportunities. The most remarkable findings were: a) significant improvements in HH compliance with respect to baseline (25% mean increase); b) sustained high level (82%) of HH compliance during intervention; c) significant increase in AHRs consumption over time; c) significant decrease in the rate of healthcare-acquired MRSA; d) small but significant improvements in HH compliance when comparing phase 2 to phase 1 [79.5% (95% CI: 78.2-80.7) vs 84.6% (95% CI:83.8-85.4), p<0.05]; e) successful use of control charts to identify significant negative and positive deviations (special causes) related to the HH compliance process over time ("positive": 90.1% as highest HH compliance coinciding with the "World hygiene day"; and "negative":73.7% as lowest HH compliance coinciding with a statutory lay-off proceeding). CONCLUSIONS/SIGNIFICANCE CQI tools may be a key addition to WHO strategy to maintain a good HH performance over time. In addition, SPC has shown to be a powerful methodology to detect special causes in HH performance (positive and negative) and to help establishing adequate feedback to healthcare workers.

Morphological, ecological and genetic aspects associated with endemism in the Fly Orchid group.

The European genus Ophrys (Orchidaceae) is famous for its insect-like floral morphology, an adaptation for a pseudocopulatory pollination strategy involving Hymenoptera males. A large number of endemic Ophrys species have recently been described, especially within the Mediterranean Basin, which is one of the major species diversity hotspots. Subtle morphological variation and specific pollinator dependence are the two main perceptible criteria for describing numerous endemic taxa. However, the degree to which endemics differ genetically remains a challenging question. Additionally, knowledge regarding the factors underlying the emergence of such endemic entities is limited. To achieve new insights regarding speciation processes in Ophrys, we have investigated species boundaries in the Fly Orchid group (Ophrys insectifera sensu lato) by examining morphological, ecological and genetic evidence. Classically, authors have recognized one widespread taxon (O. insectifera) and two endemics (O. aymoninii from France and O. subinsectifera from Spain). Our research has identified clear morphological and ecological factors segregating among these taxa; however, genetic differences were more ambiguous. Insights from cpDNA sequencing and amplified fragment length polymorphisms genotyping indicated a recent diversification in the three extant Fly Orchid species, which may have been further obscured by active migration and admixture across the European continent. Our genetic results still indicate weak but noticeable phylogeographic clustering that partially correlates with the described species. Particularly, we report several isolated haplotypes and genetic clusters in central and southeastern Europe. With regard to the morphological, ecological and genetic aspects, we discuss the endemism status within the Fly Orchid group from evolutionary, taxonomical and conservation perspectives.

Phenotypic and molecular characterization of Bruck syndrome (Osteogenesis imperfecta with contractures of the large joints) caused by a recessive mutation in PLOD2

Le Syndrome de Bruck (Bruck Syndrome; BS) est une maladie autosomique récessive assemblant la combinaison inhabituelle de fragilité osseuse semblable à celle de l'Ostéogenèse Imparfaite (0I) avec des contractures congénitales tendineuses et cutanées des grandes articulations («ptérygia»). Les cas décrits jusqu'à ce jour mettent en évidence une grande hétérogénéité du tableau clinique, liée en partie au manque d'un diagnostic biochimique ou moléculaire. Nous savons que dans le BS les gènes codant pour le collagène 1 ne sont pas mutés, mais savons néanmoins, grâce à l'étude du collagène extrait de biopsies osseuses, qu'il y a un déficit d'hydroxylation des résidus de lysine dans les télopeptides du collagène 1 qui servent à la formation des liens intermoléculaires (crosslinks) et donc à la stabilisation des fibres de collagène. Un locus génétique du BS à été mappé sur 17q12, mais le gène responsable sur ce locus reste inconnu; plus récemment, deux mutations dans le gène de la lysyl hydroxylase 2 (PLOD2, position chromosomique 3q23-q24) ont été identifiées, démontrant l'hétérogénéité génétique du ES. La proportion de ES liée à 17p22 (BS type 1) et celle liée à une mutation dans PLOD2 (BS type 2) est encore incertaine et nous manquons de données sur la corrélation phenotype-génotype. Nous avons étudié le cas d'un garçon avec des contractures et des ptérygia dès la naissance, combinées à une ostéopénie sévère de type OI menant à des fractures multiples. Ses urines contenaient une quantité élevée d'hydroxyproline, indiquant un remaniement important du tissu osseux, mais peu de produits de dégradation des crosslinks du collagène, indiquant donc une réduction de la proportion de crosslinks dans le collagène in vivo. Nous avons pu démontrer chez lui la présence d'une nouvelle mutation homozygote dans le gène PLOD2 menant à une substitution Arg598His; les deux parents du sujet étaient hétérozygotes pour la mutation et celle-ci était absente dans notre population témoin. La mutation est adjacente aux deux mutations rapportées précédemment (Gly601Val et Thr608Ile), ce qui suggère la présence d'un ''hotspot'' mutationnel mais aussi d'une région de grande importance fonctionnelle sur PLOD2 : cette observation est importante pour la création d'inhibiteurs de PLOD2, recherchés en ce moment pour le traitement de la fibrose. La combinaison de ptérygia et de fragilité osseuse, comme illustrée par notre patient est apparemment contradictoire et donc difficilement explicable mais indique que l'hydroxylation des résidus lysyl des télopeptides est importante non seulement pour la stabilité osseuse mais aussi dans la morphogénèse et la formation des articulations dans la période prénatale. Finalement, la mesure des produits de dégradation du collagène dans l'urine et l'analyse de mutation de PLOD2 permet le diagnostic du syndrome de Bruck et permet de le différencier de l'Osteogénèse Imparfaite. -- Bruck syndrome (BS) is a recessively-inherited phenotypic disorder featuring the unusual combination of skeletal changes resembling osteogenesis imperfecta (0I) with congenital contractures of the large joints. Clinical heterogeneity is apparent in cases reported thus far. While the genes coding for collagen 1 chains are unaffected in BS, there is biochemical evidence for a defect in the hydroxylation of lysine residues in collagen 1 telopeptides. One BS locus has been mapped at 17p12, but more recently, two mutations in the lysyl hydroxylase 2 gene (PLOD2, 3q23-q24) have been identified in BS, showing genetic heterogeneity. The proportion of BS cases linked to 17p22 (BS type 1) or caused by mutations in PLOD2 (BS type 2) is still uncertain, and phenotypic correlations are lacking. We report on a boy who had congenital contractures with pterygia at birth and severe 0I-like osteopenia and multiple frac-tures. His urine contained high amounts of hydroxyproline but low amounts of collagen crosslinks degradation products; and he was shown to be homozygous for a novel mutation leading to an Arg598His substitution in PLOD2. The mutation is adjacent to the two mutations previously reported (Gly601Val and Thr608Ile), suggesting a functionally important hotspot in PLOD2. The combination of pterygia with bone fragility, as illustrated by this case, is difficult to explain; it suggests that telopeptide lysyl hydroxylation must be involved in prenatal joint formation and morphogenesis. Collagen degradation products in urine and mutation analysis ofPLOD2 maybe used to diagnose BS and differentiate it from M.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Turnover and accumulation of genetic diversity across large time-scale cycles of isolation and connection of populations.

Major climatic and geological events but also population history (secondary contacts) have generated cycles of population isolation and connection of long and short periods. Recent empirical and theoretical studies suggest that fast evolutionary processes might be triggered by such events, as commonly illustrated in ecology by the adaptive radiation of cichlid fishes (isolation and reconnection of lakes and watersheds) and in epidemiology by the fast adaptation of the influenza virus (isolation and reconnection in hosts). We test whether cyclic population isolation and connection provide the raw material (standing genetic variation) for species evolution and diversification. Our analytical results demonstrate that population isolation and connection can provide, to populations, a high excess of genetic diversity compared with what is expected at equilibrium. This excess is either cyclic (high allele turnover) or cumulates with time depending on the duration of the isolation and the connection periods and the mutation rate. We show that diversification rates of animal clades are associated with specific periods of climatic cycles in the Quaternary. We finally discuss the importance of our results for macroevolutionary patterns and for the inference of population history from genomic data.

Behavioural Flexibility: evolutionary past and its role in a changing world

A través de la historia de la vida, gran parte de los organismos han desarrollado estrategias para responder a un mundo en constante cambio. Hoy en día, las actividades humanas producen cambios ambientales a una velocidad sin precedentes, lo cual se traduce en grandes desafíos para la persistencia de biodiversidad. Esta investigación evalúa las respuesta de los animales a los cambios ambientales enfocándose en la flexibilidad del comportamiento como estrategia adaptativa. En una primera aproximación a una escala evolutiva, se otorgan evidencias del vínculo hasta ahora tenue entre la cognición e historias de vida, entregando un claro apoyo a la relación entre longevidad, vida reproductiva y el tamaño del cerebro en mamíferos. La longevidad es el centro de muchas hipótesis respecto a las ventajas de desarrollar un cerebro grande, como por ejemplo en la hipótesis del buffer cognitivo y las respuestas flexibles frente a nuevos ambientes. En un segundo nivel, se abordan factores extrínsecos e intrínsecos que podrían explicar las diferencias individuales en innovación, un componente clave en la flexibilidad del comportamiento. Por medio de una aproximación experimental, se evalúan potenciales escenarios que podrían conducir a consistentes diferencias individuales en uno de los principales factores subyacentes a la innovación (i.e. la motivación), y el potencial control endocrino sobre estos escenarios. Posteriormente, con el objetivo de evaluar la respuesta de los animales frente a los cambios ambientales actuales, se explora la respuesta de los animales frente a una de las actividades humanas mas disruptivas sobre los ecosistemas, la urbanización. Por medio de un analisis filogenetico comparativo a nivel global en aves se abordan los mecanismos implicados en la perdida de biodiversidad observada en ambientes urbanos. Los resultados entregan evidencias sobre la importancia de procesos de dispersión local junto con el papel clave de los rasgos de historia de vida, pero en un sentido diferente al clasicamente pensado. Finalmente por medio de una revisión bibliográfica se entregan evidencias teóricas y empíricas que respaldan el rol clave de la flexibilidad del comportamiento en confrontar los desafíos de una vida urbana. La integración de estos resultados muestra cómo el pasado evolutivo contribuye a hacer frente a los retos ambientales actuales, y pone de relieve posibles consecuencias ante un planeta más cambiante que nunca.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completelyabsent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and byMartín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involvedparts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method isintroduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that thetheoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approachhas reasonable properties from a compositional point of view. In particular, it is “natural” in the sense thatit recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in thesame paper a substitution method for missing values on compositional data sets is introduced

Behavioural Flexibility: evolutionary past and its role in a changing world

A través de la historia de la vida, gran parte de los organismos han desarrollado estrategias para responder a un mundo en constante cambio. Hoy en día, las actividades humanas producen cambios ambientales a una velocidad sin precedentes, lo cual se traduce en grandes desafíos para la persistencia de biodiversidad. Esta investigación evalúa las respuesta de los animales a los cambios ambientales enfocándose en la flexibilidad del comportamiento como estrategia adaptativa. En una primera aproximación a una escala evolutiva, se otorgan evidencias del vínculo hasta ahora tenue entre la cognición e historias de vida, entregando un claro apoyo a la relación entre longevidad, vida reproductiva y el tamaño del cerebro en mamíferos. La longevidad es el centro de muchas hipótesis respecto a las ventajas de desarrollar un cerebro grande, como por ejemplo en la hipótesis del buffer cognitivo y las respuestas flexibles frente a nuevos ambientes. En un segundo nivel, se abordan factores extrínsecos e intrínsecos que podrían explicar las diferencias individuales en innovación, un componente clave en la flexibilidad del comportamiento. Por medio de una aproximación experimental, se evalúan potenciales escenarios que podrían conducir a consistentes diferencias individuales en uno de los principales factores subyacentes a la innovación (i.e. la motivación), y el potencial control endocrino sobre estos escenarios. Posteriormente, con el objetivo de evaluar la respuesta de los animales frente a los cambios ambientales actuales, se explora la respuesta de los animales frente a una de las actividades humanas mas disruptivas sobre los ecosistemas, la urbanización. Por medio de un analisis filogenetico comparativo a nivel global en aves se abordan los mecanismos implicados en la perdida de biodiversidad observada en ambientes urbanos. Los resultados entregan evidencias sobre la importancia de procesos de dispersión local junto con el papel clave de los rasgos de historia de vida, pero en un sentido diferente al clasicamente pensado. Finalmente por medio de una revisión bibliográfica se entregan evidencias teóricas y empíricas que respaldan el rol clave de la flexibilidad del comportamiento en confrontar los desafíos de una vida urbana. La integración de estos resultados muestra cómo el pasado evolutivo contribuye a hacer frente a los retos ambientales actuales, y pone de relieve posibles consecuencias ante un planeta más cambiante que nunca.

International analysis of the countries where Brazilian franchise chains operate

This paper aims to demonstrate which external environment factors are involved in the international commitment of Brazilian franchise chains. Our objectives herein are to understand which external country characteristics lead to international franchising operations and to ascertain the influence of such characteristics in the commitment of franchise chains in each country they operate in. The database has 54 Brazilian franchise chains with international operations in 26 countries, which implies in 157 franchises operating abroad. Regarding external environment factors, the independent variables form a group divided into market opportunity and business efficiency – trust and ease to doing business. The result of the three distinct clusters show how the market opportunity and the business efficiency (trust and ease in doing business) work as drivers to the international operation of Brazilian franchises. The paper shows that the franchise chains operating in the USA (cluster 3) have an inferior international commitment in comparison with the franchises which operate in developed countries and in small Latin American countries (cluster 2). It is also possible to notice a large number of franchises that operate in underdeveloped countries from Latin America and Africa (cluster 1) with worse business efficiency due to the advantage of learning how to operate in a country that could have some similarities with Brazil

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Stochastic stability and the evolution of coordination in spatially structured populations.

Animals can often coordinate their actions to achieve mutually beneficial outcomes. However, this can result in a social dilemma when uncertainty about the behavior of partners creates multiple fitness peaks. Strategies that minimize risk ("risk dominant") instead of maximizing reward ("payoff dominant") are favored in economic models when individuals learn behaviors that increase their payoffs. Specifically, such strategies are shown to be "stochastically stable" (a refinement of evolutionary stability). Here, we extend the notion of stochastic stability to biological models of continuous phenotypes at a mutation-selection-drift balance. This allows us to make a unique prediction for long-term evolution in games with multiple equilibria. We show how genetic relatedness due to limited dispersal and scaled to account for local competition can crucially affect the stochastically-stable outcome of coordination games. We find that positive relatedness (weak local competition) increases the chance the payoff dominant strategy is stochastically stable, even when it is not risk dominant. Conversely, negative relatedness (strong local competition) increases the chance that strategies evolve that are neither payoff nor risk dominant. Extending our results to large multiplayer coordination games we find that negative relatedness can create competition so extreme that the game effectively changes to a hawk-dove game and a stochastically stable polymorphism between the alternative strategies evolves. These results demonstrate the usefulness of stochastic stability in characterizing long-term evolution of continuous phenotypes: the outcomes of multiplayer games can be reduced to the generic equilibria of two-player games and the effect of spatial structure can be analyzed readily.

Untangling the evolution of Rab G proteins: implications of a comprehensive genomic analysis.

BACKGROUND: Membrane-bound organelles are a defining feature of eukaryotic cells, and play a central role in most of their fundamental processes. The Rab G proteins are the single largest family of proteins that participate in the traffic between organelles, with 66 Rabs encoded in the human genome. Rabs direct the organelle-specific recruitment of vesicle tethering factors, motor proteins, and regulators of membrane traffic. Each organelle or vesicle class is typically associated with one or more Rab, with the Rabs present in a particular cell reflecting that cell's complement of organelles and trafficking routes. RESULTS: Through iterative use of hidden Markov models and tree building, we classified Rabs across the eukaryotic kingdom to provide the most comprehensive view of Rab evolution obtained to date. A strikingly large repertoire of at least 20 Rabs appears to have been present in the last eukaryotic common ancestor (LECA), consistent with the 'complexity early' view of eukaryotic evolution. We were able to place these Rabs into six supergroups, giving a deep view into eukaryotic prehistory. CONCLUSIONS: Tracing the fate of the LECA Rabs revealed extensive losses with many extant eukaryotes having fewer Rabs, and none having the full complement. We found that other Rabs have expanded and diversified, including a large expansion at the dawn of metazoans, which could be followed to provide an account of the evolutionary history of all human Rabs. Some Rab changes could be correlated with differences in cellular organization, and the relative lack of variation in other families of membrane-traffic proteins suggests that it is the changes in Rabs that primarily underlies the variation in organelles between species and cell types.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.